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ON THE STABILITY OF A RELAY SYSTEM'S EQUILIBRIUM STATE 


S. D. Kinyapin and Yu. 1. Neimark 
(Gor’kii) 


A number of works have been devoted to the investigation of the stability of a relay system's equilibrium 
state. Initially, stability conditions obtained on the basis of inductive reasoning and the results of M. V. Meerov's 
investigations were cited by Ya. Z. Tsypkin [1]. Subsequently, they received a basis in the work of L. S. Pontrya- 
gin and V. G. Boltyanskii, a brief discussion of which is contained in [2]* and, in the case when the difference 
of the orders of the denominator and numerator of the relay system's linear links transfer constant is two, in 
work [3]. 


In the present paper a new derivation of the stability conditions for a relay system's equilibrium state is 
given. With this, the approach to the question is borrowed from work [3] but the investigation of the transformation's 
fixed points’ stability is carried out in a different way, being based on the relationship, established in [4], between 
the stability of the equilibrium state of the system of differential equations and the stability of the fixed points of 
the corresponding mapping. 


Preliminary Remarks and Posing of the Problem; Information About the Point Mapping 


If the transfer constant, K(p), or the admittance function, ¢(t), of a relay system's linear link, shown on Fig. 
1 with the relay characteristic shown in Fig. 2, satisfy the conditions 





lim K(p) = (0) =0, lim pK (p) = 9’ (0) <0, (1) 
p+ poo 


then, as is well known [5], the so-called sliding motion is possible in the relay system. 
If conditions (1) hold,then, on the relay switching plane in the relay 




















system's phase space, there is a so-called sliding motion lamina. The equili- 
ar brium state lies on the sliding motion lamina, and is stable if all the zeroes 
. of K(p) lie in the left hand-plane, and unstable otherwise [5]. The phase tra- 


z y jectories close to, and on, the lamina are represented in Fig. 3. 











As ¢'(0) goes to zero, the sliding motion lamina disappears, and sliding 
mode operation can no longer occur. Convergence of the phase points to the 
equilibrium state is only possible for ever accelerating relay switching, as 

Fig. 1. shown on Fig. 4. More exactly, in order that the equilibrium state be asympto- 
tically stable, it is necessary that, first of all, for the phase points which are 
approaching the equilibrium position the time interval between successive relay switchings decrease without bound 
and, secondly, the points of intersection of the phase trajectory with the relay switching plame, My, M;, Mg, ..., 
converge to the equilibrium position for those points which are in the neighborhood of the equilibrium position. 


Relay 

















* After this paper was sent to the printer, work [8] was published, introducing simplifications in the investigations 
of L. S. Pontryagin and V. G. Boltyanskii. 
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We assume that 








y 
e_—_—_—_— n 
Ay, (p) c 
at = 2 
K (p) B, 2 pti,’ (2) 
Tt 
where 
—T! 
Fig. 2. Am (p) = a) p™ + a, p™ "+... + 4m, a) 0, (3) 
By (p) = p® + yp? ++... + dn, 
and Ay, Ag, .. ., Ap are roots of the polynomial B,(p). Since we 


are interested in the case when ¢(0) = ¢ (0), then m= n~2 and 
Cy +eg+...+¢,=0. (4) 


The equation of motion of the relay system under considera- 
tion may be written in the form® [5, 6] 





H+hix = csign(4,+a%y+...+2n). (5) 





We take x1, %, . . ., X, a8 the phase variables. In the region 
G,, where x, + x+...+X,> 0, andinregionG , where x, + x + 
+...+ X, < 0, the solution of equations (5) has the form: 














s,m af gM Oe 4 % 4 — ei tH) 6 
i= He a, "¢—-¢ ), (6) 
s 
% where x} = x;(t)) and the upper sign should be chosen for region Gy 
the lower sign for region G_. 
Mp Let M(x, %, . . -, Xq) be a point lying on the relay switching 
E,! plane 

rz=%+2,+...+2,=0. (7) 


Fig. 4. 


A phase trajectory passing through point M is directed towards region G, if x > 0, and towards region G_ if 
%< 0. By adding equations (5) termwise, and keeping (4) in mind, we find that 


Bm By... +n = — yy — dg Zg—.-s —AnZn: (8) 
Let M(x;, . . ., Xn) and M(x, . . ., %y) be two successive points of intersection of phase trajectory L with the 


relay switching plane x = 0. Segment MM of the phase trajectory will lie in Gy if 4yx,+...+ ALx_ < 0, and 
in G_ for the inverse inequality. In accordance with (6), there will hold the following relationship between the 
coordinates of points M and M 


Z = ne7 4 + reat —e"), (9) 


* By definition, sign g 7 +1 for E > 0, 


~1 for — < 0. 
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the 





where T, the time of motion of the phase point from M to M, is the least positive root of the equation 


aiuaiid c ia 
diet + Tie }=0., (10) 


In (9) and (10), the upper or lower sign is chosen as a function of the sign of the quantity Ayx,+ . . . + Apxp. 
Equation (10) can be written in the form 


 (t) = 9(t), 
where 
—Ayt dé -w 
O(c) = F D) aie +" and (Fe) eny = F WH <O. 
In order that, in a neighborhood close to the equilibrium position, x, =. . . = x, = 0, equation (10) have 


positive roots which tend to zero as point M approaches the equilibrium point, the function g(t) must be negative 
for sufficiently small t. Since, for small t 


e(t)=—+t—+— tt , 
where 
d, =(— 1)" Sear, (11) 
the first nonzero coefficient d, must be positive, i.e., 
if d, = d, =... > dy—; = 0, then d, > 0. (12) 


Relationships (12) provide a necessary condition for the stability of the equilibrium state in the relay's 
switching mode. If conditions (12) hold,then, for stability in the small of the equilibrium state M, it is necessary 
and sufficient that the sequence of points, M, M, M, . . ., defined by relationship (9), tend to equilibrium state 
M® when point M is sufficiently close to M®. 


Relationship (9) may be considered as a point mapping of the switching plane into itself, translating point 
MtoM. The point M*, with coordinates x,, = x, =. . . = X_ = 0, is mapped into itself by this transformation. 
Thus, the question of the stability of equilibrium point M® of the relay system reduces to the question of the 
stability of the fixed point M® of mapping (9). 


Mapping (9) maps every point for which one of the following inequalities holds 
hy Ty + Ag det. .- + AnTn <0, yt + Agee +... Anta >, (13) 


into apoint for which the inverse inequality holds. Therefore, if relay switching plane x = 0 is subdivided into 
two half-planes, E, and E_, in which the first or, respectively, the second of the inequalities in (13) holds, then 
mapping (9) will map half-plane E4 on E_, and half-plane E_ on E,. If we now introduce the new coordinates 
X} =—x, in half-plane E_, then the transformations of E, into E_ and of E_ into E+ will have the identical form. 
Specifically, 


‘ all At | be —Ayt sini’ , ~Agts ba | ho —hy tT 
= — xe » vad et "and aj me x, e ), 


where T, and T, are the least positive roots of the equations, respectively , 2x} = 0 and Dx; = 0. 
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Therefore, the question of the stability of the fixed point M®* of transformation (9) reduces to the investi- 
gation of stability of the fixed point of the transformation 





m= — me t* — =~ (4 —e**), (14) 
i 
where T is the least positive root of the equation 
—hyt o% —Agt, 
rye i—e = 0), 15 
> i + x, ( ) ) (15) 
where, for the variables xj, the following inequality holds 
dri <0. (16) 
Investigation of the Stability of the Fixed Point 
By expanding the right members of (14) and equation (15) in series in terms of T, we find that 
in A,t 
A= —m(L—uet...J—ar(t—t+...), (17) 
A. t 
Yat—uet..)tee(1— +...) =0. (18) 
For the sequel, it is convenient to go to the new variables 
uy = >) Xi. (19) 


If we multiply the first equation of (17) by rd, the second by rd, etc. and then add termwise, we find that 
By = — uy + (— 1) ed —ii td 
U; Uy + (— 1) tds + tu jys + (— 1) > djit---, (20) 


where the repeated dots indicate terms of higher than second order in uy, ug, . . ., Uy, —y and T. In the new 
variables, the equation for T takes the form: 


—y—Fa+put...=0. (21) 
Further investigation of the stability of fixed point M®, with the coordinates wy = ug =... Up _y = 0, will 
be carried out separately for the cases nm = 2 and n-m 2 3. 
1. For nm = 2 and d; > 0, we find from (21) that 
2u; 4d; 


ie Wa tee | woes (22) 





After substitution of (22) in (20), we go to the following formulae for the transformation T under consider- 
ation by us; 
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) 


5) 


7) 


18) 


19) 


hat 


20) 


21) 


will 


22) 


der- 








Uy = uj —=— ur+..., 
3 d? 1 
= d d u2 9 d 
y= — uy — 2m — 24 yu, — 2H 42 L(g, 2B) 4 


2dy_1 2d,_4 uu, 4 (23) 


thes =e Un—) aa (-- 1)" ai uy; + (— 1)" a uy ug-—— 2 dy 
1 


— 2u? 2 dgd,_, 
+ (— 1) a (h—F a ) 














It is easily seen that the general method [4] of investigating the stability of the fixed point M of transfor- 
mation (23) does not lead to the goal, since all roots of the characteristic equation have unit modulus. In this 
particular case, investigation of stability can be carried out in the following way. We set up the transformation 


u, =U So ten SS 24 
1 - . aoe. - ae. a one Kes 
3 d? 1 3 d@ 1 ’ 











=~ - d 4 @u? 
ee — ht etiam iets @ thao 
Pe ra : sie bt 0 0} * os tA ee ee (24) 
= _ 4d uu, 
Un—y = —Un_y +... = Una + (— 1)” = Wy Uy + 4—* + 
4d, d 
+ (19+. 


In accordance with §8 of work [4], the fixed point of transformation T* (and, consequently, of T also) is 
stable or unstable depending on whether the equilibrium state M® is stable or unstable in the linear approxi- 


mation of the system of differential equations 








du, 4 dy uy, 
we" tate 
4d? u 
dug d, uy Us Uy 241 uy 
at = a 2s +4 d, Tt + 3d8 t ’ 
es 6 a 0 - 82°61 * . . > . . . . . . . . . . . . . . . . . . (25) 
du, 4d Us u 4d, d “i uy 
a ia 7S ee +... 


- 


obtained from formula (24) for transformation T® by replacing by = , The stability of the equili- 


brium state of system (25) can be investigated by the ordinary means. Specifically, by noting that, according to 





(22) 


on we oe Pe 26 
p.. 2 (' a qtr} (26) 


and that, moreover, the following relationship holds (cf, Appendix 1) between the variables uy, ug, .. . U, 
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Un — by Un—) +...+ (-- 1)"-1 ba—1 uy +- (-- 1)" ba Uy = 0, (27) 


we write the characteristic equation of system (25) of differential equations in the form 











—Z +P 0 00... 00 
ad? 
7 H+d 10 0 0 
1 
— 1 dads d 
A (p) 3g 7 OS can 2.8 (28) 
Nd dsd,_, n@n—1 
(— 1)" a") di OD ox. Pp 1 
(—1)"-1b.. (—1)"4,4..... —b, 1 
By expanding determinant (28) (cf, Appendix 2), we find that 
id 
A (p) =(—3- +P) Am(p), (29) 


where A,,(p) is the numerator of the transfer constant of the relay system's linear link. In accordance with what 


has already been said, and in correspondence with (29), the equilibrium state for the relay system for n—m = 2 
will be stable if 


d,>0, —d,>0 (30) 


and if, in addition, all the roots of Am(p) lie in the left half-plane. By expressing d, and d, in terms of the coef- 
ficients of polynomials A,,(p) and Bp(p), we find that (30) is equivalent to the conditions 


A <0, a, —a,b, > 0. (31) 


2. We now turn to the consideration of the case when n—m = 3, Forn—m = 38, d; = 0 but dg #0. We 
introduce the new variable v, = u,/T and write equations (20) and (21) in the form 


Hy = — 0, + ty — dy — 5 g— dy + 2%, - 


ty = — uy + ty + ty + dy + (...), 


(32) 
Uyn—1) = —Un—1 + (— 1)"tdn_y + Un + (— 1)" dy + Pl. ...% 
—+3—Zd,—+(...)=0. (33) 


We find from equations (33) that 


By substituting the value of T just found in (32), we are led to the following formulae for the transformation 
T being considered by us: 
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ion 





“9 i 
¥, = 2, —zue+..., 


us = —6y,+u,+..., 











(34) 
x ae e 
Un—; = —Un—, + (— 1)" _ ( 1 3)+ 
The characteristic equation for transformation (34) is written in the form 
Seem mae 20 0 0 
’ 2 , ‘ 
—6, 2—z; 0, 0, ‘ns 0 
X@«| pats, ON Oe (35) 
il eae th ae +3, —i1-—z F 0 
Pak oe i 0, 0, »—i—s 
It has the (n —3)-fold root —1 and, in addition, the two roots defined by the equation 
2?—4z4+1=0, 
One of the roots of this quadratic equation is greater than unity in modulus and, consequently, the fixed 
point vy = ug =... = Uy _y = 0 of transformation (34) is unstable. From the instability of transformation (34) 


(since, with this, not only v, but also u, increases) flows the instability of the fixed point of transformation T and, 
consequently, the instability as well as of the equilibrium state of the relay system being considered. 


In the case when n—m = k, for any k = 3, we proceed analogously. Specifically, as a result of the replace- 


“8-1 and considering that d, ==... = = 0 and d, #0 
9 ng at d; d pe dy —4 k ’ 


- 





u Us 

ments = 47» V2 = a = 

we write formulae (20) and (21) in the form 
y= — y+ y— a +e + (k— yl ed 


- v (—1)*-» Ux td, 
y=—%+%—>+---+—EoH TEST: (36) 





a td 
De = — Pra $ Ve + (— 19 = + cosy 


Uk = —u,+(—1)* td +. ee 9 


Un—) = —Un—1 + (— 1)" tdi + - ey 


(—1)*—* uw, td, 


y— stage tee. + ki +agpit::- =? 


(37) 





By substituting in (36) the quantity 
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defined by (37), we write transformation T in the new variables in the form 


(— 1)'** », (k +1 —i) (—1)**u, 





¥, = ky, —S— D4 + 7 +...+-——“ t+ ---> 
* - (— 1)*+? », oi 
Vg = —(k + i)ko,+. oot ay UR + IR — EGE — 1) + 


(— 4)*+2 2u, 


+.--t+ gap te 





es — 4)*+—-1 
Pra = (— ir etotng 45 a “t i(k +41)k...3—i(i—1)... (38) 
? (k — 1) u, 
0+ {i—(k—2)}]+...—-—y— eee, 


aoe (— 1)* ++ », Sot 
uy, = (—1)*41(k 4+ 1)lyu+... 4+ —j Wk + tt —ti— 1)... 


oo (6—k+ 41) 4+... + haet..., 








k—1 (— 4)*+ 
k—z ees me ee TS 0 0 
k+i)k—2 — 1)*+2.2 
—(k+1yk &t% —:... GS 0 0 
X (zs) = ee nee + Sa ae eo CS & & #6 8 8 © Be 6.4) (O' F058! £6 ice oS @ Oo ed (39) 
_4)h2 
spe + pi et ks | 0... 0 
Rete ee oN. 4% wae ose a aan. @ | 
Te A ee ew ee os 5 of ate 6 bie 68 0, ..—i—g 











has the (n ~k + 1)-fold root ~1, with the remaining m roots defined by the equation 





k k—1 (— 1)*-2 (—1)*#" 

sg oie ee | =e |) a 
k+1)k—2 — 1)*+2.2 
—(k +1)k ( + ) i atin. @ -ti6 a ee Saar 


X*(s) = 
(2) — 


*eeveeeefeetgfefeskeeee#ee#e#see#ee##e#ee#e#e#*te#eee#e#se%e*# ee *# # #* 














(ay (kq sy a ykk—s 


which has a root greater than unity in absolute value, since the coefficient of z* 


any k = 3, there is instability. 


~! is greater than k. Thus, for 
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Appendix 1 








From the joint holding of the equations 


Mm +tat. - - + oy + Ty to = 0, 
@yAq + Seg + ey +H Eghy + Fq4)% = 9, 
wrP+oart+. . Hah + Fq5 Hn = 0 


it follows that 


; 
44...4 te 
Ay Age » AQ 
Ag=| AP AR. . AR ws | =O. 
AP AT. + AR My 








By expanding A, by the last column, and by introducing the notation 














ae 4 
ef { Ay Ms An 
Ay As An oe 68 0 © ot’ @ uae 
,°-1 23-1 . oa e—1 
of fag PA A» |= As, 
+: & oer aaa note oe, wage 
aoe | 
a® oo sort 0 0.0 &.0. eam 
ee 


we have 
Ae = ugh d— ws A2 +. . . +(—1)"wARt = 0, 
We note that AD = 0, _ An —4 and that O, —. is the sum of all possible products of s elements of the n 


quantities Ay, Ag, ..., Ap. Therefore, by bearing in mind that 0, —, is the coefficient b, —,, , of polynomial 
B,(p), we finally obtain 


tied, + (— 1)aexdy_, +. (— 1), = 0. 


Appendix 2 





4 (p) = (p —-5-t) 0). 
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where 


(—1)""%,_, (—1™-%,_, (—1)"*b,_,..-—b 1 
dy 
pP+z ‘ 0 0 0 
d 
A! (p) = -7 p { Oe ee (1) 
i syne 0 0 p 1 
1 








By expanding the determinant Ap) by the first column, we find that it equals 


























(—1)"-*6,_, (—1)""*b,_,. - - hf 
P i 0 Oo 
d 
(—1)"-%,_,— (r+ 9) 0 Se 0 0 |— 
0 0 ie 
(— 1)"~%,_, (—1)"*—*,,_ —b, 1 
a i 0 ~~. © 
BJ 
Zz 0 p $s 22@ Oi—... (1) 
0 0 coca PD i 
(—1)°7%,_, (—'1)"7%,_ 4. by 4 
dys i 0 mee” @ 
*-. —a P i —s. wo oi= 
Ze Oe | Ep wcome ap 
dy d d 
= (—1)°-%,_,—(p + F) oo — Fa. . — Ay 
where 
(—4)**o,_, (—1)**,_ 4 - ht 
P | 0 O 
0 0 p i 
In accordance with [7] 
By = (— 1)" [ogg + Op ng H One aP t+ st bip F144 prt, (iil) 


After substitution of (III) in (11) and dropping of the nonessential factor (~1)"~* dj *, we get 


Ae =dyp** +... +p (did, ppt dedp pig t+ - +4, ilt--- 


sa + PUdiby gt dady gt «+ +Gg ight Gide gt Op gat... +dy_y- av) 


By expressing the coefficients bj and dj in terms of cy, Cg, . . ., C, amd Ay, Ag, . . ., Ap and by using 
expression (4) in expressing the c;, we get that polynomial A(p) coincides with polynomial Ap)(p), written in 
the form 
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iil) 


(iv) 


in 





An, (P) = 1 (p+ >a) (P+ As) « - » (P+ AQ) + (P+ Ar) Ca(P+ As) - - « (P+A,) 4 
&® . beds 1 ee PRES @a + (p+ 1)(p+A2)- - -(P+AQ_ 4) ey: 
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ON THE ANALYSIS OF THE STABILITY OF THE PERIODIC 
MODES OF OPERATION IN NONLINEAR CONTROL SYSTEMS 
WITH MANY DEGREES OF FREEDOM 


V. A. Taft 
(Moscow) 


On the basis of a generalized corollary of Hill's equation, a deriva- 
tion is given of the characteristic equation (in finite form) of a system with 
many degrees of freedom whose parameters are periodic functions of time. 
The results obtained permit the use of the well-known Mikhailov criterion 
for analyzing the stability of the periodic modes of operation. 


The question as to the stability of periodic modes of operation is quite essential for various technical 
problems. For example, in electrical engineering, an analysis of the existence of periodic modes turns out to be 
necessary for the investigation of the modes of operation of assemblies transmitting dc and ac energy (the stability 
of operation of rectifying and inverting units, the generation of subharmonic oscillations, etc.), for the investi- 
gation of the parametric generation of ac currents, and elsewhere. In automatic control theory, the investi- 
gation of periodic modes is important for the analysis of operation of vibration regulators with account taken of 
nonlinearities and other factors, and in radio engineering for the analysis of oscillation generator operation. 


The investigation of the stability of periodic modes of operation leads to the investigation of the stability 
of solutions of differential equations with periodic coefficients. A large number of works are devoted to this 
question, but there is no exact solution for the general case. Relatively recently there were made attempts to 
apply, to the solution of this problem, well-known graphical methods based on frequency methods [1, 2]. 


In particular, the problem of obtaining an approximate solution based on the assumption of the presence of 
a linear element with the properties of an ideal filter in the system being investigated, posed by M. A. Aizerman 
[3]. was considered earlier [4]. However, such an approximate solution does not give the capability of directly 
obtaining an answer to the question of taking into account the actual, nonidealized, frequency characteristic of 
the system's linear portion. 


In the present paper we make the attempt to solve the problem exactly without having recourse to the 
assumption of the presence of an ideal filter. 


The solution of the problem leads to the transformation of the characteristic equation of a system of equations 
with periodically varying coefficients, written in the form of an infinite determinant to a finite form, by using 
a method analogous to that employed in deriving the characteristic equation for the Hill equation. 


Initially, we limit ourselves to the consideration of a system with many degrees of freedom with one non- 
linear parameter, and we dwell on the transformation to finite form of the characteristic equation for the equation 
with small deviations from periodic motions. After the characteristic equation has been transformed to finite 
form, we can apply the well-known frequency criteria for analyzing stability. The methodology described is 
then generalized to the case of several equations with periodic coefficients. 


The results obtained can be directly used in the solution of the enumerated problems. 
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Derivation of the Characteristic Equation 








We consider an N-loop circuit in which the inductance of the principal arm of the vy ‘th loop is a nonlinear 
function of the current in the v ‘th loop. In this circuit, let there exist a periodic mode with fundamental frequency 


Q. For the given system, the equation for small deviations from the periodic motion can be presented in the fol- 
lowing form. * 


N 
P [Lu (t) iv (t)] + >} Mus (p) is (t) = 0, (1a) 


s=1 


N 
> Me, (p) ig (t) = 0 (E=1,2,. he ,v—i,v+1, > Se » N), (1b) 


t=] 


where M,,, and Mg , are functions of the operator p = d/dt of the form ~ 


i 


Sav P 


M,, =L,p+Rw+ 





For the nonlinear inductance of the v'th loop in the periodic mode we can write*® [5] 


Ly (t) = 2199 + 22% cos Qt +... .+ 218" cosmQi+...= >) Lae™, (2) 


(0) _ 0 
where ui) 21 s ; 
We will seek a solution of (1) in the form 


co co 
ig(th=et Dy IP Fer oe DY HP e™ (e=1,2,...,%,...,N). (3) 
k=—00 k =—o 


If, in (1b), we transfer all terms containing i,, to the right side, we obtain a system of (N~1) equations for 
the currents ig + (for s & v). In system (1b) of equations, by expressing the currents ig (s # v) in terms of iy, 
we get 


i, = Zev (p) i, (4) 


A 
where Z,, = — » A4,=Ayt+-.-.-+ Amys, A is the determinant of the system of equations obtained 


from (1b) after transferring iy to the right side, Ag, is the minor of determinant A obtained by striking out the 
€ "th row and the sth column. 


By substituting the values of i, thus found in (1a), we get 


N 
plow (t) in(t) + >) (Mw + Mve Zw) iv(t) = 0 (5) 


*By using electrical engineering terminology, we can consider that the inductance or capacitance is found in the 
corresponding principal arm of the circuit. 


m 
** Here, as is usual, we assume that the sum > lLwl is a finite quantity. 
m=0 
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P (Lw (t) iv (t)) + Zvin i, (t) = 0, 


where Z y jn is the input impedance measured between the variable inductance terminals connected in the 
principal arm of the v ‘th loop. 


If, in (5), we substitute the values of L,, ,,(t) and iy(t) from (2) and (3), we get 





oc n : ; 
D {0 + prey [ie + Sin CHAO) jm goto 


kee— 00 


oo co 
+ SY DD RAF (E+ m) Qy LOT ed+i +m at — 0. 
k=—oo oom 


After replacing (k + m) by k in the second sum, we get 





2. {a + jkQ) [x 4 = ye) (+580) ‘ 4 


o © (6) 
+ JD QB A+ Fk) LEIE-™ cosines — 0, 


hk =—o M=—oo 
m+0 


If the coefficients of like powers of e are equated to zero, we obtain the system of equations 


(184. Seem + 5 Li) je) _ 9), 


m=—co 





(7) 
(k=—oo,...,—1,0,41,...., co). 


We now introduce the notation*® 


A+ jkQ F (K+ 7k Q) * 


vain 





By substituting (8) in (7) we get 


(LY Fre in (+ fk Q) + Py inQh + FRAT + 
+ Fain(h}+ jkQ) YY LP” =0 


mM=— oo 








Zr Feain (A + 18 0) e) , Fin(A+IRQ) Ss ony pce—m) 
[1+ vi in (A + 7k Q) jis + Fo, + ikD) > 1Y-™ = 0. (9) 


* Here, and below, we let ag FO oo for to 14.) ~@. 
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7) 


8) 


(9) 





We let (A + jkQ) = My, and take oy, Oy, ..., a, as the roots of the polynomial Fy, in (* + jkQ) = Fyys inf iy). 


Then, (9) can be rewritten in the form 


+ Be ma) *% (u, — a,)* gs ry (p, — a) +..A08 + 
Att) nr git) n—r—t gi) (10) 
a Lae a en (k—m) 
+ 3 Ge a) +2 @,—apt + >» aay t fs zs = 0, 
= Fad | 


where n is the number of different roots of polynomial Fy; in(,), r is the number of simple roots, t is the number 
of double roots, etc. 


We now write the determinant of system (10), equal to A( A). 
By substituting yy, = A + jkQ, we transform the determinant A(A). 


Without computing the determinant A(A), we limit ourselves to writing, in general form, its diagonal and 
nondiagonal elements. For this we shall initially assume, for simplicity though without destroying the generality 


of the reasoning, that all the poles are either simple or double. A diagonal term of this determinant will have 
the following form: 


A n—r Aw 
+ Sarnia + Dope ayi (11) 
=] 








and a nondiagonal term will equal 





n At), n—r Aw s 
p K+ jk Q—a, v 2 Xtina (12) 
=) 


Thus, the difference between the determinant obtained here and the corresponding determinant for the Hill 
equation is defined by the number of poles and their multiplicity which, in the general case, here, can have 
arbitrary values, and by the presence of the term a in the denominators of the determinant's elements. 


The determinant A(A) is a meromorphic function with poles at the points (a, —jkQ; a, — jkQ,; 


tn —JkQ;...). If we expand this determinant in a Laurent series in the neighborhood of the point a, ~jkf 
we obtain 


A (A) = 4104) + 430) = 4,0) + 2 AP + & AP, (13) 
iQ {amp 


where A,( A) is the integral part and A,(2) is the fractional part of the Laurent series expansion which, in the case 
of a simple pole, equals 


co 
AP a= > CWA —a,+ fkQ)* G@—1,2,...,2) (14) 
k==—0o 
and, in the case of a double pole, equals 


APO) = y) CRA—a+jkQl* (G=4,2,....8—P). (15) 
k=—oo 
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The coefficient cl) equals (in the case of a simple pole® 
Cie = lim [A(A) (& — a4 + 7 Q)) (16) 
d-vi(ay — Jk Q) 
and for a multiple pole of multiplicity two 


Ci? = lim [A (A) (A — oy + fk Q)%). (17) 
) 


A+ (ai—jka 
We take into account that, for all values of k, Cj, = Dj, where D, is the determinant analogous to the Hill 
determinant for the case of a simple pole, equal to : 


The elements of the center row of this determinant and of the row shifted to k rows below the center row 
will have, respectively, the forms 








4 A; i A;, 0 . A; =i. 
*?% 1 ’ 1 1 , 
and 
Aj, 1 — Ai,o = i,—t 
E=1 a=1 E=1 


(for double poles, the denominators of the determinant elements must be raised to the second power, etc.). 


The sum of all the component fractional parts of the Laurent series of the function A()), corresponding to 
a given subscript i, for the case of a double pole equals [6] 





n—r n—r pe n? 
> Ay? (A) = >} ee TT" (19) 
i=l i=) Q* sin* [a —a,)x a 


Correspondingly, the component sum of the fractional parts of the function for the subscript i in the case of 
a simple pole will equal® * 


n n 
(1) Sham (1) a 8 i (i) 
Ber mi= Zhe = 2a [FO—a) | DP Fr (20) 


By using expressions (19-20) we can obtain the characteristic equation in finite form. 


* The coefficient Ce , corresponding to a simple fraction expansion, is found by well-known rules for multiple 
poles. For example, the coefficient cf) from a second-order pole will equal 


F as d h pe. . 
P= {a AM A—a+jRor} 


** By taking into account that d cot v/dv =~-sin* v, if there are poles of third and higher order in Z(p) we will 


find, in the expression for A(A), linear combinations of the functions * xt Fe (A — a,) | etc. 
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In correspondence with (11-14) we can write 


AQ) = O10) + 420) = 34n0) + FAP A+ F aWey. (21) 
im 


i=] i=] 


In correspondence with (19-20), as ReA->ao, A,(A) is a bounded quantity. 
We note that the function 4(A) has period jQ, in correspondence with which we limit ourselves to a consider- 
ation of the function in the band | ImA| <* . With this, lim A,(A)=0,alim A,(A) = lim A(A) 
ReA co ReA + co 


Re A-0o 
has a finite value equal to unity, i.e., 


>) Ax: (A) = 1. 


i=] 


In the particular case of the Hill equation, we obtain from (19) and (21) that 


sin? [= (A — a) a = rp? = . (22) 


After the characteristic equation has been put in finite form, one can use one of the well-known frequency 
criteria for stability analysis. It should be taken into account that the use of frequency criteria has its own peculiar 
features here. These peculiar features are determined by the fact that, in the given case, the characteristic function 
is a periodic function of frequency. In correspondence with this, it suffices here, as in the analogous case of pulse- 
amplitude modulation, to construct the hodograph of A(jw) as wvaries from — Q/2 to 9/2. 


We consider briefly the application of the Mikhailov frequency criterion to the system of equations of motion 
for which the characteristic equation has the form of (22). If a, is a pure imaginary number, then the Mikhailov 
hodograph of the left member of equation (22) will be a segment of a line coinciding with the real axis. 


Correspondingly, if a, is a real or complex number, the Mikhailov hodograph of the left member of (22) 
will be an infinite number of ellipses, superimposed one on another, and passing by turns through all quadrants. 
For small values of the left member of (22), i.e., for small values of pf), the Mikhailov hodograph will pass by 
turns through all quadrants, and the system will be stable. As the determinant pf?) increases, the right member of 
(22) increases, and the ellipses corresponding to the left member of (22) are shifted to the right and, for a certain 
limiting value of determinant pf), the system becomes unstable. The limiting value of determinant of) , cor- 
responding to the boundary of stability, is found from the condition [3] 


: , i 
Re [sin* (jo, — 2%) a = Di? =, (23) 
where py? depends on the modulation depth and w, is the most negative root of the equation 
Tm sin® [(jo — &;) a: 


On the basis of what has been presented here, we can solve the problem of contructing a control system 
which satisfies the conditions of stability for some periodic mode with fundamental frequency Q. As an example, 
we consider the control system described by the equations 


D(p)X =—K(p)¥, Y=F(X). (24) 


By means of a transformation analogous to those given above [3, 4], we arrive at the following infinite 
system of equations: 
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{140K (p) + D (py e**"') Xe + [ AamiK (pe!) Ky, + 
+ [x AomK (pe?) Kass te... + [5 Age K(p) OF] Kept mie 
+ [z Aym,K (p) gies Xngr +... =0, 


where Amy, is the coefficient of the kth harmonic in the Fourier series expansion of the function 





A, (t) = [Ar]. “e 


By setting up the determinant of system (25) and equating it to zero, we obtain the characteristic equation 
of the system, 


By dividing the k'th row of the determinant thus obtained by D(p), we obtain the determinant with diagonal 
terms of the form 


K (A + jQ) 
1+ Aap, + jQ) 
and nondiagonal terms of the form 
4 K(A + 
3 Ae TE ay 
We make the analogous assumption that 


Ne —»0 for Re(k)—> co. 


Let the multinomial D(p) have n different roots, of which r are simple and t double, and then the deter- 
minant of system (25) can be given in a form analogous to that of the determinant of system (10). 


By transforming the infinite determinant thus obtained to a finite form in correspondence with (14-20), just 
as this was done in the case considered above, we obtain the characteristic equation in the form 





n n—r (2) 9 

i D's 
scot | x(A— oy) —— | Df? — 3 : ~== i, 
2 ie = Osint| (Qa, x 


Here, the difference from the previous case is determined only by the denominators of the elements of 
determinant D. 


In conclusion, we consider the more general case when the variable parameters of L(t) occur in b loops, 
where b < n. In this case, instead of system (1) of equations, we can write 


b N 
P (Low (t) iv (t)] + SY) Mve(p)is(t) + >) Ma (p)iz(t) =O (v=1,2,...,b) 
s=1 E=b+1 


: 26 
> Mes(pig(t} =9 (E=6,..4,2) (26) 


s=1 
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where N is the number of loops, v is the number of a loop containing the variable inductance, € is the number 
of a loop containing only invariable parameters. 


In contradistinctlon to the cases considered above, we here assume that the periodic function Ly ,(t) can be 
given in the form of a finite series, containing m, terms, The solution, as above, will be sought in the form of 
(3). 


If we express all the ig in terms of the i,,, and carry out all the transformations given above in the tran- 
sition from (1) to (6), we obtain a system of b equations in the form 





K+ikQ 
k=—0o 
foe) ™ : 
FSD RATE+ my QL? Btierman 4 
k=—CO M=—ny, 
m+0o (27) 
oo b 
+ SD A+ fRQ) May (i + FQ) TL’ OTH! Giga... by. 


k=—0CO 8=] 
8+yv 


On the basis of (27) we can, analogously to what was done previously, write 





Z, in(A + 7k Q)7 a ; ae Sica 
(K+ jkQ)| LO + FE EP + + QL 1 4 
> 1 
+ Si (A+ jkQ) My, (h + fRQ)IS =0 (28) 
rs 
(k = —oo,...,—1,0,4,...,c), (w=4,2,..., B). 


System (28) of equations, for each value, contains b equations in b + 2m, unknowns, In this system, the 
unknowns are the variables 
19°, LP, IES ck tr ot Pee oa 


’ 


By using a method, analogous to [4], of replacing k by (k + 1), (k~1), etc,, we increase the system of equa- 
tions to b + 2m, equations, after which we eliminate from them all variables, ), if** ™ etc., for all v except, 
for example, vy = 1. Asa result, we obtain expressions for if®), if** 1)... which are analogous to (7), after 
which the necessary result can be obtained by repeated carrying out of the above operation, 


The author wishes to express his gratitude to I. M. Smirnova for reading the manuscript and making a 
number of valuable remarks. 
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COMPUTING THE TRANSIENT RESPONSE IN LINEAR SYSTEMS 
BY THE METHOD OF REDUCING THE ORDER OF THE 
DIFFERENTIAL EQUATION 


A. V. Kalyaev 


(Taganrog) 


The paper considers the computation of the transient responses in 
linear systems with constant parameters by the method of reducing the 
order of the differential equation. Examples are given of the application 
of the method developed. 


The question of calculating the nonstationary-state modes in linear systems has been developed in great 
detail. There exist such powerful instruments for analyzing these processes as the operator method. Its dis- 
advantage is the necessity of determining the roots of the characteristic equation which, for high orders of the 
differential equation, is a quite difficult task which can only be solved in numerical form. In order to avoid this 
difficulty a significant number of methods were suggested for the approximate construction of the transient 
responses in linear systems which allow the computation of the transient responses to be carried out without seek- 
ing the roots of the characteristic equation [1-4]. The development of computers, both analog and digital, 
allowed one to analyze, quite rapidly and simply, the transient responses both in linear and in nonlinear dynamic 
systems. 


The mathematic models and digital machines of today forced the approximate methods for computing 
transient responses into the background. The value of approximate methods today is retained only in so “ar as 
they can be used for setting up the simplest algorithm for solving the problem which is necessary for programming. 


However, it might still be urged that the value of approximate methods of analysis is far from being ex- 
hausted. Indeed, neither mathematical machines nor the existing contemporary methods of approximate calcu~- 
lation permit the solution of the problem being considered in a general literal form. It is unnecessary to prove 
that the importance of obtaining such a solution is in no way diminished as a result of the vigorous development 
of computational technology. If one succeeds in analyzing, at least approximately, the transient responses in 
high-order dynamic systems in general literal form, then this will be an important step in the development of 
the theory of nonsteady-state modes. 


In this paper we consider the question of seeking an approximate solution of a linear differential equation 
with constant coefficients by reducing the order of the equation, which in many cases permits one to obtain the 
solution in the most general literal form. 


We assume that there is some dynamic system which, within certain limits of its behavior, can be described 
by a linear differential equation with constant coefficients 


k=n 


>) oxy = f(t). (1) 


k=0 





In the majority of cases one can set up a sufficiently accurate representation of the properties of the linear 
system to be analyzed if one allows to act on the latter, not a signal of arbitrary form, but one of the standard 
signals. In particular, for analysing linear systems, wide usage is made of a standard signal in the form of a unit 
function, f(t) = 1(t). In this case, equation (1) takes the following form: 


k=n 


>) cxy™ = 4 (2). (2) 


k=0 


If the right member of the linear differential equation is a constant, then a simple transformation of the 
(y, t) coordinate system, translating it along the y axis, will take equation (2) into an equation with a zero right 
member. In the sequel we shall assume that such a coordinate transformation has been carried out and that we 
are given the following equation to be analyzed 


k=n 
> cxy™ == 0 (3) 
k=0 
with given initial conditions (for t = 0): 
y) =y = [k=0, 1, 2,...,(n—4)]. (4) 


We limit ourselves to the investigation of stable systems, the transient responses in which play an important 
role in practical cases. 


If there are no multiple roots, the solution of equation (3) has, as is well-known, the form: 


k=n 
y(t}= >) a,e*™. (5) 


k=] 


The presence of multiple roots, for which the solution of equation (3) is written in the form 


y(t) = Dj Py,-r(t)e**', 
k 


where Py» _4{*) 1s a polynomial of degree (7, ~1) in the variable t, gives nothing new to the theory and will not 
enter into the discussion which follows. Thus, we shall consider below only the solution in cases where there are 
no repeated roots. 


It frequently turns out that some of the terms of expression (5) are either small in absolute magnitude or 
are rapidly damped. Under these conditions it becomes possible to lower the order of differential equation (3) 
and, thereby, to simplify its solution. 


In fact, it is possible to approximate the exact curve y(t) by some approximate function 


k=m 
Yap (t) = >) ane™*", (6) 
k=] 
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for which m < n, i.e., the number of terms in the right member of expression (6) is less than the number of terms 
in the right member of expression (5). On the other hand, expression (6) corresponds to some differential equation 


k=m 
>i bp = 0, (7) 
k=0 


whose order is lower than the order of the original equation (3). 


We assume for definiteness that the initial conditions which correspond to equation (7) are the same as in 
the case of equation (3) (for t = 0): 


pad [k= 0, 1, 2,..., (m—4)}. (8) 
The coefficients b, (k = 0,1, 2, ..., m) are unknown. If we succeed in determining them by some means 


or other, then we can seek, not the exact solution of equation (3), but the approximate solution, Yapt*) y(t), the 
exact solution of equation (7) which has a lower order than the initial equation (3). Thus, we set ourselves the 
problem of determining the coefficients by (k= 0,1, 2,..., m). Without loss of generality we may choose one 
of these coefficients arbitrarily. For simplicity, we set b_, = 1. The remaining coefficients must be so chosen 
that the approximate function Yap(t) differs as little as possible from the exact function y(t). We require, for 
this, that the error which arises as the result of substituting the exact function y(t) in the approximating equation 
(7) be minimal: 


k=m 
8(t)= >) dey (2). (9) 
k=O 


We now set up the integral of the squared error 


I= (e (t) dt = Ves by (o] at, (10) 


0 keep 
which has a minimum when the following conditions hold 


al 
%, =9 {i = 0,1,2..., (m—4)]- (11) 


By substituting the value of Is from (10) in (11) and differentiating, we obtain the following system of 
equations: 


k=m—1 © - 
be \ y® yd dt = — [yim yi de [i =0,4,2,...,(m—1)s (12) 
k=0 0 0 


which allows us to determine the coefficients b,. The coefficients obtained lead to a minimum of the error, 5(t), 
but do not provide a minimum of the error 


e(t) = y(t) — Yap (2). (13) 
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To decrease the error €(t) one can require, in addition, that the following condition hold 


oo co 

1, =\e(t)dt =\ ty) — yap ()) dt = 0, (14) 
0 0 

while simultaneously discarding one of equations (12). 
It follows from relationship (14) that 


( vap(t a = (y (t) dt. (15) 


eo 
The integral \ Yep (t)dt can be determined by using equation (7). Indeed, by integrating this equation 
0 


between the limits of 0 and o, we obtain 


k=m oo 
Do Vy (t) dt = 0, 
k=0 P 


from whence, taking the initial conditions in (8) into account, we find that 


k=m 


be yan (a! = 5 b,yt”. 


k=] 


Thus, by taking (15) into account, we can write 


k=m 


bo\ y(t)dt = S) b,yt. (16) 
0 


ka} 


As the result, for determining the coefficients by [k = 0, 1, 2,...,(m-~1)], we obtain the following system 
of linear algebraic equations: 


k=am—1 © = 
> by yy dt = — {ym y de [i= 0,4,2,....(m—2)], 
k=—0 ’ 0 0 (17) 
Y byt =— yo), 
k=z=0 


where we have used the notation 


ow 8 


1144 








@ co 
The constant quantities of the form \y (t)dt and [yy dt, which enter into this system may 
0 0 


be easily determined by using the initial differential equation (3) and the corresponding initial conditions (4). 
Indeed, if we multiply equation (3) termwise by the quantities ys J=0,1,2,...,(0 —1)] and integrate between 
the limits of 0 and op while taking initial conditions (4) and final conditions y “) = oft = 0, 1, 2,.. . (n—1)} into 
account, we obtain the system of algebraic equations which is necessary for computing integrals of the type 


\y yt) de: 


0 


k=n 


> ex\ y™yOdt=O0 fi —0,1,2,...,(mn—4)]. (18) 
0 


k=0 


This system of n linear algebraic equations contains [n(n + 3)}/2 definite integrals: 


& 
Yn -| 1) ee Veer | 


However, in actuality, only n integrals of the following type are unknown 


Yur = (ty ae [ek = 0,4, 2,...,("—4)]- (19) 
0 


The remaining [n(n + 1)]/2 integrals are either removed by integrating by parts, if the sum (k + i) is odd: 


> p=='/, (k—i—8) t=co 

\ yy de = [ S (—Aypyte y+] + 

(— 1) et— 
{ES 





1) ==00 
[y%s+t—0] | (k >i), 


t 
2 

t= 
or are reduced to integrals of the form of (19) if the sum (k + 1) is even: 


- p= '/;(h—f—2) 


\y~ yOdt = [ SY (—1)Ppyt-» yoo] + 
0 p=0 


i yer) \ yrs dt (k>i). 


We note that system (18) of equations obtained above is compatible, since its determinant, as is easily 
shown, reduces to the Hurwitz determinant for Equation (3) and, consequently, is different from zero. 


co 


The integral \y(t)de, which enters into system (17) of equations is simply computed by integrating 
9 


equation (4) between the limits of zero and infinity: 
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ea k=n 


\y(eae = >) * yt; (20) 


0 k=] 


Thus, by computing integrals of the types \y dt and \y” y dt by means of expressions (18) and 


0 0 


(20) and then, from the system (17) of linear algebraic equations, determining the coefficients by [k= 0,1, 2,..., 
(m —1)], we take the problem of analyzing the transient response in the system described by equation (3) and trans- 
form it to the simpler problem of determining its approximate solution by means of a differential equation of 
lower order, namely, equation (7). The size of the error engendered by this can be estimated by computing the 
integral 


In =\8*(t)at, 


which presents no difficulties in theory if the approximate solution, Yap(t), is determined. Obviously, for the very 
best approximating function Yapt*) to the exact function y(t), it is necessary that the following inequality hold 


= le (Qa< (y (d) dt, 
0 0 


where the integral \y (t)dt is determined from system (18) of equations. 
0 


The results of carrying out the calculations given here for concrete linear dynamic systems operating in the 
nonsteady~state modes support the correctness of the method presented of lowering the order of linear differential 
equations with constant coefficients. 


In particular, the transient responses in several actual 
automatic control systems were computed by means of the 
method herein developed (Appendix I and II). With this, 
the order of the original differential equation was reduced by 
a factor or two or three. As the results of these computations, 
\Z shown in Figs. 1 and 2, demonstrate, the approximate curve, 

‘ Yaptt). does not differ much from the corresponding exact 
curve. 











\ 

N It should be emphasized that all the transformations 
‘ carried out in decreasing the order of an equation are linear 
» ones and, therefore, may be implemented in general literal 

| h. form. With this, if the lowering is carried out down to the 
pian, 


























* second order then, in certain cases, it is completely pos- 

4 0 tsec sible that the problem of investigating the transient response 
will be solved in general literal form. As an example of 

Fig. 1. 1) is the exact computed curve; 2) is such an analysis of transient responses, Appendix III gives 

the approximate solution; y = 47.6 goat the computation of the transient response in a third-order 

sin 0.42 t. system. Naturally, the method presented is not universal 

and, in certain cases, may lead to significant errors, but 
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for the most widely-used automatic control systems it is completely effective and can give satisfactory results. 
In all cases the method given can be recommended for investigating the quality of the transient response. The 
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Fig. 2. 1) is the exact computed curve; 2) is the approximate 
solution, y = [8.78-6.67 &~*-°6t_3. 0 ¢-8-72t sin (38.5t + 0.755)] 
107°. 


basic quality indicators — the duration of the transient response, the size of the overshoot and the number of oscil- 
lations — may, by this method, be computed with sufficient accuracy. 


Appendix I 








As an example, we consider the computation of the transient response in the system “airplane with auto- 
pilot” [1]. Such a system can be described by a sixth-order linear differential equation: 


y® + 16.4y + 107.4y + 364.2y + 1146.5y” + 774.2y’ + 292.4y = 0 (21) 
We assume that the initial conditions are given (for t = 0): 
Yo = 0, Yo = 20, Yo = 0, y\?) = 0, y = 0, y® = 0. (22) 


The exact solution of equation (21), obtained in [1], looks like: 


y (t) =0.164e—7-8* cos (3.35¢ + 0.732) — 1.232e—* cos (3.744 — 1.312) + 


P (23) 
+ 59.6e—°-9?* cos (0.428¢ + 1.57), 


We now determine the approximate solution of differential equation (21) by the method of lowering its 
order. By using system (18) of equations, set up for original equation (21) and initial conditions (22), we determine 
the integrals: 


co co co 
Yu \ vat = 1076, Yu= Sw dt = 356, Yo= \ (oy? at = 250. (24) 
0 0 0 
Equation (20) allows us to compute the integral: 
co 
\y (t) dt = 78.5, (25) 
0 
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We approximate original equation (21) by the second-order linear differential equation 


y” + diy’ + doy = 0, (26) 


where, for t = 0, 
yo=0, y= 20. 


To compute the coefficients b, and by we use system (17) of equations, which will have the following form 
in the given case: 


co 


"\2 
ny, -oe bye —by | ydt = — yi, (27) 
0 
From this we have 
b,.356 = 200, —b,-78.5 = — 20 


b, = 0.562, by = 0.255. 
Consequently, the approximating equation looks like: 


y” + 0.562y’ + 0.255y = 0. (28) 
Finally, we determine its solution with the given initial conditions 


y = 47.6¢—°8"! sin 0.428, (29) 


The curves corresponding to the exact solution, (23), of initial equation (21) and to its approximate solution, 
(29), are given in Fig. 1. As follows from the figure, despite the lowering of the differential equation's order by a 
factor of three, the solution of the approximating differential equation (28) coincides quite satisfactorily with the 
accurate solution, (23), of initial differential equation (21). 


Appendix II 





We now consider still another example of the computation of the transient response by the method of lower- 
ing the order of the differential equation, this time for an electro-hydraulic servo system [1] described by the 
equation 


y® +4 103y + 2065y" + 149250y’ + 1081500y = 9.5, (30) 


where, for t = 0, 
Yo=0, yo=0, yy=3.510-%, yy = 9.5°10-*. 


In the steady-state, we obviously have 


9.5 
Ys = [0R1500 = 8-78-10-*. 


Consequently, for the transient component, y,(t), of the function sought, y(t), taking the relationship y(t) = 
= yp(t) + y,{t) into account, we obtain the differential equation: 


yi + 103y%) + 3065y;, + 149250y;, + 1081500y,, = 0 (31) 
with initial conditions (for t = 0) 
Yno =— 8.78-10-* y,,=0, yng = 3.510%, yn, = 9,5-10-*. 
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By using system (18) of equations and equation (20), we find the integrals: 


co co 
\ y,dt = —8,68-10-7, \ (y,,)*dt = 8.16- 10-2, 
0 
© = (32) 
\ 2dt = 3.9-10-1%, \ (y,,)2dt = 8,34-10-7, 
0 


By now setting up system (17) of linear algebraic equations and solving it, we determine the coefficients of 


the approximating equation. If, as the approximating equation, we choose a third-order linear differential equa- 
tion, then the computation just cited yields 


bp = 1.22-108, b, = 1.605-108, by = 19.5. 
Thus, the approximating equation will look like: 
y,, + 19.5y,, + 1605y,, + 12200y,, = 0; (33) 
for t = 0 


Yno = 8.78-10-*, yn4 =0, Ying = 3.5-10-*. 


Its solution has the following form: 


Yn = —6.67-10-%e— _ 3.09. 10-te-*™ sin (38.5¢ + 0.755). (34) 
Hence, the approximate expression for the function sought, y = y, + y, is: 


y = [8.78 — 6,67e—%™ — 3. 09¢—*"# sin (38.58 + 0.755)]-10-*. (35) 
The exact solution of initial equation (30) is given by [1] 


y = [8.78 —6.85e~** + 0.06585" — 2.53¢-*4” sin (38.5¢ + 0.913)]-10-*. (36) 


The curves for the exact and the approximate solutions are shown on Fig. 2. As is clear, the deviations be~ 
tween them are very significant. 


Appendix III 











In conclusion, we consider the lowering of the order of a linear differential equation in general form. We 
assume that we are given the equation 
y” + esy” + ey’ + coy = 0 (37) 
and the initial conditions (for t = 0) corresponding to it 
y=0, y=y oy’ =0. 
We will assume that the system is stable, i.e., that the following conditions hold 


Co > 0, c, > 0, eg >0, C103 > eo. (38) 


We now lower the order of equation (37). For this, by means of system (18) of equations and equation (20), 
we compute the integrals: 


) P co 
- C2Yo > sal 3 + Co ° 
c cb+e om 
, ee. na YY 
rae = AT Wy 
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Then, by using system (17) of equations, we find the coefficients of the second-order approximation equation; 


C1Cq — Co 


As the result, the approximating equation takes the form: 
eo, Me—% ,, 2 
: Fla +a)" tan (4) 
where, for t = 0, 
y=9, y=% 
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ON DESIGNING CORRECTING CIRCUITS FOR AUTOMATIC CONTROL 
SYSTEMS IN ACCORDANCE WITH THE MEAN SQUARE ERROR CRITERION® 


V.ti. Kukhtenko 
(Moscow) 


A method is presented for finding the optimal (in the sense 
of minimum mean square error) transfer function, the degree of 
whose numerator is a given number of units less than the degree of 
the denominator, and a discussion is given of the relationship of 
this question with the question of physical realizability of the cor- 
rection circuit. 


An example is provided. 


The problem of designing control systems which are optimal in some sense is divided into two parts. First, 
it is necessary to determine the transfer function which is optimal from the point of view of the criterion chosen. 
Second, a method should be developed for synthesizing the optimal transfer function obtained by means of avail- 
able actual elements. These two parts are different in their physical and mathematical content and are most 
frequently solved in isolation from one another. 


In the makeup of automatic control systems there frequently enter elements whose presence is determined 
by requirements of design execution, character of operation, form of the energy sources provided, questions of 
reliability, etc., i.e., by conditions which relegate to the background of the design the dynamic properties of 
these elements. The dynamic quality of these “invariable” elements may turn out to be incompatible with the 
dynamic qualities of the optimal transfer function, defined by some criterion, in the sense that, for implementing 
the system as a whole, there might be required a physically nonrealizable correcting device. 


In this paper we consider a method of obtaining an optimal transfer function in accordance with the criterion 
considered in [1] and essentially identical with the criterion of Zadeh and Ragazzini [2]. 


A control system as a whole, or its individual elements, are considered below to be physically realizable 
in those cases wherein the transfer functions corresponding to them have the form of rational fractions in which 
the degree of the numerator is less than, or equal to, the degree of the denominator. This limitation is imposed 
in addition to the ordinarily adopted condition of physical realizability given by the equality 


where k(t) is the impulsive response of the system or its elements. The holding of equality (1) alone does not 
bespeak the possibility of realizing the system or the element. For example, the impulsive response of the ideal 


differentiating link, k(t) = 5(t), satisfies (1), nevertheless, in view of the unavoidable coordinate limitations and 
the inertia of all real devices, such a link in its pure form cannot be implemented. 


* Presented at the All- Moscow Seminar on probabilistic methods of automatic control, held at the IAT AN SSSR 
on April 24, 1958. 
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We now consider the influence of the structure of the “invariable* 
portion of the system on the form of the transfer function of the closed 
system, bearing in mind that, as correcting links, we shall use physically 
realizable (in the sense given above) elements. 





Figures 1 and 2 show two of the most frequently appearing schemes 
for correcting automatic control systems; on Fig. 1, by means of a cor- 
recting link in series and, on Fig. 2, by means of a feedback circuit. 
The transfer function of the object to be controlled (the “invariable” 
portion of the system) is denoted by W,(s), the transfer function of the 
correcting circuit by W,(s): 





No(s) P;(s) 
Mo(s)’ Hy, (s) x 


No(s) = ans* + an, 8"! +... + aS + do, 
M,(s) = bms™ + bm is™! + ... + dys + Do, 
P,(s) = Cps? + Cp_, SP +... +048 + Cp, 
H,(s) = d,s" + d,,s""+ ...+ds+d, 
(nm, pSh), 


Wo(s) = W, (s) = 




















We denote by K(s) the transfer function of the closed-loop control system, and present it in the form 


L(s) 


K(s)=F5> 


where L{s) and R(s) are polynomials ins of degrees, respectively, 7 andr. 
For the system with the correcting link in series, 


L (s) My (8) p 
No(s) [R (s)—L(s)] ” 





W,(s) = =—l+m, h=n+r. 


In this case, the correction link is realizable when h—p = 0, or 
r—l>m—n. 
For the system with the correction link in parallel, 


Rs) Mo(s) 


Wi (8s) =F (s)  No(s)° (6) 





“In the left member of (6) there are two rational fractions, whereby the second allows an integral part to be 
isolated 


M;(s) 


M,(s)_ °m nn 
Wey =o he Ee” 


where the degree of M,(s) equals n. 


In order that the function W,(s) be realizable, it is necessary that it not contain an integral part, which 
must be cancelled by the corresponding multinomial which is separated from the fraction R(s)/L{s). 


























ible* 


ed This gives the following condition for physical realizability in the case of feedback: 
call 

y r—l=m—n. (5a) 
lemes Consequently, the limitation imposed on the form of the optimal transfer function in these two very common 
-or~ cases consists of this, that the degree of its numerator must be (m~—n) less than the degree of its denominator. 
J We now form the integra}: 
he 


T= \ (k(t) ]® dt. | 
0 { 

(2) 
By the definition of the impulsive response, 


 8(¢—+)8(t)dt = 8(5) for t>0 | 
(3) F 


and 


[8 (t)8(¢) dt = 00, 


If it be required that the integral T, be finite, the function k(t) will then not have in its makeup a pulse 
constituent, the degree of whose numerator is one less than the degree of its denominator. Since differentiation 
in the t domain corresponds to multiplication in the s domain, the finitude of the integrals T, and 


T, =\ (k (rat ! 
(4) ; ! 


will mean that the degree of the numerator of the function sK(s) is one less than the degree of its denominator; 
correspondingly, the degree of the numerator of the function K(s) is two less than the degree of its denominator. 


| 
(5) Analogous reasoning leads to the conclusion that, for the degree of the numerator of K(s) to be q less than | 
the degree of its denominator, it is necessary and sufficient that the following inequalities hold: | 


[k (t)]® dt < 00, | 


T= 
(6) | 
T. =-{ [k’ (t)]® dt < 00 | 
tt to be : (8) 
To = { [k*—” (t)}* dt < 00. | 
(7) 0 | 
For the circuits shown in Figs. 1 and 2, let n(t) be the noise, given as a stationary random function of time, 
‘ich and let it be known of the signal, g(t), that its Laplace transform, g*(s), is a rational fraction. We note that 


everything presented below is extended, by analogy, to the case when the signal contains a stationary random 
component, m(t). 


We shall call an automatic control system optimal if [1} 









1) the mean square error 










72 = M [22 (t)] ={ \ Ry(t —9) k(t) k (0) de d0 (9) 


(where R,(t) is the correlation function of the noise and k(t) is the impulsive response of the closed-loop system) 
is minimal; 






2) the instantaneous value of the error Component corresponding to the forced motions [3] engendered by 
the signal g(t) equals zero: 










e¢ (t) = g(t) — ag for (t) = 0. (10) 
















Condition (10) means that, after termination of the transient response, the signal, g(t), passes through the 
system without distortion. 


In accordance with [1], in order to satisfy (10) the impulsive response k(t) must satisfy the following equali- 





ties:* 


K (s, + Ai) = YR ee tae im 1 a 


0 


om ms ” —(8; +A; )t ee -. 
Os be +4; \ tk (i)e dt Eri, 
0 
mi—1 © A 
k x00] me (AyD aT (tye tO dt oon, 
as” - 8 = 8; +4; 


0 
(i <= 1,2,... 5%), 
> — 
Re [K (ax + iB, + Ax)] = 0) e hte)! cos Bal dt = 1 — kore. 
0 


Im [K (ox + iB, + Ax)] = —{k (t) e(*x+40! sin Bat dt. eoor, 


(11) 


| == (—14 ihe tk h(t) eat aw! y 
s= ap-+iBy +A, 


xX cos B, tdt = Seny—1) 1» 


0 





aes yn ( t™*—* h(t) eat andt y 
0 


Im i K 0] 
Os Joma tiBy ty 


x sin By tdt = em, 4.2% (ke =1,2,...,p, ¥+p—a). 


Here, Re denotes the real part, Im denotes the imaginary part, v is the number of real poles, sj, of the 
image g*(s), u is the number of complex poles Sx = oy + ig, of the image g*(s), m; and my, are the multi- 
plicities of the corresponding poles, A; and 4A, are arbitrary positive numbers, so chosen that the following con- 
dition holds 






Re ($j; + 4;) > 0. (12) 


*The equation corresponding to distortion indicator €3k was omitted in [1]. 
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(9) 


em) 


by ' 


(10) 





quali- 


(11) 


the 
iti- 
g con- 


(12) 





The quantities €¢ , the so-called distortion indicators, determine the accuracy of transmission of the 
controlled signal g(t), distinguished from g(t) by the fact that its image, §*(s), has poles s;+ Aj. The meaning 
of equations (11) consists in this, that they establish some still undefined measure of distortion of the fictitious 
signal Q(t) in its passage through the system K(s), which can then be so chosen that the basic signal, g(t), is re- 
produced in accordance with (10). For this, it is necessary that one set ay = €, pe = 0 in (11). 


As is obvious from the discussion above, the problem of determining the optimal impulsive response consists 
in minimizing the functional in (9), with additional limitations of the form of (8) and (11). Problems of this 


type were often considered in the theory of optimal systems (4, 5]. Investigation of the first variation of the 
functional 


n m—1 ql @ 
Kiko) =7—-23 Yo wK@)+ x Nn r\ [e” (e)]? de, (13) 
i=i j=0 0 


where K\(s,) conventionally denotes the right member of (11), leads to the following equation, of the Wiener-Hopf 
type, for determining k(t): 


n m—t 
| k(t) Ro(t—t)de—S) HY (— Ay eh ee tty 
i=1 j=0 (14) 
X [1,; 008 Bt — %,j8in B.t]=0 for ¢>0, 
q—1 (15 
Ra (t) = Rn (t) + Sd” (2). ) 
r=0 
Here, A, and y4j are Lagrange multipliers. 
The solution of equation (14) is carried out by the ordinary methods [5, 6]. 
Application of Parseval's formula to the second variation gives a sufficient condition for a minimum 
S,(w) >0, (16) 


where Sg(w) is the spectral density corresponding to the correlation function R,(t). 


The expression for the mean value of the square of the noise of the output provided by an optimal system, 
has the form: 


nm-—tl q-1 


=3 Bw O—-Be" nef ROA (a (17) 


i=1 j=0 


The Lagrange multipliers, yj, from the expression for the optimal transfer function are excluded by 
equations (11), in which it is necessary to set A; = foe 


In the general case, decomposition of the function §, (w) into two complex conjugate factors is impossible 
since, for this, one would have to solve an algebraic equation of high degree with unknown coefficients, which 
is impossible. Therefore, instead of A,, it is necessary to introduce new unknown coefficients, b,, by the formula 


P (w) = dgy (Ap) © + dgqi—y (Ar) OF F—D + 
«+ Gg (Ap) @® + Go (Ar) = Gyr (Ay) I] (b; + iw) (b, — iw), (18) 


where P({ w) is the numerator of the rational-fraction function S,(w). 


1155 

























After determining K(s) by formulae (4) aad (6), one can find the transfer function, W,(s), of the correcting 
circuit, which will depend on the undetermined parameters b, . The values of these parameters should now be so 
chosen that the coefficients of the integral part in the expression for W;(s) become equal to zero (with the exception 
of the coefficient of s°), Then, the transfer function W,(s) will have a numerator whose degree is not greater than 
that of the denominator, and, consequently, will satisfy the criterion of physical realizability adopted above. 


Example. Let the executive element in a control system be an electric motor with the transfer function 


K 1 
s ’ 


4 =. = . 
then AK a - volt 





The signal has the form: 
8 (t)= 80 + Bit, 





and the correlation function of the noise is 


R(t) = R,e*""!, Ro=1, a=! soles 


In accordance with the general theory for obtaining physically realizable correcting circuits, the degree of 
the numerator of the optimal transfer function must, in this case, be one less than the degree of the denominator, 
which leads to the equation 
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fs ~<a eat_\ R, (t—1) k(t) dt =0, 
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where Ro(t) = Roe I el + A, dt). 


If we use, for the solution, .the formula 





K (iw) = - i | B® (iw) | 
Y(iw) LY (— iw) |4+ 


and take into account that 













To tT S.(w) = do (b + iw) (6 — iw) 
io +A ~ (iw + A)?’ ’ (a+ iw)(a—iw) ’ 








B* (iw) = 





As*+- (Aa + FE) s+ Ea 
K (8) = S720 4b) s* +A (2b4-A) s+ A" 





The conditions for undistorted reproduction of the signal g(t) become transformed to the conditions for 
obtaining a system with second order astatism, which gives 
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The transfer function of the feedback circuit (Fig. 2) is 


A At 
be RA +) — > b+ Aye (20-4-4)4— F- dJs4-0% 
Ko ~X¥ “ar As*} (2b +- A) Os + A% 
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Wi (8) = 


The conditions of physical realizability have the form: 


A=K, 
A A? 
(24 + 6) — | (2b + A)] s* + [(2b + A) A—d] s+ A% 
As* — (2b + B) As + A*% 


a Ka— 
“A Je—b ° 


W;(s) = 





b 





For choosing the value of A it is necessary to take three factors into consideration: the duration of the 
transient response, the satisfying of inequality (16) and the magnitude of x, [1]. 


Computations provide the formula for the mean-square error 
2 = 4°0*B + A*al — 2aACD — 2ACa* +- Db*a + B*a* — C* +- D*aC + Da*C +- DBC 
2 (aB + C + Da* + a*) (DB —C) 
B=A(2b+A), C=A%, D=2A+5. 





x(t) 








Ls me ] By setting A = 3, we find that b = 2, From the expres- 
3 wecyptit, 2 ees sion for the optimal transfer function it is obvious that the 
duration of the transient response in this case is of the order 
of one or two seconds. 
A» > 0, which leads to the inequality b > a which, for a = 1 
: ; and b = 2, is satisfied. 
2 
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al A sufficient condition for a minimum has the form: 
\ The final expressions for the complete transfer function 
T 





and for the transfer function of the feedback path will be 








3s*+ 21s + 18 
K (s)= S784 Dis + 18° 
. s? + 15s + 18 
\ Wi (8) = Ga pds + 18 











The transfer function of the feedback link can be realized 


"7 by means of a passive RC-circuit [7]. 
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For the parameter values chosen, the mean value of the 


- , tibess squared noise at the output is 
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0146 /t-1) z= 1,02. 
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The impulsive response corresponding to transfer function 


Fig. 3. K(s) is written in the form 





A—b—A AM*—2AbA+2Ab* _,,  B(A—2A) — 
K (t) = ae + a ent +e ot. 18te8 4. 156-8 120-2, 





The graph of K(t) is shown on Fig. 3 (curve I). For purposes of comparison, Fig. 3 also shows the graph of the 
optimal impulsive response (II) as determined by the method of Zadeh-Ragazzini [2], giving, for the same condi- 
tions, a value a which is approximately equal to x. On the basis of [5], for T = 1, xXpnz © 0.984 and k,(t) = 
= 0,492-0.317t + 0.8 5(t)—0.14 5(t—1). It is clear from Fig. 3 that the durations of the transient responses for k(t) 
and k,(t) are approximately identical. 


It should be mentioned that the attempt to obtain a realizable transfer function of the correcting link 
approximately on the basis of an optimal transfer function which does not possess property (5) leads to significant 
difficulties. For example, for the circuit of Fig. 2, let the optimal transfer function have the form: 


As* +- Bs +- 1 


K(s)= Cs? + Ds+1 





Then, the transfer function of the correcting circuit is written as follows: 


— As® + (KC — B)s* + (KD—1)s4+K 
(As* + Bs +1) K ; 





W,(s) = 


This transfer function is unrealizable. If it is approximated to by a realizable function in the natural way, 
i.e., if we introduce into its denominator the term xs* with a small coefficient x, then the new transfer function, 
(8) becomes equal to the improper fraction 


Ww’ (s) -As§ 4+- (KC — B) s* + (KD—1)s+K 
ov om zs* + KAs* + KBs +K f 





and the transfer function of the entire system becomes equal to the proper fraction 


7 P+ AP + Bs +h 





K’ (s) = 


ar st + Cs* 4 Ds +1 


It is clear from this last expression that the system with transfer function K(s) is unstable. 


Generally, with the introduction of small parameters with increase the degree of the polynomial in the 
denominator of the transfer function of the correcting circuit, a stable system can become unstable [8], which 
makes the method of natural approximation inapplicable, independently of the order of smallness of the para- 
meters introduced. 


In conclusion, I wish to thank V. S. Pugachev for his interest in this work. 
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THE EFFECT OF NOISE ON PHASE-LOCKED OSCILLATOR OPERATION 


V. I. Tikhonov 
(Moscow) 


The action of external noise and of inherent fluctuations on 
an inertialess phase-locked oscillator are considered by means of 
the Fokker-Planck equation. Formulae for the mean and dispersion 
of the synchronized generatc:'s frequency are found. The analogy 
is established between the problem considered here and the problem 
of synchronizing a self-excited oscillator in the presence of noise. 


Today, in various devices of radio engineering and automation, widespread usage is made of phase-locked 
oscillators for automatic frequency control (A.F.C.), the theory of operation of which can be elucidated using the 
example of the simplest functional schematic (Fig. 1). 


The harmonic oscillations, u(t) of the fundamental generator 












































ree? (FG) and u,(t) of the synchronized generator (SG), act on phase detector 
FG alt PD Unyy LFF PD at whose output there appears a voltage which depends on the dif- 
uftt) ference of phase of the sinusoids ug(t) and u,(t). This voltage, by means 
Jusftl Ju(t) of low-frequency filter LFF and reactance modulator RM, changes the 
SG RM frequency of the synchronized generator, causing it to coincide with 
the frequency of the fundamental generator. With this, there is estab- 

















lished a constant difference of phase between the sinusoids of the 
Fig. 1. A.F.C. functional schematic. fundamental and the synchronized generators, which provides synchro- 
nized operation of both generators [1]. 


The operation of an A.F.C. circuit is considerably complicated if the presence of noise is taken into account. 
Under practical conditions of A.F.C. operation there occur both external noise and fluctuations which are organically 
inherent in the circuit elements themselves. As is well-known [2-4], anode current noise in the generator tubes, 
random variations in the voltage supply and other causes result in fluctuations of the generators’ frequency (phase). 
Moreover, together with the sinusoid from the fundamental generator, additional noise, n(t), may act on the phase 
detector, particularly if the fundamental generator is spatially remote from the rest of the A.F.C. elements. 


The problem of the effect of noise on A.F.C. operation is, in general, very complicated mathematically 
because of the specific nonlinearity of the A.F.C. differential equation. Although several works [5-8] have been 
devoted to analyses of particular cases of the action of noise on A.F.C. systems, it is impossible to consider, 
even from a qualitative point of view, that the problem has been solved with sufficient completeness. 


In the sequel, we shall limit ourselves to the consideration of a simplified, inertialess A.F.C. system for 
which the block schematic of Fig. 1 applies with the differences that the PD is replaced by a multiplying device 
and the LFF is considered to be an ideal low-frequency filter (its transfer constant equalling unity for low 
frequencies and zero for high frequencies). With this, we shall take into account the fluctuations of frequency 
of the fundamental and synchronized generators and the presence of the noise, n(t). This noise is assumed to be 
norma! with zero mean; the spectral density of this noise is symmetric with respect to the mean frequency of the 
fundamental generator. 
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Posing of the Problem 








If we ignore amplitude fluctuations, we can write the oscillations of the fundamental and synchronized 
generators in the following forms, respectively: 


ug (t) = A; cos D, (t) = A, cos [wet + 9, (t)}, (1) 

Ug (t) coe Ag sin ®, (t) = A, sin [wt a 6, (t)], (2) 
where A, and Ag are constant amplitudes, w,) and w are the mean frequencies and 6 ;(t) and 09(t) are the random 
phases of the generators’ oscillations. 


The normal fluctuating noise n(t), with spectral density symmetric about frequency w), can be presented in 
the form 


n(t) = E(t) cos @(t) = E(t) cos [wft + 4(t)], (3) 


where E(t) is the envelope of the fluctuations and 9@(t) is the random phase. 


The total voltage of the fundamental generator and the noise which acts on the multiplier block can be 
written as 


us (t) + n(t) = A, cos D, (t) + E(t) cos ®, (t) — Z, (t) sin ®, (t), 


where 
E,(t) = E(t) cos (@—9,), £y(t) = E(t) sin (0 —9,). (4) 


By using the same methodology as in [9], we may easily show that the sinusoid E(t) and the cosinusoid E,(t) 
composing the envelope E(t) have normal probability densities with zero means. If we denote the correlation 
function of the noise n(t) by 

ky (t) = 9° (t) cos wt, 


where 0 is the dispersion of the noise n(t), then the auto-correlation and cross-correlation functions of E, and E, 
are defined by the relationships 





h(t) = Ey (t) By (t + *) = By (t) Bg (t + t) = o°r (2), 


5 ae A carr (5) 
h (=) = B (t) & (t + *) = — Bs (t) & (t + *). 





We denote the multiplier's transformation coefficient by ». Then, the voltage at the multiplier's output 
will be 


um(t) = u [u g(t) + m(t)] us (t) = +A, {(A; + £,) [sin (O, — ®,) + 


+ sin (, + ®,)] — Ey [eos (, — ®,) — cos (®, +). m 


Since the LFF is ideal, i.e., its transfer constant is unity for low frequencies (w —w)) and zero for high 
frequencies (w+ Wp»), the voltage at the input to the RM will be 


u(t)= sind, [(A, + £,) sing —E,cosq], p= O,— D,. (7) 


Let the average frequency w of the synchronized generator depend linearly on the voltage at the reactance 
modulator's input, i.e., 


® = w,.— su (t), (8) 


ao 


1161 








where Wap is the average frequency of the synchronized generator for u(t) = 0 and s is the slope of the linear 
portion of the reactance mcdulator's modulation characteristic. 


From expressions (1), (2) and (7) we obtain 


d . . 4 
Fe = G =D, — Dy = (Wyo — «) — su (t) + (6, —4,) 
or 
o= A,— Asing —7-A (Essing — Etcos¢) + 4, (9) 


where A, is the initial tuning of the generators, 


Ag = @m—®, A= < wsA,Ag, ¢ = 6,—6,. (10) 


To compute the statistical characteristics of g from the nonlinear differential equation (9), in which the 
random functions E,(t), E,(t) and ¢(t) enter is, in the general case, quite difficult. By means of the substitution 
x = tan (¢/2), this equation can be reduced to the Ricatti equation and this can then be transformed into a linear 
homogeneous, second-order equation with variable coefficients. However, these transformations do not give any 
simplifications in principle. 


In the sequel we shall apply to equation (9) the apparatus of Markoff processes, which is valid if the fol- 
lowing inequality holds 


2 
isAa(A FE) > eT Taje (11) 





where T;,,, and Tj are the correlation times on n(t) and 9(t) respectively. This inequality holds in many practical 
cases, for example, for radio communication equipment. 
Solution of the Fokker-Planck Equation 
If we have the first-order differential equation 
9 = F [p, E(¢)], (12) 





where €(t) is a stationary random function then, given a condition analogous to (11), the following Fokker-Planck 
equation (10) corresponds to it, defining the one-dimensional probability density for the random function ¢(t): 


10 (9) = — 3 1Ks (9) w(9)] + 5 gee [Ka (9) w(@)1- (13) 


The structural numbers, K; and K,, are computed from the formulae 





K, (9) = Fig, &(t)], Ka(y) = \ (Fie, (1 F Ig, E(t + t)] — Ki @)} de. (14) 


j—s 


For the statistical averaging, denoted by the superscript bar, y is considered as a fixed parameter, not as 
a random variable. 


We now introduce a physically graphical interpretation of the Fokker-Planck equation (13) which we shall 
use in the sequel. 
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In each concrete realization of the process described by equation (12), the trajectory g(t) can have a very 
complicated form. The representative point moves along the t axis like a Brownian particle (i.e., like a particle 
engaged in Brownian motion). If we take a large number of realizations of the random process, we shall have a 
large number of representative points carrying out random wanderings. The points form, as it were, a “gas” under- 
going diffusion, where the density of this "gas" at any point is proportional to the probability density, Each indivi- 


dual representative point constitutes a “molecule” of this “gas.” It moves, but cannot be created or annihilated. 


In correspondence with equation (13), the average velocity of a “molecule” at a chosen point of the ¢ axis 
equals K;(g). By multiplying this by the relative density of the "molecules," w(y),, we obtain the convection 
flux of the "molecules," equal to Ky(y)w(y). The presence of the concentration gradient, 8 w/8 g, also leads to 
the appearance of a diffusion flux, which is proportional to the concentration gradient. The total flux of “mole- 
cules,” or the probability across section ¢ in unit time, equals 


G(p) = Ky (9) w(9) — 5 Ge [Ke (9) w()- (15) 


Thus, Equation (13) can be written in the form 


(9) + F> G(g) =0 (16) 


and may be considered as the equation of probability conservation. 


In the sequel we shall only be interested in stationary distributions, characterized by the fact that the 
probability distribution, w(¢), does not depend on time, i.e., W(g) = 0. Therefore, we obtain for it, from (13), 
the following equation 


|% 


5 Fal Ka (9) w (9) — J (Ki (9) w(@)] = 0. a7) 


Qa 


9? 


By using formulae (5) and (14) for computing K,, one can easily show that, for equation (9), the structural 
numbers equal 


Ki (9) =4,—Asing, K,(9) = \ (4) ko + ky (t) | de. (18) 


—oo 


Thus, as applied to the problem under consideration, equation (17) can be written in the form 


0*w d 7 
oe —_— Fry [(D, — Dsin 9) w) = 0, (19) 
where 
Dy = 240 = 2a) 2A wera (20) 


~ in ° 


It should be mentioned that, in the given case, Ky does not depend on ¢, the fluctuations of the generator 
frequencies and the noise entering additively into the fundamental signal channel, the sole difference being that 
the noise is multiplied by a factor which takes into account the translation of the noise from the multiplier input 
to the input of the synchronized generator. In view of this, we obtain a complete analogy of the problem under 
consideration with the problem of synchronizing a self-excited oscillator in the presence of noise [11]. 


The general solution of equation (19) can be written in the form 


4 
w (p) = CyeDietD 008 » \ e—D—D 008 + diy, 
Cc 
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It contains two constants of integration, Cy and C,, which are defined by two conditions: 1) the condition 


of periodicity, w(¢ + 2%) = w(y); 2) the normalization condition for the probability density in each period, 
Qn 
\ w (9) dp = 1. 
0 
By carrying out the necessary computations we get [11] 


eter 
w (9) = y eDve+D cos @ \ e—Did—D 008 4 di, (21) 


? 
N = 4n%e—"»| Ip, (D) |’. (22) 


Here I ;,,(z) is a Bessel function of imaginary order and imaginary argument. 


With zero initial tuning, D, = 0, the probability density for phase is determined by the expression 


i 
w (¢) _ jale(D° 


and is symmetric with respect to the values gy = 22k (k= 0, 1, 2, .. .), corresponding to stable operation of the 
A.F.C. and the absence of noise. For D, + 0, the form of the curve for w(¢) is asymmetric. It can be found by 
using tables of Bessel functions of imaginary order and imaginary argument [12] and carrying out the numerical 
integration of (21). 


Statistical Characteristics 








For the qualitative behavior of the phase difference, y, in the presence of noise, one can give the following 
graphic explanation, suggested by R. L. Stratonovich. We first assume that the noise, n(t), is absent, i.e., E, = E, = 0. 
Then, equation (9) takes the form: 


. , 23 

og = A,— Asing+ 9. - 
t 

We shall formally consider the quantity s = \o(2)dz as the path followed by some point during time t 
0 


and correspondingly, @ will be considered the point's acceleration. Then, equation (23) can be interpreted as the 
equation of motion of a point in a field of regular external forces which depend nonlinearly on position, with 

the presence of the random disturbance y. If A)- Asin ¢ is considered as the strength of the external force field, 
then the potential of these forces equals 


? 


U =—\(a0.— Asin 9) dp = — Ayp— Acosg. 
0 


This potential has the same form as the potential of a mass point moving on an undulating sloping surface 
(Fig. 2). This surface, for A, > 0 and A > A), has “pockets” in the neighborhoods of the points gy = 42k (k = 
= 0,1, 2,...), these being stable positions of the point in the absence of random disturbances. At some initial 
moment of time, let the point be at "pocket" My. A random disturbance, first of all, moves the point in the 
neighborhood of the "pocket" in a random fashion, causing small fluctuations of frequency and phase, and, secondly, 
more or less frequently takes it from the “pocket,” after which it gradually rolls to the following next lower “pocket,” 
Mg. Less frequently, a random shock will throw the point into "pocket" Mg which is placed at a higher level. The 
. first transition corresponds to an increase of the phase by 2m, the second to a decrease by the same amount. In 
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general, there is a sliding off of points from the sloping surface which occurs 
the more rapidly, the larger is the slope (the initial tuning), the shallower are 
U the “pockets” (the lower is the amplitude of the synchronizing voltage) and the 
| larger are the random disturbances. This means that there occurs a regular 
L change of the synchronized generator's average frequency, causing it to approx- 
7 ~ ? 
, imate to the frequency of the fundamental generator. 





If the random disturbance, %, is very small, then the transitions from Mj, to 
M, and Ms will have low probability and, in practice, need not be taken into 
account. With this, the magnitude of ¢ will experience small fluctuations in 
the neighborhood of stable state My. For small disturbances, independently of 
Fig. 2, Potential contour. their structure, and if condition (11) holds, one can limit oneself to considering, 
instead of equation (9), a linearized equation with respect to the fluctuating 
quantities [13-15]. 





For 4 < A), the “pockets” vanish, due to the overly large slope and, with this, the point uninteruptedly 
slides down, and a synchronized mode of operation is impossible. 


If one takes into account the presence of noise in the channel of the fundamental generator (E, # 0, E, ¥ 0), 
then the slope of the undulating surface will experience fluctuations of the standing-wave type, with random 
amplitudes distributed normally with zero mean. If condition (11) holds, these random pulsations occur very 
rapidly. They also make possible a systematic translation of the point to the average slope of the surface. 


We note that the idea of jumps of the representative points in systems similar to the ones we are consider- 
ing, is contained in [16]. 


We now compute the mean frequency of the synchronized generator. From equation (9) we have 


. : A ; 
o = (A4,— Asin 9) — 7, (E,sine — 2, cos¢) +4. 








According to condition (11), the random functions E(t) and }(t) vary rapidly in time in comparison with 
a(t). Therefore, despite the correlation between E(t), #(t) and g(t) which, in theory, can be taken into account, 
we shall, with an admissible approximation, first average the rapidly varying quantities separately, and then the 
slowly varying ones. We note, by the way, that the same approach can also be used in computing the dispersion 
of the frequency. 


By taking into acc ount that E; = E, = 0 and = 0, we can write 





ss 2n 
9 = (0) — Ssing) =| (A, — Asin) w (¢) de. (24) 
0 


It follows from formulae (15) and (18) that 


: oe d 
(4) —A sing) w() = Ki (@)w(@) =G + 7K. (25) 


If we substitute (25) in (24) and bear in mind that, in the stationary state, the flux G is constant but the 
probability density is a periodic function with period 27, we get that 


¢ = 2nG. (26) 


If we use expression (21) for computing the flux with (15), we obtain the final formula 


- nK, 


on one ee sh tDoy 
7 ae 


(1 —e™) = Ay— D, | Tio, (P) | (27) 
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The results of computations carried out by means of 
this formula by I, G. Akopyan are shown on Fig. 3. 
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We now consider two limiting particular cases. 
1. As is obvious from formulae (20), to the case when 
D-> 0 and D, ~ 0 there corresponds very considerabie noise. 
By using the well-known equality, 1(0) = 1, we get from 
ff 8 formula (27) that 
i+ DxJ 
| ) Thus, in the case of a very high noise level, the synch- 
De§ V/ Dae ronizing signal has no effect, and the average frequency of the 
4 1_| synchronized generator remains equal to its initial value wa. 
ae D=” “3 
0 re Zapf 2. It can be shown that in the opposite limiting case, 
when the noise level is very low, the following dependency, 
Fig. 3. Difference of the average frequen- obtained by Yu. B. Kobzarev and R. V. Khokhlov [17], is valid: 


cies as a function of the circuit parameters 
and the characteristics of the noise. 


0 for Ag <A, (29) 
(A3— A*)” for A> A. 


It follows from this that, for a very low noise level, a normal synchronized mode of A.F.C. operation is 
established if the initial tuning, A,, does not exceed A, i.e., the hold-in range in the given case equals A. If the 
noise level is not very low then, if even the inequality A, = A holds, a synchronous mode of operation is 
impossible in the inertialess A.F.C. system under consideration, and the mean frequency of the synchronized 
generator will differ from the frequency of the fundamental generator. 


© — of 


We now compute the dispersion of the frequency. If we square the right member of equation (9) and average 
the rapidly varying quantities, we find that 


(do— Asin 9)* + (A) (Bt sin*e + Et costg +# = 
= (A, — Asin 9)” + 07 + 203, 
since, as a consequence of (4), 
EF, sin*@ + Ef cos? = +E == o%, v =k; (0) = 205. 
We may therefore write 
v A s e" 
g = (z- 2) +202-+( (do —A sing)" w(¢) de. (30) 


By using formula (25) and computing the integral, we obtain the following result: 


¢° = a 2) + 2o5 + 2rA,G — 
tw aT A | Tip, (P) |"? Yi-av, ((D) Jin, (iD) — Jiziv, (@D) Jin, (iP)1, 


where J; 4 jy (iz) is a Bessel function of a complex order and an imaginary argument. 
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Consequently, the dispersion of the frequency equals 


ot = g*—(g)* =(A =) + 203 + 2nG (dy) — 1) — 


> 
Fr Wap, (D)[? Hha, (D) Lip, (D) —Jiyan, (D) Lan, (DI, (31) 





where 


sinhxD 1 
Bu do -ape | Tip, (D) |? 


SUMMARY 


In the analysis of the action of noise on an A.F.C. system, several assumptions were made which permitted 
the problem to be solved (the phase detector was replaced by a multiplier, the low-frequency filter was assumed 
to be ideal, the noise was assumed to be a stationary Gaussian function with a spectral density symmetrically dis- 
tributed about wy, etc.), The solution of the problem allowed a number of interesting facts to be established, for 
example, the presence of phase variations causing divergences of the average frequencies of the fundamental and 
synchronized generators. 


If these assumptions are not made, the mathematical solution of the problem is significantly complicated. 
Thus, for example, if an actual low-frequency RC-filter is taken into account, the process in the A.F.C. will be 
described by the following equations: 
9 = Ao — su + $, 
RC = +u= = As (A, + Es)sing — Ercos 9]. 


These equations can be reduced to one nonlinear, second-order differential equation. If condition (11) 
holds, its solution can then, in theory, be obtained by a special transition from the second-order equation to two 
first-order equations, with the Fokker-Planck equation then being applied to them [18]. 


If the noise, n(t), and the fluctuations, ¥(t), are slow, so that the inequality inverse to (11) holds, then a 
quasi-static approximation may be used [4]. If inequality (11) holds then, for low noise level, independently of its 
structure, one may use the linearlzation method in the neighborhood of a stable state. 
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THE EFFECT OF BALANCING LOOP COUPLING ON THE DYNAMIC 
PROPERTIES OF AC BRIDGES AND COMPENSATORS® 


V. Yu. Kneller 


(Moscow) 


This paper investigates the dependence of the critical gains 
of the balancing loops of automatic ac bridges and compensators 
on the amount of coupling of the loops. A knowledge of this 
dependency is necessary for choosing the tuning of such instruments. 


In the design of automatic ac bridges and compensators with two balancing organs (Figs. 1 and 2), these 
forming an important class of automatic control devices, it is of great practical importance to make the proper 
choice of the reference voltages’ phases for the phase~sensitive indicators (PSI) by means of which the balancing 
organs form the control signals. As was shown in [1, 2], such devices, in the general case, are two servo systems 
with mutual coupling, and may be presented by means of the equivalent block schematic shown in Fig. 3. As the 
PSI reference voltage phases change, the amount of coupling between the balancing loops also changes and this, 
to a significant degree, changes the system's dynamic properties. In order to provide a tuning of such devices 
which would guarantee the best dynamic properties, one must know how the magnitude of the coupling between 
the loops affects the system's dynamics. The present work is devoted to the investigation of this question. 
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Fig. 1. 1) Phase-sensitive indicator; 2) phase-shifting Fig. 2. 1) Phase~-sensitive indicator; 2) phase- 
network ; 3) servoamplifier; 4) servomotor; 5) preampli- shifting network; 3) servoamplifier; 4) servo~ 
fier. motor. 


* This work was presented at the Sixth Scientific-Engineering Conference of the junior members of the IAT AN 
SSSR on Automatic Control, held January 21, 1959. 
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Fig. 3. 


We run into similar questions when dealing with the design of systems for controlling objects with coupled 
quantities and in designing multi-dimensional servo systems [3, 4] but the existing results do not allow an answer 
to be given for the case we are considering. 


Here, we shall consider systems whose block schematic can be represented as in Fig. 3, where the W's are 
transfer constants whose magnitudes can be chosen at the designer's dption. Among such systems are automatic 
ac bridges with PSI, for any form of bridge scheme, including specially connected balancing organs [5], and also 
automatic ac compensators [2]. Special cases of such systems are systems with antisymmetric couplings [6]. 


In choosing a criterion for judging the dynamic properties of the class of systems under consideration, we 
shall start from the first basic requirements imposed on any automatic measuring device — the requirements of 
providing a given sensitivity to the quantity to be measured and providing a given speed of response.* The giving 
of these quantitiesdetermines the magnitude of the total gains which must be realized in the balancing loops. 
With this, the larger the given sensitivity and speed of response, the larger must these gains be. However, stability 
considerations impose limitations on the size of the balancing loop gains. The critical values of the gains, Kg,, 
are determined by the parameters of the balancing loops and by the magnitudes of the couplings. It is clear that 
the larger are the magnitudes of the critical gains, the greater will be the sensitivity and the speed of response of 
the system which can be provided under essentially the same conditions. It is therefore natural to adopt, as the 
criterion for judging the system's dynamic properties, the magnitudes of the balancing loops’ critical gains. 

Thus, we shall investigate how the magnitudes of the couplings, Wy and Wt, , affect the magnitudes of the 
critical gains, bearing in mind that systems with large K,, are better. A similar problem for one class of bridge 
systems was partially solved in [5]. In this case, it was assumed that the balancing loops were completely identi- 
cal and that there was a 90-degree phase shift between the phases, ¢ and 0, of the PSI reference voltages, thanks 
to which the system could be investigated by the method of complex coordinates. However, with this,one lost 
sight of the interesting, from the practical point of view, modes of operation with different gains in the loops and 
with asymmetric couplings between the loops. The analytical method adopted in the present work not only allows 
these modes of operation to be investigated, but also makes it possible to include a far wider class of systems than 
the class considered in [5]. 


With changes in the reference voltage phases, ¢ and 6, there occur changes both in the coupling coefficients, 
Wye and W‘,, and in the transfer constants, W!, and W) and, consequently, the sensitivity of the instrument to 
changes in He measured quantities also changes. However, a comparison of systems with different g and 9 for the 
purpose of finding the system with maximal critical gains must be carried out under essentially identical condi- 
tions and, in particular, for one and the same value of sensitivity. Consequently, it is necessary to compare systems 


* Here, speed of response is judged, as is customarily done with measuring instruments with automatic balancing, 
as the time required for a full-scale deflection. 
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which are obtained from the original one, not only by changing the coupling coefficients, but also by simultaneously 
changing the gains of parts of the paths from the inputs to the servomotors in such a way that these gains, despite the 
variations in W) and Wy, (with changes in y and @), remain as they were originally. It is possible to take this 
into account by replacing our consideration of the system whose block schematic is given in Fig. 3 by consideration 
of the system whose block schematic is given by Fig. 4. The latter design differs from that of Fig. 3 only in this, 
that before the motors in the balancing loops there are introduced links with transfer constants Why ar, and Why are: 


where kyary and kyarg are respectively, the factors of Wx and Wx which change with changes in g and 0. 


Thus, the problem reduces to the investigation of the changes in the critical gains of the loops of the equi- 
valent scheme of Fig. 4 as y and 9 vary. The characteristic equation of this system can be written in the form 





k e I » ° 
Ww, yo — v2 |_ wW,W,, =0. 1 
[ VU + kK; (s) iif UV: + k Ks (s) | Va" WU ( ) 


In the class of systems under consideration there can be two 
j types of functional dependence of the transfer constants W on ¢ 
’m\—"| Wari | “9 and 6. In the sequel, we shall analyze these two groups separately. 



































by The first group contains ac bridges with balancing active and 
4 passive impedances, placed on one arm, and rectangular-coordinate 
compensators (Fig. 2). 



















































































L The second group contains ac bridges with balancing resistors, 
Wy Wy, distributed on different arms (Fig. 1). 
My t lay, Analysis of Systems in the First Group 
| For systems in the first group, the following relationship is 
“. nod K,(s) valid {1}: 
Fig. 4. W,, = —kcos p, 
Wy, = —ksinp, kya, = cosp. 
Wi, =-+ksine, hAyg, = cose, (2) 
W,, = —kcose, 


where, in the case of rectangular-coordinate compensators, » and € are the angles between the voltages on the 
balancing organs and the corresponding PSI reference voltages (for the circuit of Fig. 2, w = LUM ox, 

e= /UWUoy). In the case of an ac bridge, up = 9 + $,e = 9+ 4—90°, where» is the angle between 
the supply voltage of the bridge circuit and the voltage on the bridge arm adjoining the arc to which the control- 
ling organ is connected. 


By substituting (2) in (1) we obtain 








+ k cos p + ie | [- kcose + saan + sin p sin ek? = 0. 
ri 


By assuming that cosu +0, cose+0, and by using the notation 





wiptlinin’ ane Sea dn — 
k kK, (8) Ki " kk K,(s) ants 
where Ky > 0 and K, > 0, we obtain 
[D, (s) + Ki) [De (s) + Kol + Ki Ka tgpige =0. (4) 
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It is obvious from (4) that a change in the phase of a PSI reference voltage has an effect only on the free 
term of the characteristic equation. Thus, the problem reduces to the investigation of the effect of changes in the 
characteristic equation's free term on the system's stability. For given D,(s), D,(s), Ky and Kg this is easily car- 
ried out by using the well-known methods of automatic control theory. However, our prehlem is to obtain a 
series of general recommendations which will be valid for the class of systems under consideration for different 
values of these parameters. These recommendations must give the possibility of choosing y and @ from the point 
of view of providing a maximal K,, without having recourse each time to concrete numerical calculations. 
For the problem thus posed we impose the following limitations which take into account the properties of the class 
of systems under consideration: 


1) we shall assume that both balancing loops have approximately identical hodographs, i.e., Dy(jw) © Da(jw) = 
= D(jw); 


2) we shall assume that the Mikhailov hodographs of the balancing loops take the form of unwinding spirals 
which do not intersect themselves. 


The first limitation is valid insofar as both balancing loops of automatic bridges and compensators generally 
are made up of mono-typic elements. The second limitation is also not contrary to experience, since balancing 
loops of such automatic measuring devices ordinarily consist of aperiodic, astatic and well-damped oscillatory 
links, and with this the hodograph D(jw) does not intersect itself [7]. 


We now elucidate the qualitative character of the dependence of the balancing loops’ critical gains on pu 
and €, in considering them as Mikhailov hodographs. It is clear from (4) that, to construct the Mikhailov hodo- 
graph of the complete system, it is necessary to multiply the Mikhailov hodographs of the individual loops (curve 
1, Fig. 5,a) and thereafter to translate the resulting hodograph (Fig. 5,b) by the amount KyK, tan p tan € parallel 
to the real axis, to the rightiftan py tane > 0 and to the leftiftan u tane < 0. With this, the system re- 
mains stable if (Fig. 5, a and b),* fortan p tane > 0, 


_ OA - VA 5 
|'gulye| < pay and) tgp tge| < KK,’ @) 


and, for tan p tane = 0, 


; al, 
jlgulge| < rm and| igpige| <1. 
142 (6) 


Here we shall consider only those portions of the Mikhailov hodographs of the balancing loops which include 
the points a;, a and as since, by virtue of the second assumption made, it is precisely these points which form the 
points A;, A, and As which are situated the closest to point 0. As is obvious from (5) and (6), the values of the 
balancing loops’ gains are bounded by the quantity tan yw ande€ . If the lengths of the segments OAj (i = 1, 2, 3) did 
not depend on K;K,, then a single answer would derive from (5) and (6), namely, the larger tan » tan € , the smal- 
ler the critical gains of the balancing loops. But the lengths of the OA; also depend on K,K, and it is therefore 
necessary to explain to which side the points Aj are translated when K,K, increases. 


We assume initially that Kj = K, = K. With this, as K increases, segment OA, always decreases, and, con- 
sequently, OA, decreases, i.e., point A, is translated to the left. 


Thus, the right member of inequality (6) decreases as K increases and we obtain an unequivocal answer, 
namely, that for tan » tan € <0, as yu and € increase the critical gains decrease, and this means that if it is 
desired to increase K it is then necessary to choose tan y tan € as small as possible. The law of change of OA, 
with increasing K depends on the form of the Mikhailov hodograph. Thus, for example, if the hodograph has the 
form of curve 2, Fig. 5,a then, with certain limits of K's values, as K increases the segment OA, can increase. 
However, for further increases of K, beginning with a certain value of it, OA, again decreases with increasing K. 


* On Figs. 5-12, the left side of the figure (labeled "a”) represents the Mikhailov hodographs of the individual 
loops, while the right side (labeled "b") represents the product of the Mikhailov hodographs of the individual 


loops. 


1172 









































Thus, we arrive at the conclusion that also for tan u tan € > 0, if it is desired to increase K then tan y tan € 
‘ should be chosen as small as possible. 
1 the 
r It is not difficult to see that, for K close to 
K,,, when K increases there is a decrease both in 
t OAg and also in OA, and OAs, while for K = Ker, 
oint as OA, reduces to zero so will at least one of the 
; quantities OA, or OAs, i.e., the points A; go to 
class zero from two sides simultaneously. Consequently, 
for K = K,y, the over-all system cannot be made 
stable for any values of y and €. 
De( jw) = 
Thus, for the case considered we obtain the 
: unequivocal answer: as the coupling of the balancing 
irals loops increases, the loops’ critical gains decrease. 
The maximum possible K,, in the loops are realized 
rally Fig. 5. when the system is completely uncoupled. 
ing We now consider the case when Ky # Ky. We 
y show the Mikhailov hodographs of both balancing 
mee loops on one figure (Fig. 6,a). The balancing loops 
u 4 rs oe b j with different gains are represented on Fig. 6,a in 
xdo- such wise that the axes of the hodographs’ ordinates 
urve a do not coincide. As the gain of one of the loops 
rallel a v< increases, the corresponding axis of ordinates on 
re- ii Dal ; fp AX Ol YA Fig. 6,a is translated to the left. Here, OA = Oya - 
. ti" eine 3 , ; A, *Ogag, and with an increase in gain of one of the 
; ‘ loops, for example, Ky, segment OAg is decreased. 
7 : Consequently, for tan » tan € < 0, the larger is 
5) \ty | tan p tan 4 , the smaller is the critical value of 
v Kis. 
NA The frequencies for which the hodograph 
J 
6) M (jw) = [D, (jw) + Ky) [D3 (jo) + Ko] 
, Fig. 6. crosses the real axis at the points Ay and As (Fig. 6,b) 
are found on Fig. 6,a as the frequencies at the points 
nclude of intersection of the individual loop hodographs 
m the ; with the straight line parallel to the axes of ordinates 
he going through the point n which bisects the segment 
2, 3) did 0,0). 
qnety We now consider how OA, and OAs change as 
oe the gains in the balancing loops change. Here, 
several cases are possible. If the gains simultaneously 
son- increase in such fashion that their difference remains 
invariant, then OA, and OAs, as in the previously 
considered case, decrease starting with some values of 
- K; and Kg. Therefore, also in the given case, for 
das increasing K,,, | tan y tan €| should be decreased. 
Ory However, if only one of the gains increases, the law 
na Fig. 7. of change of OA, and OAs is determined by two 
a phenomena: the decreases in the quantities na, and 
.* nag as one of the axes of ordinates is translated to the left, and the change in the magnitude of 0,0). 
al Thus, if the lesser of the loop gains is increased, there occurs a simultaneous decrease both of ayn and agn 
l and of 0,03. Consequently, OAy = Oya; + Ogay and OAs = Oyay + Ogas decrease as rapidly as in the previous case. 
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Therefore it is again desirable to decrease |tan » tan €|. If now the larger of the gains is increased (for example, 
K, on Fig. 6,a) then, simultaneously with a decrease in na, and nag, there occurs an increase of 0,0, and a con- 
sequence of this can be an increase of OA, and OAs. 


In considering the process of constructing the Mikhailov hodograph of the over-all system from the hodo- 
graphs of the individual systems, we may note the following. If the points O, (Fig. 7,a) and ag coincide, i.e., in 


the case when OA; = O (the gain in one of the loops equals a critical value), points A, and Ag will still be found 
to the right of point O, since 


O.a,:a,0, + 0 andQ),a,+O,as; + 0. 


This means that, although the critical value of the gain is realized in one of the loops, the over-all system 
can be made stable if the hodograph Mj(jw)M,(jw) (Fig. 7,b) is shifted to the right parallel to the axis of abscissas 
in such fashion that point O falls inside the loop (inside segment A A ). For this, it is necessary to choose the cor- 
responding magnitude of tan y» tan ¢ > 0. Thus, in the case considered, the system with coupling between its loops 
permits a gain greater than the critical one to be realized in one of its loops. 
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= 4 Ly 
| _ OT 
a, | 0, 

















Fig. 8. 


We now explain to what value it is possible to increase the gain in one of the balancing loops while simul- 
taneously so choosing the value of tan y tan € that the system remains stable. The system can be made stable, 
at the cost of changing tan » tan €, if the hodograph M(jw) = My(jw)M,(jw) has a loop inside of which it is pos- 


sible to place point O at the cost of only one parallel transfer of M(jw). With this, the hodograph will pass through 
the necessary number of quadrants. 


To determine the conditions for which M(jw) will have a small loop, it suffices to consider how, on the plane 


of hodograph M(jw), the points of hodographs M,(jw) and M,(jw) lying between the straight lines Oycy and Oc, are 
imaged. We consider the possible variations. 


1. Point a lies to the left of point n (Fig. 8,a). The point of tangency, k, of the line from 
O, to the hodograph lies to the left of ayn. As the point m moves on the hodograph (Fig. 8,a) 
from point 1 to point a;, point M (Fig. 8,b) will move in the second quadrant and, at point A,, will 
fall on the real axis, since y =~, + 2 = 180° —5 + v2 = 180° (since at Az, 5=y 2). As m moves from 4 tok, 
point M moves into the third quadrant, since y = 180°—5 + gy, but g,>65 and ¢ increases with increasing w. 
After point k , yg begins to decrease and, consequently, point M approaches the real axis, falling on it at point Ag. 
OA, < OAj, since Oya; + Oza, > Oyn-Oyn. As m moves from ag to ag, |v q|>[5] but, since g = 180° + | 5] —| yl . 
M is translated into the second quadrant, then falls upon point A; on the real axis. OAs > OAg, since O,a3 + O24 > 
> Oya + O,a,. With subsequent motion of m from point as, 15] > | gel, and the variation of ¢ is determined by the 
variation of |5|. The latter quantity increases with increasing w. 


Thus, for O,a; > Ogn and O,k, > O2n, there will always be a small loop which intersects the real axis in 
such fashion that, by translation of M(jw), one can always place point O inside the loop. 
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2. Point k lies to the right of line ayn (Fig. 9,a). As m moves from a4 to ag, point M (Fig. 9,b) does not 
fall in the third quadrant since ¢q > 5, but yg decreases as m moves from a; to ag. As w increases, the hodo- 
graph is changed in such manner that after passage through point a, angle o increases in modulus, but retains 
its negative sign. Therefore point M still does not fall in the third quadrant after its passage through point A, (to 
point As). Consequently, there is no small loop. If point n is found to the left of point a, then the motion of M 
with increasing w will always be determined by the character of the variations of 6 and the hodograph M(jw) will 
have a monotonic character but will contain the origin of coordinates, O, only so many times as the hodograph of 
the individual balancing loop contains point O,. It is also impossible to make such a system stable by means of 
a proper choice of tan w tan €. Thus, a loop appears only in the case when point k lies to the left of line ayn. 

















Fig. 9. 


We have thus shown that, in the system considered, it is possible to choose the gain for one of the balancing 
loops so that it will exceed the critical value. 


It was assumed by us that 


Oe ST Ds(s) st 
Ki kkK,(s) ’ Ka kK, (8) 





We now change the polarity in one of the balancing loops, for example, in the first loop. Then, 





D, (s) 1 ; 
=> ( 
Ki ak kK, (8) ies 
and the characteristic equation takes the form: 
[D, (s) — K,][Dz.(s) + K,.] — K,Kztgutge =0. (7) 


The necessary condition for stability is 


tan p taneé < —1. 


As before, we shall assume that D,(jw) = D,(jw) = D(jw), and we shall consider several cases. 


By proceeding as we did above, constructing hodograph M(jw) from the hodographs of the individual systems 
and using the analogous reasoning, we can show that, in this case, there will be a segment of the axis of abscissas 
which will be contained by the hodograph M(jw) the necessary number of times for stability (on Fig. 10,b, this 
segment is AjAy). 
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Fig. 10. 
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Fig. 11. 


In order that the system be stable, it is necessary so to choose the magnitude of tan y tan € that point O 
will lie on the interior of this segment. 


OA, < K,K,|tgutge| < OA, (8) 
4 2 
K,K, < K,K,|tgpige| < h® + aes . (9) 


Thus, in the case under consideration, it is possible, by choosing the proper magnitude of |tane tany| 0, 
to make the system stable for arbitrarily large gains in the individual loops, if these gains satisfy the conditions 





10,> 24" > k,, (10) 


where I is the projection on the axis of abscissas of the hodograph point with the maximum positive ordinate 
(located on the first turn). 
2, Ky < Ker. O,n > Ky. 


As follows from what has been constructed (Fig. 11, a and b),in this case the system can also be made 
stable by choosing the proper value of tan tan € 
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Fig. 13. 


3. Kg < K.,, Oyn< Ky. 


As point m moves from point O (w = 0) (Fig. 12,a) the phase angle immediately begins to increase and point 
M (Fig. 12,b) moves into the third quadrant. With this, a small loop is not formed, and the system cannot be made 
stable by properly choosing tan p tan €. 


The results obtained can be presented graphically on one figure (Fig. 13) if the signs before Ky or Kg in (4) 
or (7) are taken to refer to Ky or Kg themselves and, in correspondence with this, if the values of Ky and K, are 
laid off on the negative or positive semiaxes in the K,K, plane. Thehatched region in Fig. 13, excluding the 
boundaries and the coordinate axes, is the region of values of Ky and K, for which the system can be made stable 
at the cost of properly choosing tan y tan €. Boundary c-a-c of this region, for large K, asymptotically approaches 
the line e-e which is parallel to line d-d which bounds the region on the left. The magnitude of 1m depends on 
the form of D(j w and is proportional to ag (Fig. 10,a). It is clear from Fig. 13 that the system can be made 
stable for arbitrarily large gains in the balancing loops if these gains are related by (10). For that, it is necessary 
to choose the magnitude of tan » tan € which satisfied (9). Thus, the maximum gains are realized in coupled 
syste ins. 
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However, system operation with gains exceeding the critical values is attended by a number of disadvantages 
which are so essential that, in practice, this mode turns out to be undesirable. As is obvious from Fig. 13, with 
operation in this mode (for example, at point A of Fig. 13) the values of the loop gains are bounded both above 
and below. With this, the range of permissible variations of the gains constitutes a small part of their absolute 
values, when they are sufficiently large and, consequently, for relatively small oscillations of the gains the system 
can become unstable. More than that, with operation with very high gains, it is necessary that p and € be strictly 
constant in order that (9) be satisfied. Computations show that when it is desirable to realize loop gains which 
exceed the critical values by factors of from 5 to 7, it is necessary to establish, and to maintain, the values of 
uw and € with an accuracy of several minutes. When we bear in mind that, in actual systems of the class con- 
sidered, there are always oscillations of uy, € , Ky and Ky relative to the established values, we are led to the con- 
clusion that operation with K; and Ky larger than the critical values turns out to be impossible in practice. Thus, 
the drawback is not the presence of coupling between the balancing loops, but the inconstancy of the couplings. 
Operation with gains less than the critical values, in the case when one of the gains is “negative,” has no super- 
iority in comparison with operation with both gains "positive." Besides, in the first case, one of the gains is also 
bounded below. Therefore, it is advantageous to operate in the first quadrant of Fig. 13 with “positive” gains less 
than the critical value. With this, the maximal loop gains (close to Kcr) are realized in coupled systems. 


Thus, by taking into account purely technological facts (the accuracy of establishing » and € and the 
possible oscillations of Ky and Ky), we conclude that the optimal (in our understanding) system turns out to be 
the uncoupled system (tan ytan« = 0). The most advantageous choice is to take yp = 0 and e = 0 since, with 
this, the value of tan y» tan € will be minimal for the possible oscillations in 4 and €. Moreover, with this the 
sensitivity of the measuring portions of the system to variations in the measured parameters will be maximal. 


Analysis of Systems in the Second Group 





For these systems, the following relationships are valid [1]: 


i -—kcosp, 

W,,, == k3 cos(p — &), kyari= cos p, 

W,, = ksine, kvarg= COSE, (11) 
W,, = — kgsin (e — &), 


Ve 


where € is the phase angle of the measured complex impedance. By substituting (11) in (4) we obtain 





[— k cos + st. [- kg sin (e — &) + oo | — kgksin e cus (p — €) = 9. (12) 
k Ky (s) k Ke (s) 


By assuming that cos » = 0, cos € = 0, and using the notation 








4 ae D, (s) 4 a Dz (8) 
kk Ky (s) Ky kegk Ke (8) Ks 
where K; > 0, K, > 0, we obtain 
{Ky 4- Dy (s)) | Ke BRE —*) + yy (8) | + KK, OSU) tee = 0, (13) 
f cos € ’ COS 2 





From a comparison of (4) and (13) it is clear that, if we consider the term A, sin(S—€) as the gain of 
cos € 
the second balancing loop, denoting it by K,', then all the previously made conclusions relative to the magnitudes 
of the gains in the balancing loops which can be realized with the proper choices of coupling coefficients remain 
in force for bridges of the second type as well. All remarks previously made anent the choice of the quantity 
cos (uu —€) 


K,K, tan p tan € now relate to the magnitude of KyK, a ¢ tan €, 
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On the basis of the relationships given before we can draw the conclusion that, in the given case also, it is 
advantageous to operate with “positive” Ky and Kj, not exceeding the critical values, and that to realize, with 
this, the greatest loop gains, it is necessary to try to decrease [cos (u —€ ) tan €]/cos y as far as possible, i.e., to un- 
couple the system. Whence come immediately the recommendations as to the choice of € and y. 


It is interesting to note that, in the case considered, it is possible, with the proper choice of the coupling 
coefficients, to uncouple the direct feedback loop in one balancing loop, choosing € = € [cf, (13) and Fig. 3] 
and, with this, the bridge will continue to operate normally, since there exists a second coupling via the first loop 
and the inter-loop coupling. 


Thus, the analysis given allows us to draw the conclusion that, in the design of automatic ac bridges and 
compensators, it is necessary to try to decrease as far as possible the couplings between the balancing loops. 


Received March 30, 1959 
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EXPERIMENTAL DETERMINATION OF AUTOMATIC CONTROL 
SYSTEM LINKS’ TRANSFER FUNCTIONS BY MEANS OF STANDARD 
ELECTRONIC MODELS 


L. N. Darovskikh 


A method is suggested for determining the coefficients of transfer 
functions by means of standard electronic models (analogs). A table is 
provided which allows one to set up a system of linear algebraic equations for 
computing the coefficients of the transfer function with, or without, a poly- 
nomial in the numerator, and also when there is only a polynomial in the 
numerator. Examples of the application of this method are given. 


The existing methods for determining the transfer function of a link, or of a set of links, of an automatic 
control system (ACS) [1-7] require a great amount of processing of experimental results, or are applicable only in 
special cases [1-3]. The method presented here is characterized by a minimum number of computations subse- 
quent to the experimentation, reducing to the solution of a system of linear algebraic equations. 


Posing of the Question 





The idea of our method consists in the transformation of the transfer function 


B 
W (pr) = 5, (1) 


where 
A(p)=14+4-a,p-+-agp?+...-+anp", B(p)=1+ byp+ bap? +... + Omp™, 


by means of a constructed scheme which includes a solving block of standard electronic analogs. This trans- 
formation renders it possible to obtain a system of linear’algebraic equations for the computation of the unknown 
coefficients of transfer function (1). 


The transformation has the goal of obtaining a new transfer function with constant terms in the numerator, 
which will be a linear function of the coefficients, aj and bj, of transfer function (1). The size of a constant term 
is determined by the indications of a voltmeter connected to the output as the gain of the circuit obtaining the 
investigated portion of the ACS. 


Setting m < n, we consider the expression 


W, (p) = — <1 (p) — 11. (2) 
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We now carry out the transformation: 


(by — a4) + (bz — a2) p + (bs — a) p? 4-. .. + (—,) p™™ 
A(p) é 





Wp) = - 
A circuit with such a transfer function has the gain 


ky = — (dy) — a). 


A circuit with the transfer function 


Vs (p) = — <1, (p) — I 


has the gain 
Ky = kya; —a@ + by, etc. 


The construction of a circuit for carrying out the transformation cited presents no difficulties. 


If the transfer function contains a gain k which differs from unity then, with the introduction of a nomial 
unit input signal, 1° = 1/k, the link’s transfer function can be considered as a transfer function with unit gain. 
Therefore, the gain k does not enter into consideration; the nominal unit, 1°, will be called simply the unit in- 
put signal. 


Determination of the Time Constant of the Simplest Transfer Function 





Let a link of an ACS be characterized by the simplest transfer function 


: l 
W (P) = 77>" (3) 
In this case, expression (2) gives 
T 
Wi (P) = 79° 


When a unit step input voltage is applied to a circuit with twansfer function Wp), a voltage equal to T is 
established at the output. The operations of algebraic addition of W(p)—1 and multiplication by —1/p can be 
implemented by means of operational amplifiers in the set of standard analogs. 


Figure 1 shows the block schematic for determining the time 
constant, T, of a link with transfer function (3). To the circuit's 
2 input there is applied a voltage of 1 volt (in the general case, m volts, 
“I where m is the scale factor). 


The output of the link under investigation is connected to the 
input of operation amplifier 2 which operates in the summing mode. 
By means of compensator 3, a null voltage is established at the out- 
put of amplifier 2. Verification of the compensation of the signal at 
the output of the investigated link is implemented by switching block 
2 to integration. 





























Input 





Fig. 1. 


When zero is attained at the output of block 2, the input signal is taken off. Block 2 is switched to inte- 
grating. A unit step voltage is applied to the circuit's input. Obviously, a voltage is established at the integrating 
block's output equal, in volts, to 


ky =T. (4) 
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Determination of the Coefficients of Transfer Functions of the Form W(p) = 1/(1 + 
Tt 
+ Q)p + agp ) 


Figure 2 shows the block schematic for finding the coefficients of the transfer function 








i 
WIR Precres (5) 





For this there are required two operational amplifier blocks and two compensators. 


A one-volt signal is applied to the circuit's input. 
A voltage is established at its output, also equal to 1 volt. 
i> Block 2 is in the summing mode. At its output, by means 
of compensator 4, a null voltage is established, which cor- 
K, responds to the passage of the input signal through a link 
with transfer function 





















































Input wip) J re 





W, (p) ois i | Pp (a; + asp) 


i+ap+ap = | 1+ a:p+asp*° 








4, % 
Fig. 2. 
The input signal is then taken off and block 2 is switched to integrating. When a unit step voltage is applied 


to the circuit's input, a voltage is established at the output of integrating block 2 which corresponds to the passage 
of the input signal through a link with transfer function 








— _ _Plaitasp) (_ 1\_ 41+ %p 
Ws (P) = Cyr al >)= 1 + aip + agp* 
and which equals 
ky = a. (6) 


This voltage is applied to the input of block 3, which is in the summing mode. By means of compensator 5, a 
null voltage is established at the output. This corresponds to the passage of the input signal through a link with 
transfer function 


_ i + Gap __  (@p — a, ky) p — kyagp* 
We) = taptar = “Thapar 








Then the voltage is at the circuit's input is taken off and block 3 is switched to integrating. At its output, 
when a unit step voltage is applied to the circuit's input, there is established the voltage 


ke = kya, — a2, (7) 


corresponding to the passage of the input signal through a link with transfer function 


_ ___(@g — a)ky) — ky agp 
Wa(P) = 1+ aip + agp" 





Equations (6) and (7) form the system 


k, = 4, k. = kya, — ap, 
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whose solution gives the values of the unknown coefficients of the transfer function: 
ay = k,, Us = ke — kp. 


Consequently, for determining the coefficients a; and a, of a given transfer function (5), only two voltage 
measurements need be made. The subsequent work reduces to the solution of linear algebraic equations. 


Obtaining the System of Algebraic Equations for Determining the Transfer Function 
Coefficients in the General Case 











Let 


1+ bsp +bep?+.--+bnp™ Bip) 


Ww : = ‘ 
(P) 1 +-a)p + agp? + ...+4,p" A(p) 


(8) 





In solving the problem of obtaining a system of algebraic equations for determining the coefficients of a 
transfer function of the type cited, we shall understand by addition, and by multiplication by —1/p, the operations 
which were cited out in the foregoing samples. 


Letting m < n, we consider the expression 








1 
W, (p) = — > |W (p) —1]. 
We carry out the transformation: 
W, (p) = — tO et (bs — as) P+... 
. + (b,, — a,,) pt 4. Aa 2 + (— a.) p> 
A (p) : 


The constant term, (b, —a,), in the numerator is determined as the circuit's gain, ky, when the transfer 
function is w,(p): 


— ky = b; —a,. (9) 


We now consider the expression 


1 r 
W(p) = — P [Ws (p) — Ay). 


After transformation, we have 
R (by — ag + ayky) + (bs — ag -+ aek;) p+... + (Oy, —4y_ + 2m—, a) = 
W,(p) = Ap) es. 


(— a, + an) ky) y= + a,,ky a 
A(p) , 





oot 





The constant term in the numerator is determined as the gain, kp, of the circuit with transfer function 
W2(P): 


k, — be > ag -{- a, k,. (10) 


We next consider the expression 


W3(p) = — ; [W. (p) — fy). 
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We carry out the transformation: 

















ey - » ~-} _ & Pe die oielies " k of ee 
Ww, (p) a ol (by - @y -+- Agky thy) 4 we ag +- agk, -— agke) p + - 
tba, ee ky 1 An—g ke) sg 
a A (p) ; 3 ieee 


(—@,+¢,_, hy an» ke)p” . 4- (a4, hy as 


™ + (a ~_1 ke) iin a, kop 
A(p) 














The constant term in the numerator is determined as the gain, ks, of the circuit with transfer function W,(p): 
ties ks = hy — a3 -}- A, ky — aykg. (11) 


By noting the regularity of the formation of the coefficients kj, and by applying mathematical induction to 
equations (9-11), we can set up the following system of algebraic equations: 


for m <n 


ky + by = ay, 
k, — be — kya, — ag, 
ky -|- by = kya, — kya, + as, 
ke, — by = kegny — heag 4- hyag — a, 
(12) 
Ken + (— 1) lam = hem—1 yy — m—2 eg + mg 3 — . - 


..- +(— 1)"k,am—, +-(— 1)™"—lan, 
hy, = a a, — Repies Ag -+- kn—s ag eee! 4. (— 1)"k,an_; -- (— 1)"-! an; 
Knam = Kn4m—1 _ Kins m—24g 4-.+.+ (— | Baio Kimi @n + (—1)"+"—'kpay 


form> n 


ky +b, =a, 
ky — by = kya, — a, 
ks + bg = kyay — kay + a3, 
ky — by = kya, — kya, + kyas — ay, 
ky + bg == hyd, — egy -| gg — Kyay + 45. 


Ken 4 (— 198! by = Ry, @y — hn 2 Ae -| kn-nQg—... 
+ (—1)"k,a,_, -|- (—1)"-"' a, 
ky, +(— 1)" bin , Heyyy — km. a, +- km—34s ees 
+ (— 1)" Km—ng-1On—1 4° (— 1)"—'km—nn, 
Keng = Kin, — htm, dg-| kKm—2dg—.-. 
-}- (— 1 a Key, n424yn— -|° (— he fon. 1-1 An, 


(13) 


km j-n = km -n—1! ay — kn i-n--2 Gg ees 
|- (--- 1)" j-n km $1 Mn. : | (-—— 1)" j-u—1 k.., Ay. 














By taking as shown the coefficients ky, ky, ks, . . . and by solving system (12), for m < n, or system (13), for 
m> n, each of these systems consisting of n+ m algebraic equations, we can find the unknown coefficients of the 
transfer function of the ACS link under investigation. The coefficients ky, ke, ks, . ... are determined experimental- 
ly by the methodology explained above. 






If it was necessary to have two compensators and two operational amplifiers for the determination of the 
transfer function coefficients when there was only a second-degree polynomial in the denominator, then in the 
case of m + n unknown coefficients of transfer function (8) it is required to have m + n compensators and m + n 
operational amplifiers. The circuit in which this apparatus is connected remains essentially as it was before. 


We now consider two particular cases of formula (8). 





): Case 1. B(p)= 1. The system of algebraic equations has the form: 
. ky =a, 
k, = ka, — a, 
ee . ssid (14) 
Ken = ken, — hyo, 4 .. . -} (— 1) Ayan, 4- (— 1) oan. 


In considering this system of equations, we can draw the conclusion that, generally speaking, in the deter- 
mination of the coefficients of a transfer function with no polynomial in the numerator, a knowledge of the order 
of the polynomial in the denominator is not necessary, since the coefficient a, is obtained by means of the first 
measurement of voltage, the coefficient a, is obtained after the second voltage measurement, ag is obtained after 
the third measurement, etc. Thus, the process of experimentation can also elicit the order of the denominators 
polynomial (of course, within the limits of experimental accuracy). The method for determining the coefficients 
when there is a polynomial in the transfer function's numerator necessitates either previous knowledge or eliciting 
of the degrees of the polynomials in the numerator and denominator, which must be considered a disadvantage of 
the method presented here. 


Case 2. A(p)=1. The system of algebraic equations has the form: 


k, -4 b, = 0, 
k, —b, =0, 
ks + bs = 0, 


<@: ee «6 6.9 


km +- (— 1)" "bm, = 0. 





The case when the transfer function has no polynomial in the denominator is of great practical interest, 
since by measuring the voltages ky, Kg, ks, . . . it is possible to determine in the most direct way all the unknown 
coefficients of the transfer function. From this the result derives that it is advantageous in experimenting to 
replace a link with transfer function W,(p) = 1/A(p) by a link with transfer function W,(p) = A(p) which can be 
implemented by fully automatic means.° 


For convenience in using the method presented here for determining the transfer function, we provide a table 
below which allows one to set up the system of algebraic equations for the determination of the coefficients of 
the transfer function of the link under consideration, independently of whether or not the operator p is contained 
in the numerator or denominator of the transfer function. The table is made up for transfer functions containing 
m+n = 10 unknown coefficients. The basis of this was the determinant set up from the system of algebraic 
equations in (13). We now consider two examples of the use of the table. 


Example 1. It is required to set up the system of algebraic equations for the transfer function of the form 


1 
1 + ayp + agp* + aap - 





W (p) = 





*It is well known that an operational amplifier of an electronic analog which is shunted by a negative feedback 
path which is a link with transfer function K(p) has the transfer function W(p) * 1/K(p). 








The transfer function contains three unknowns; consequently, it is necessary to set up a system of three 
equations. For setting up this system, only three rows of the table are required. We obtain one equation from 
each row, resulting in the system 


a, == ky, 
hyay — Gg == keg, 
haa, — kyaq + a3 — key, 


Example 2. The transfer function has the form: 





1 + ap + agp* + asp" 


The number of unknowns is five. Consequently, for setting up the system it is necessary to use five rows of 
the table: 


from the first row 


hi 
SS 
5 
es 
we 
‘ 
i 
, 
is 


ay = ky + by, (16) 
from the second row 
ky aq —@2 = ky —by, (17) 
from the third row 
Kpay —kyag + as = ks, (18) 
from the fourth row 
Kg —Kgag + kyas = ky, (19) 
from the fifth row 
Kgts —Kgag + leas = Ke. (20) 


Equations (16-20) form a system of algebraic equations whose solution gives the unknown coefficients of 
given transfer function (15). 


For setting up the systems of algebraic equations, m + n rows of the table were used. The equations set up 
from the following rows reduce mathematically to identities. For the purpose of checking the experimental results 
it is advantageous to introduce still another equation, the (m+ n+ 1)'st. For transfer function (15), the check 
equation is 


kay gag + Kgag = Ke. 


To obtain coefficient kg it is necessary to connect one additional operational amplifier in the experimental 
setup. In practice, the control equation must hold within the limits of accuracy of the given experiment. If the 
left member of the equation differs too sharply from the right member, it is necessary to repeat the experiment 
more carefully, using other input voltages for this. 


Errors in Practical Application of the Method 





The accuracy of the determination of the unknown coefficients of transfer functions depends on the magni- 
tude of the zero drift of the operational amplifiers, on the error of the integration process and also on the error of 
the measuring instruments (the voltmeters). 


The error due to operational amplifier zero drift depends on the type of dc integrator and, in many types of 
electronic analogs, is practically zero. The error is small and can easily be taken into account. The error in the 
parameters found by the method presented will depend, to the greatest degree, on the accuracy of the voltage 
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measurements by the dc voltmeters, This accuracy will, in practice, not be greater than 0.5%. With an increase 
in the number of unknown coefficients, the number of working units also increases, and some error is unavoidable. 


Since the errors will have a fundamentally random character, repetition of the experiment several times and 
averaging the obtained values of the coefficients can increase the accuracy of the results, and their reliability can 
be estimated by statistical processing. Verification of the degree of identity of the check equation might serve as 
a sufficient criterion for judging the accuracy of the results. 


Received June 30, 1958 
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MINIMIZATION OF THE BOOLEAN FUNCTIONS CHARACTERIZING 
SWITCHING CIRCUITS 


M. A. Gavrilov 
(Moscow) 


The paper considers a method for minimizing the Boolean 
functions which characterize switching circuits, the method being 
based on an analysis of the conditions to be realized by the oper- 
ation of the switching circuits. 


One of the essential problems of the modern structural theory of relay (switching) devices is the problem of 
relay circuits minimization, i.e., decreasing the number of elements in them. Since switching circuits can be 
characterized by Boolean functions, the aforementioned problem coincides with the problem of minimizing the 
latter but, however, has its own specific peculiarities. 


A large number of works have been devoted to questions of minimization. These questions were already 
touched upon to a significant extent by Claude Shannon in 1938 [1] with the publication of one of the first works 
in the domain of switching theory. The first regular method was described in [2] in the form of the so-called 
“minimizing charts." Work in the minimization domain progressed considerably with the publication of the papers 
of W. Quine [3, 4] who first formulated the minimization problem in its present form. A number of subsequently 
published works [5-9] were based on Quine's ideas. The majority of them consider the problem of minimizing 
Boolean functions as the problem of finding those of their terms which are “neighbors” [10], i.e., which contain some 
letter both nonnegated and negated. By using the well-known relationship 


ju + fu=f(u+2) =}, (1) 
and also by successively applying the well-known operations 


r+4+-o(z, y,...,w) = 27+ /(0,9,... w) (2) 


11a + fot + fife = fit + fox, (3) 


to the terms of the Boolean functions, one reduces the number of terms in the functions and the number of letters 
in the terms until the point when the operations just cited can no longer be applied. The Boolean functions thus 
obtained are considered minimal. 


In the works of various authors [10-13], various methods have been suggested for placing the terms of Boolean 
functions in such fashion that it becomes possible to determine graphically the terms’ “neighbors.” There have 
been described various types of tables and their transformations, based on relationship (1), subsidiary means in the 
form of minimizing charts or stencils allowing one to find and compute the number of “neighbors,” etc. An 
original method of minimizing was suggested by K. K. Voishvillo [14] which solves the minimization problem by 
finding the so-called “closing sets.” Most of the works cited deal with the problem of minimizing the structure 
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of nonsequential switching devices and do not deal with the minimiza- 
tion of sequential switching devices although, as will be shown in the 
sequel, taking the special features of sequential devices into account 
significantly simplifies the minimization procedure. 


Sena 





In the present paper we consider a somewhat different approach to 
Fig. 1. the solution of the minimization problem, namely, an approach based 
_ on the analysis of the conditions of operation and nonoperation which 
the switching circuit is to realize. The same approach is applied both to sequential and nonsequential circuits, 
and makes it possible to significantly simplify the process of minimization itself.* We understand by the latter, 
as do a number of other authors, the problem of obtaining, based on the conditions of operation of a device as 
described in the form of a switching table or by Boolean functions given in normal form, the same form but now 
containing the least possible number of letters. 


A special feature of Boolean functions which characterize switching circuits is that, in contradistinction to 
ordinary Boolean functions, they may have added to them certain nominal terms, deriving from the operations 
specified for the switching device. To explain this, we considerin more detail the operation of a switching cir- 
cuit. 


In its generalized form the function of a switching device is to transform a combination of discrete stimuli 
applied to its input to a definite combination of these same stimuli at its output (Fig. 1). 


For nonsequential switching circuits, i.e., those in which no facility is provided for any temporal sequence 
of stimuli application at the input or obtaining a stimuli at the output, the circuit's conditions of operation can 
be completely given if each combination of stimuli applied to the input is put into correspondence with a definite 
combination of stimuli at the output. 


We shall denote the presence of a stimulus (for example, the closing of a circuit, or the appearance of a 
high voltage) by unity, and the absence of a stimulus (for example, the opening of a circuit, or the presence of a 
low voltage) by zero. Then, the conditions of operation of a switching circuit can be given in the form of a 
table of binary numbers (the so-called “switching table”) in which, to each binary number which corresponds to 
a combination of stimuli at the inputs or, as we shall say in the sequel, a combination of input “states,” there is 
set up on the same row a binary number which corresponds to the combination of stimuli at the outputs (the com- 
bination of output "states”). 


Such a table of binary numbers is shown in Fig. 2,a for a particular case of the conditions of operation of a 
switching circuit with four inputs and two outputs. For each of the. circuit's outputs all of the combinations of the 
input states may be divided into two groups: one which corresponds to ones at the output and the other correspond- 
ing to zeroes at the output. Figures 2b and 2c give the tables of binary numbers for each of the outputs, Z; and Z,, 
separately. 


We shall call the combination of the input states which correspond to ones at a given output as its “working” 
states, and shall denote the set of them by the symbol Fy. The working states define the conditions of output oper- 
ation (conditions of circuit closure, or the appearance of a high voltage), They can be described by a Boolean 
function whose terms characterize the working states at the inputs. Such a Boolean function is easily obtained if, 
for example, we write each of the combinations of working states (inverted for a zero, not inverted for a one) in 
the form of a product of letters corresponding to the circuit's inputs and then take the sum of such products. Thus, 
for example, for the operating conditions of the table in Fig. 2 we shall have:* * 





* what is presented below is a generalization, development, and the mechanical basis of, work previously published 
by the author on the minimization question. The method described, and its application to sequential circuits, was 
first mentioned in [15]. It was then further developed, as applied to the same devices, in [16] and, finally, as 
applied to nonsequential circuits, was mentioned in the author's paper, "Developmental work on the structural 
theory of switching circuits in the Soviet Union," given at Harvard University in April, 1957 at the International 
Symposium on Switching Theory. In a generalized form, this method was presented by the author at a joint 
seminar on switching theory of the IAT AN SSSR, the Laboratory for scientific problems of line communications, 
and the Mathematical Inst. AN SSSR, in April, 1958. 

** Below, as is done in many works on the structural] theory of switching circuits, the conjunction sign is replaced 
by the sign for algebraic multiplication, and the disjunction sign by the sign for algebraic addition. 
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for element Z; 
Fiz) = abed + abed+abed+abed; | (4) 
for element Z, 
Fyz, = abed + abed + abed. (5) 
It is easily seen that, with such an algorithm for making the transition from the switching table to an analytic 


form of writing the operating conditions of a relay device, one obtains expressions in the form of Boolean functions 
written in the so-called “complete normal disjunctive form.” 
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Fig. 2. 


The second group of conditions, namely, those combinations of states which correspond to values of zero for 
the given output of the switching circuit, will be denoted in their entirety by the symbol Fp, and will be called 
"forbidden." Just as for the working states, the Boolean function characterizing the forbidden states will be obtained 
by writing each of the forbidden states in the form of a product of letters, corresponding to the circuit's inputs, and 
then taking the sum of these products. 


Since, as was shown in [15, 16], the operating conditions of the output can be expressed in terms of the for- 
bidden states ty taking the inverse of the Boolean function characterizing the latter, the Boolean function for the 
combinations of the forbidden states, if used in writing analytically the conditions of operation, will be obtained 
in the so-called “complete normal conjunctive form." 


For a complete grasp of the operating conditions of a switching circuit, it is necessary to know the nature of 
the remaining possible combinations of input states which, when added to F, and Fy, make up the total of 2" combi- 
nations (where n is the total number of inputs). These combinations can be divided into two categories. One of 
these is that of the combinations known to be impossible for them to exist, i.e., for no conditions of operation 
of the objects supplying stimuli to the switching circuit can they occur. Obviously, any term of the Boolean function 
corresponding to one of these states will, when written in the disjunctive form, always have the value zero, and when 
written in the conjunctive form, the value one. We shall call such combinations of states “unused,” and will denote 
the set of them by F,.. 


To the other category of states belong those of which it is known that it is indifferent for the operating con- 
ditions of the switching device whether a one does or does not appear at the output when they appear at the input. 
Such combinations of states will be called “indifferent,” and the set of them will be denoted by Fp: 


It is obvious that terms of a Boolean function which characterize any unused or indifferent combinations of 
states can be added to the Boolean function characterizing the set of working combinations of states without dis- 


turbing the operating conditions of the switching circuit: when the Boolean function is written in disjunctive form, 
these terms will enter in as factors. 
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The addition of terms corresponding to unused combinations of conditions will not change the initial 
switching table, since each such term equals zero, and the addition of terms corresponding to indifferent combi- 
nations of states will give one of the permissible variants of the switching table. The term of the Boolean 
function which correspond to unused or to indifferent states are the nominal terms which were mentioned earlier. 


All further reasoning will apply to Boolean functions written in the normal disjunctive form called, for 
brevity, simply the “normal form.” However, all the minimization methods considered will apply equally when 
the operating conditions are expressed in terms of Boolean functions written in the normal conjunctive form, 
it being necessary only to interchange F; with F, and to interchange zeroes and ones. 


When Boolean functions are written in the normal disjunctive form, all sets of combinations of input states 
corresponding to the operating conditions of some output of the switching device can be written in the form: 


Fiz) = Fi + (Fat Fp), (6) 


where the terms in the brackets are nominal ones, i.e., those which can,at will,be included or not in the Boolean 
function characterizing the operating conditions. 


It is clear that the switching table is realizable if, for each of the outputs, no combination of input states 
is simultaneously in Fy and in Fy.* 


We shall say that a Boolean function, written in complete normal form, “realizes* sme one of the combi- 
nation of input states of a switching circuit, corresponding to the conditions of operation, or nonoperation, of some 
circuit output, if the combination of letters in one of the terms of this Boolean function, when the noninverted 
(i.e., nonnegated) letters are replaced by ones and the inverted letters are replaced by zeroes, is contained in the 
row of the switching table which corresponds to this combination, 


For a Boolean function written in normal form, a combination of input states will be “realized” by one of its 
terms if, when the aforementioned replacement of letters by zeroes and ones is made, the combination is contained 
in part of the row, or the complete row, of the switching table which corresponds to this combination. We note 
that a term of a normal Boolean function can realize simultaneously several combinations of input states. 


By the “length,” 72, of a term in the normal form of a Boolean function we shall mean the number of letters 
entering into it, and by the “length,” L, of the normal form itself we shall mean the total number of letters in it. 
Then, the process of minimizing a Boolean function or, as we shall say in the sequel, of obtaining its “minimal 
form" can be reduced to decreasing the length of its normal form. It is obvious that such a decreasing of length 
can be carried out only up to the point where the Boolean function still corresponds to the operating conditions 
of the switching circuit, i.e., will not realize any of the combinations of states entering into Fy. 


Thus, the minimal form of a Boolean function is that normal form of it which realizes all the combinations 
of states in Fy and does not realize any combination of states in Fy, whereby a decrease in length of any of its 
terms by at least one letter, excluding its complete expunging, would lead to a condition wherein the Boolean 
function would realize at least one combination of states which is simultaneously in F; and Fy. 


It is obvious that if a Boolean function is written in complete normal form its terms f will have maximum 
length. Let these terms lie in some set M{ 2°}. Then, the operating conditions of some output Z, including 
unused and indifferent states, can be presented as some subset of the set M, NC M. With this 


Fin =N = Dit [Diet Dhl. (7) 
where 
Dh= Fi. > fn = Fnand >} fp = Fp. 


*This condition for realizability was formulated some what differently in [16]. 
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It is obvious that any term ¢ of the normal form of a Boolean function will have a length 7 9 equal to or 
less that the length, I, of the terms of the complete normal form: 1 y= If. 


We shall say that a Boolean function in either its normal or its complete normal form, F(z) Doe or F(z, = 
Zuf + realizes the operating conditions of output Z if the totality of its terms realizes all the working states and 
realizes none of the forbidden states. 


We shall consider that a term of the normal form of a Boolean function has minimal length, or is "(minimal” 
(¢% min). if it contains such a number of letters that decreasing it by at least one letter (excluding the very last one) 
would lead to a condition where this term realizes at least one combination of states which is simultaneously in 
Fiz) and in F(z), 


Theorem 1. The minimal form of a Boolean function which realizes the operating conditions of one of the 
outputs Z of a switching circuit is that form which contains only minimal terms.* 


F2min= Di Pmin. (8) 


Indeed, if at least one of the terms ¢ has a length/ greater than the minimal length, then tl.is means that 
the length L of the expression F7, can be decreased. If at least one of the terms ¢ has a length less than the 
minimal length, then this term will simultaneously realize at least one of the states in Fj¢z) and Faz), i.e., the 
normal form of the Boolean function containing it will not, by our definition, by minimal. 


On the basis of Theorem 1, we can present the following algorithm for obtaining the minimal form of a 
Boolean function realizing the condition F;;7): we shall choose terms of the normal form which realize states of 
Fiz), starting with terms containing only one letter and, gradually increasing the number of letters, we shall 
verify, for each of these terms, whether or not they realize some state from Fz). The shortest terms ¢ which do 
not realize states in F(z) will be entered in the minimal form of the Boolean function which realizes Fiz.” * 


We note that each of these minimal forms will automatically include the necessary number of terms of the 
complete normal function from f;, fy and fp. 


We now consider the application of this algorithm to the example given above in Fig. 2, starting this con- 
sideration with the operating conditions for output Z;. If we consider for this the elements of fy, we immediately 
convince ourselves that c = 1 is absent from them. Therefore, the minimal term c can enter into the minimal 
form of the Boolean function which realizes Fy(z,). It is obvious that all terms containing this letter as a factor 
will not be minimal, and these, therefore, need not be tried in the remainder of the process. The two-letter 
combinations which do not give realizations of states in Fy) are the terms ad and bd. 


Thus, the minim] form of the Boolean function which realizes the operating conditions of output Z will 
have the form: 


Fizy =ce+ad+hd. (9) 


By determining the minimal form of the Boolean function which realizes the operating conditions for output 
Z, in the same manner, we obtain 


Fiz) = ac + ad + be + be + bd. (10) 


* A theorem analogous in its formulation was proven by W. Quine in [3]. However, in their essentials, the theorem 
stated above and Quine's theorem are different, and from them there flow different algorithms for obtaining 
minimal forms of Boolean functions. Quine's theorem is based on the concept of the so-called “prime implicants” 
which are obtained with the use of formulae (1) and (3). A special procedure is thereby necessitated for taking 
the unused states into account (Cf., [4]). Theorem 1, formulated above, is based on the concept of “minimal 
terms” which are obtained by the method, previously given by the author, of analyzing the conditions of realiz- 
ability of a switching table (Cf., [16]), which automaticaly takes into account the necessary unused states. 

** As applied to nonsequential circuits, this algorithm was first used by the author in the previously cited paper 
given at Harvard University. It was subsequently applied in [17] for the case of minimizing the number of diodes 
in structures with diode nets. 
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A minimal form which contains all the minimal terms will be called a “general” minimal form, Such 
general minimal forms are those given by expressions (9) and (10). In many cases, however, the conditions f; 
can be realized by means of some, not necessarily all, the minimal terms. Obviously, it is sufficient for this 
that there be enough such terms to realize all the conditions in fj without exception. We shall call such a minimal 
form of a Boolean function which realizes all the operating conditions of some output of a switchiiig circuit, and 
which contains the minimal necessary number of minimal terms, a “particular” minimal form. 


To determine, from a general minimal form, all the possible particular forms, it is convenient to use the 
methodology presented by McCluskey [6].° In this methodology, a table is set up for the general minimal form, 
wherein it is noted, for each of the minimal terms, which of the individual states of Fj,from the switching table,it 
realizes relative to the given output of the switching circuit. 
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We number the individual conditions for Fiz ,) and Fiz, as indicated on Fig. 2, we denote the minimal 
terms by capital letters and we denote by a cross the fact that an individual condition has been realized, this 
cross being placed at the intersection of the row representing thiscondition and the column representing the correspond- 
ing minima] term. Then, for the example considered above with the two outputs Z; and Z,, we obtain the tables 
shown in Fig. 3. It is obvious that, from the entire set of minimal terms contained in the general minimal form, 
it is necessary, for each of the partial minimal forms, to choose those of their combinations for which the condi- 
tions of Fj would be realized without exception. 


We write this in the form of formulae, wherein we use the logical connective “or” to denote the conditions 
for realizing the individual f; and the logical connective “and” to denote the condition, cited above, of the 
necessity of realizing all the f;. 


Then, by using the conventional notation for conjunction and disjunction, we obtain: 


for the table of Fig. 3,a: 
D,=P(QV P)R(Q\V R); (11) 
for the table of Fig. 3,b 
O,=(TVQ)(PVRVS)R. (11a) 
By using the well-known formula for simplifying Boolean functions 
af(z.y,...,w) =a2f(i,y,..., w) (12) 


and by transforming the formulae obtained to the normal disjunctive form, we shall have 
*This problem can also be solved by means of the method of closing sets, suggested by K. Voishvillo. 
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O,=P(@V P)R(QYV R) = PR, (13) 
O,=(TYVQ(PVRYSR=(TVQR=TRyY QR. (14) 
In correspondence with this, the minimal forms of the Boolean functions for the outputs Z; and Z, are written 
in the following way: 
Fiz) nin = 4 + 24, (15) 


= ac + be = be + be. (16) 


Fit) min 


The method just presented gives a general solution to the problem of minimizing the Boolean function 
characterizing a switching circuit, but requires the execution of a sorting operation which, in certain cases, can 
be quite awkward. 


If it is necessary to obtain only one of the partial minimal forms, then the sorting operation is not necessary. 
To determine the algorithm for obtaining such forms, we consider Theorem 2. 


Theorem 2. A Boolean function, expressed in the normal form F’ which has one or several terms ¢' of 
length 2 * less than the minimal length, i.e., which realize simultaneously, not only all the f;, but also certain 
fy, can be taken to a normal form which realizes only f, if, for the conditions f, mentioned above, one sets up 
one of the partial minimal forms y, and then takes the product F' ¥. 


A term of minimal length is formed from a certain set of terms of fj, fy and f,, of the complete normal 
form as the result of employing operation (1). If its length is less than the minimal length, this means that at 
least one f, enters into the set of terms f forming it. 


= Dalit digin + Dy fo + dislo: (17) 


To remove this contradiction, it suffices to multiply ¢' by such a function that the resulting product goes to 
zero when there is present any of the combinations of input states corresponding to each of the f, entering into 
py In the most general case, this can be the inverse of De he. 


However, in correspondence with the definitions given above, this multiplier can also be inverse of any 
normal form » of the Boolean function which satisfies the condition: y realizes 2. /,. 


Thus, 


e=9'9 = (DLA + Dp fn +d), fo + De fo) >. (18) 


These arguments are also valid for the entire function F', since F' = D9’. 
Thus, indeed, 
Fy= F'®. (19) 


The theorem is thereby proven. 


We note that by carrying out the operations corresponding to Theorem 2 we can obtain, in addition to the 
terms realizing f;, terms which only realize f,, and fp- One can convince oneself of this by multiplying out the 
parentheses in (18): 


9 =F Yuht F (Ypfn t+ Dy tr) - (20) 


The second term of this expression can give terms which do not enter into the set of minimal terms for f;, 
but which only realize some terms of f n and f p: Therefore, the solution obtained in this case must be compared 
with the conditions for f;, and those terms which do not realize the latter must be discarded. 
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On the basis of Theorem 2 we can advance an algorithm which allows one to obtain particular minimal forms 
without the use of sorting. 


It is necessary for this to take a minimal set of certain terms g* with lengths less than the minimal length 
(for example, containing only one letter) and which realize the conditions of fj without duplication, such that 
each of them realizes a minimal number of conditions of f,, and then carry out on these terms the operations 
cited above. 


Let : 
F’ = 9 +924 .--4+9n- (21) 


For each y* we determine the corresponding function » and then carry out the termwise multiplication of 
each ¢’ by the inverse, J: 


Fi = 01 1 + G2 Gat ++» + on On (22) 


The function obtained must be the minimal form of the Boolean function which realizes all the conditions 
of fj. 


In the general case, each of the » is a sum of minimal terms € realizing Fy. If any of these terms has a 
length greater than the minimal length, i.e., has more letters composing it, or if one of the expressions » has a 
superfluous term, i.e., » is not a particular minimal form, then the corresponding product, gy), will form a 
superfluous term of ¢ or a term of ¢ with a superfluous letter, i.e., which is not minimal. Thus, in both these 
cases, the Boolean function obtained will either not be a minimal form or not be a particular minimal form. 


Now, let a term € with a length less than minimal] enter into one of the expressions »,. The function p, 
will then realize, not only terms of £, but also some terms of f;. This means that the corresponding product, 
?k¥ Will go to zero for some combinations of input states, corresponding to f; and, since the g* were so chosen 
by us that they realize all the conditions fi without duplication, this then means that F; will not realize all 
conditions f;, i.e., will not be a minimal form of the Boolean function. 


We now illustrate the application of the algorithm considered here by the example of the switching table 
of Fig. 2. 


It can easily be seen for output Z, (Fig. 2,b) that all the f; can be realized by means of the term ¢' = d, 
whereby this term also realizes a minimal number of fy. We now set up the table of conditions of f; and f, 
realized by y’ = d (Fig. 4,a). It is easily seen that the minimal term which satisfies the condition that it only 
realizes Fy will be € = ab. Thus, » = ab and 


F, = 9'$ = d (ab) = ad + bd, 


which corresponds with the result obtained previously. 
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We now take, as the ¢’, the terms a and 5, each of which realizes two of the fy. The corresponding tables 
for conditions f; and fy are shown in Fig. 4,b and c. As may easily be noted, for b the conditions of f, are 
realized by the terms £ , equal to d and a, and for a_ these conditions are realized by the € equal to b and d. 
Therefore, 





95 = b(d +-a) = abd and gq = a(b + d) = abd. 
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The terms obtained are not minimal, since b and a do not satisfy the conditions formulated above. 


If we now choose, as the terms ¢', the terms @ and b, each of which realizes only one fy (Fig. 4, d and e), 
we obtain 


; = ad = ad and » = bd = bd. 


In this case, since the terms ¢‘' were chosen in correspondence with the conditions formulated for the Boolean 
function realizing f;, they are minimal. 


For the output Z,, the conditions formulated for the ¢" are satisfied by terms b and c. By analyzing the cor- 
responding tables of conditions f; and f, (Fig. 4, f and g), we obtain 


P = b(c) = be ani, = ¢ (b) = eb. 

The algorithm being considered is not universal, since the results obtained with its use depend heavily on 
the proper choice of the terms g’. However, even when their choice is not completely satisfactory, the algorithm 
significantly speeds up the process of determining the particular minimal forms, since the solution obtained is 
close to the minima! solution and can always be quickly taken to the minimal solution by testing each of the terms 
obtained by means of the first algorithm presented here. 
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Fig. 5. 


We note that the algorithms which derive from Theorems 1 and 2 can be applied, not only when the con- 
ditions F; and Fy are given in the form of a switching table, but also even when they are expressed by Boolean 
functions in normal form. On the basis of these latter formulae, one can set up tables of binary numbers in which 
the letters lacking in the corresponding terms of the norma] forms are noted by means of lines. After this the 
determination of the minimal terms is carried out by means of the algorithms described above, but with it being 
taken into account that a line can mean either a zero or a one. 


We now turn to the consideration of sequential switching circuits. These circuits differ from the non- 
sequential ones in that they contain connections between the output states and the input states. This property is 
@ consequence of the fact that for part (or for all) of the output there is a feedback path to the inputs (Fig, 5,a). 


It is convenient, in this case, to present the switching table in the form of a sequence of combinations of 
circuit elements’ states, marking the moment of switching in by a plus sign, and the moment of switching out by 
a minus sign [15, 16]. An example of such a table is shown in Fig. 6 (in its number of inputs and outputs, this 
table corresponds to Fig. 5). 


If the feedback path can meaningfully be sundered, a sequential circuit can then be transformed to a non- 
sequential one (Fig, 5,b), and its operating conditions can be described in the form of the so-called nonsequential 
equivalent (Cf. the table in Fig. 7).* 


This circumstance permits the reduction of the problem of minimizing Boolean functions for sequential cir- 
cuits to the same problem for nonsequential circuits, and allows the methods presented earlier to be used. How- 
ever, the presence of a definite temporal sequence in which the circuit's input and output states follow each other, 
allows, in this case, a significant reduction in, or complete avoiding of, the sorting operation. 





* This concept was introduced by P. P. Parkhomenko. 
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As was shown in [16, 18], in the Boolean functions characterizing each of the outputs of a sequential cir- 
cuit, there must be contained, as a minimum, the terms characterizing the elements which change state when 
the circuit of a given output must be closed (in a sequential switching table, at the cycles designed for the 
switching in of the elements) and when the circuit must be opened (at the cycles designed for the switching off 
of the elements). If the switching table which contains the sequence of switching on and off of these elements is 
realizable, then the Boolean function for the given output can not contain terms which correspond to any other 
element. If the switching table for the aforementioned elements is not realizable, it is then necessary to add to 
it, from the original switching table, a minimal number of sequences of these elements which would make it 
realizable. This operation still does not give a minimal form of the Boolean function, but does permit the 
immediate discarding of the letters which will not enter into all its terms which, certainly, significantly simplifies 
the minimization problem.* 


If we carry out this operation on the table of Fig. 6, then, for the individual elements, we obtain the tables 
given in Fig. 8,a,b,c,d ande. As may be seen from Fig. 8, the switching table obtained is realizable only for 
element Z,. Since one sequence of switching of element Z, occurs in this table, the minimal Boolean function for 
Z, will then have the form: 


Fizanin = 4 


The remaining switching tables can be made realizable if sequences of actions of the very elements which 
comprise them are added to them (Fig. 9). 


The method just considered thus allows one to decompose an ordinary switching table for a sequential cir- 
cuit into a number of sub-tables with fewer numbers of elements. The minimization problem is thereby signi- 
ficantly simplified. 


However, the special features of sequential circuits allow an even greater simplification of the minimization 
problem to be made. 


As was shown in [16], the transition from a sequential switching table to the Boolean function realizing con- 
dition fj can be carried out by means of the formula 


Fuz) = fop + 2Fi, (23) 


where f,. is a sum of the terms which correspond to the combinations of circuit elements states in the first cycle 
of each of the switched-in periods of element Z, but without the state of this element itself being taken into 
account, and F,” is the sum of the terms which correspond to combinations of circuit element states in the other 
cycles wherein the Z circuit must be closed, again with the state of element Z being left out of account. 


We shall call f.. the operated state, and zF;" the switched-in state. 
op pe i 


Theorem 3. A particular minimal form of the Boolean function which realizes the switched~in states 
accurately with the possible exception of the letter z may be obtained by choosing such a minimal number and 
such factors which, in their totality, would overlap all the cycles of the switched-in state, exactly coinciding 
with the succeeding cycles following each of the switched-in periods of element Z. 


Indeed, all the terms forming zF;" have the factor z which equals one only for fj. Therefore, this factor 
in conjunction with any other term ¢’ which satisfies the conditions given above, gives aterm g. It must thus be 
proved that this term will be minimal, or differ from the minimal by the letter z. 


If the element Z, in the switching table, has one operated period, i.e., is switched in and switched out once, 
it then suffices, for obtaining overlap, to choose a term ¢‘ which contains one letter. It is obvious that an exact 
coincidence with the following cycles of F;" can be guaranteed if, as one of the ¢', a letter is selected which cor- 
responds to an element which changes its state after the circuit of element Z is opened. Since an exact coin- 
cidence with the beginning of F," is not required, it is then always possible to find such letters which, with a 
definite overlap, would cover the remaining cycles of Fj" without redundant repetitions. 





* An analogous operation can also be carried out on switching tables for nonsequential circuits: in them can be 
stricken on individual columns, thus decreasing the number of letters involved in the sorting, until the point is 
reached when at least one f; coinciding with an f, is obtained. 
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Fig. 10. With these conditions: 


a) none of the terms ¢' can be eliminated, since 
this calls for going counter to the condition of overlapping 
all cycles of Fj", i.e., at least one of the conditions fj entering into Fy" would be unimplemented; thus, the 
number of terms in this case is minimal, i.e., corresponds to one of the particular minimal forms for zF;”. 


b) the terms ¢ all turn out to consist of two letters, one of which will be the letter z, i.e., will in fact 
differ from the minimal by no more than z. 


For example, we consider the table for element Zs in the example given above (Fig. 9,d). 


The term fop. in correspondence with the element states at cycle 3, will equal 


hop = ZalZy. 
Switching off of element Zs occurs after operation of element Zs. Therefore, to provide exact coincidence 


with the cycles following the switched-in state, it is necessary to choose term Zs as one of the ¢’. To overlap all 
the switched-in state cycles, it is necessary to add ¢y" = z, to this. Thus, we shall have . 


F(z, = Za23 + 25 (22 + 2s). 
If we set up the enumeration of conditions f; and fy (Fig. 10,a) for the table of Fig. 9,d and then deter- 
mine the minimal implicants by means of the first algorithm, we shall obtain the same results, 


Now, let element Z have several operated periods. In theory, it is possible to overlap all the switched-in 
states in the same manner as before. It may turn out, however, that a term ¢’, consisting of one letter and 
providing exact coincidence with the succeeding cycles of the switched-in state may, in different operated 
periods of element Z, give rise to contradictory conditions. Thus, for example, in the example considered 
above for element Z, (Fig. 9,a), the term which provides coincidence with the succeeding cycles in the first 
period is the term ¢;' = zg and, in the second period, the term 9} = zg. However, in this same period, element 
zg is switched off at the limits of F;" and, therefore, the term zz will disrupt the given sequence of actions of 
Z,, retaining the value zF;" = 1 even at the limits of F;". 


This is indicative of the fact that, in this case, the term ¢', which is one letter long, has a length less 
that minimal. For such terms we seek a set of letters, of minimal content, which would liquidate the contra- 
diction noted. 


The expressions thus obtained, realizing the switched-in state, will have in each terma minimal number 
of letters necessary for the conditions of realizing the switched~-in state, and a minimal number of terms them- 


selves, and will also have the letter z making them up, i.e., will differ from the possible minimal terms by not 
more than z. 
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Thus, for example, for Z;', in correspondence with the switching table of Fig. 9,a, we will have: 
for the first operated period 
hop = 292425 and Fi. = % (25 + 24); 
and for the second operated period* 
hp = 23242, amd 2; = 2,2,. 


Since, as was stated, the term ¢' = zs contradicts the succeeding cycles of the switched-in state for the 
second operated period, we add to it an additional letter which would differentiate the first operated period from 
the second. Such a letter, as may be seen from the table of Fig. 9,a, is zs. Thus, 


F(z.) = 4232425 + OZ 5242, + 2, (2523 + Ze + %) nae (24) 
= GZ42Z, + 2, (24 + 2325). 


If, in constructing the table of conditions for fj and f, for Z;, we employ the first algorithm considered 
(Fig. 10,b), we obtain a general minimal form of the Boolean function in the form : 


Fz.) = 2424 + 242925 + 2y24Zy + 242925 + 02425. (25) 
By using McCluskey's method (cf, the table in Fig. 11), we obtain 
M@=E(CV AAV DYE )A(AV CV E)(A\Y B)B= ABE. 
Thus, the general minimal form has the form: 
Fiz, = 224 + 24252, +- A242Z,, 


i.e., coincides completely with the result obtained above. 


To obtain the general minimal form, it is necessary to take, for each of the outputs, all the variants of the 
switching table to be realized which have a minimal number of elements and, for each of these tables, to take 
all the variants of overlapping the cycles of F;”. 


We considez one more example of the sequential-nonsequential circuit. Let the switching condition of 
element Z be the application of stimuli at inputs A, B, C and D in the following order: +A, +B,—A, +C, —B, +A, 
and +D. Any disruption of this sequence must cause element Y to switch in. For element Z the operating condi- 
tions are sequential, since it must be switched in after a definite unambiguously given temporal sequency of 
stimuli applications to inputs. For element Y these conditions are sequential-nonsequential, since this element 
must not operate for the sequence of stimuli given for Z, and must operate for all other sequences, whose temporal 
orders are not given. 


If we consider each cycle of the correct sequence of application of stimuli, and after these if we enumerate 
all the possible incorrect stimuli, for which element Y must operate, we obtain the table of conditions in f; and 
fy which is given in Fig. 12 (the unprimed cycle numbers denote cycles of the given sequence while the primed 
cycle numbers denote cycles in which incorrect stimuli can occur). 


The complete switching table, with weighed states, for Z and Y is given in Fig. 13,a. For element Y, this 
table is unrealizable. An analysis of the conditions in it [15, 16] shows that three additional elements must be 
included in it in order to make it realizable. The switching table with these additional elements included is 
given in Fig. 13,b. 


Element Z is switched in after a stimulus is applied to input D. However, to obtain a realizable switching 
table, one such element turns out to be insufficient. For this it is necessary to introduce a sequence of elements 
A and C in the table (Fig. 13,c). 


* The third period (cycles 19 and 20) is a repetition of the first. 
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Fig. 14. 
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For element Y we apply the algorithm flowing from Theorem 1, since its operating conditions are non- 
sequential. The table of conditions f; and f, for it are given in Fig. 14. By carrying out analysis of the minimal 
terms, we obtain the general minimal form in the form 


Fiy) = da + db + de + abe + 24d + z,ba + 2,ba + 
+ ca + x,cb + rgba + tycb + Lyx, + 7y7,b + 
+-xgba + x,ca + 24cb + 22,0 + 22,5 + xo2,¢ + 
+ aqd + rgd + xyca + rgcb + ryxyc + FyrQb + 
+ Tyee + Lge + aba + tyra + TgX,b + X3%_a + Taq). 


The table for realizing the conditions f; according to McCluskey are given in Fig. 15. The formula for 
overlapping the conditions f; deriving from it will have the form: 


=(, VI) L(AVCV PV LV £5) (Ci V Is) (La V Ls V 
VXVE)(CV L,) (Ks V Ly) (Fs V Ky (AV BVCV FLV 
V K,V4)(Ki V AsV Ky V L,)(DV G3 V Ka V Ke V Lal Lu) X 
x (BY CY K,VL)MiUi V In V Ia VV My) (AV Fi) MX 
x (GV G.V 13 V M;)(DV G,V Ie VV Mz VV My VV My) M, x 
X(DV F.V Ki V Kau V MaV M,)(AV FV K;). 


By transforming this formula, we get 


@ = 1,M,(AV Fy) (i V 14) (Gi V Js) (C V Ly) (Ks V Ls) x 
x(F, V K,)(IDV Gs V (Ka V Ko V Lo V La) (la VV Ma VV Ma V M,)). 


We can obtain all the particular minimal forms from this formula. One of them, containing the least 
number of terms, is obtained if we take the conjunction consisting of the first terms of the bracketed expressions 
in the logical formula, i.e., LgMyAlyGyCKgF,D. The minimal form of the Boolean function which corresponds to 
this will have the form: 


Fy) = 2p %¢ + 23¢ + da + 2,ba + x, ba + de + 
+ yb + 2, ba + abe = d(a + ¢) + ab(c + ,)+ 
+ ¢(23 + 2,6) + 2, ba + Xg2qC. 


Thus, the minimization methods considered above allow a rapid determination to be made of the minimal 
forms of Boolean functions for quite complicated conditions. 


Received August 12, 1958 
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POSSIBILITIES OF INCREASING NOISE STABILITY BASED ON 


_ 


THE USE OF A PRIORI PARAMETER PROBABILITIES 


V.A. Kashirin and G. A. Shastova 
(Moscow) 


The noise stability of an ideal receiver is considered, the first dis- 
tribution law of the a priori parameter value probabilities being taken 
into account. An analysis of noise stability is carried out for the case of 
a nonlinear parameter transformation prior to modulation. The optimal 
nonlinear transformation of parameters prior to transmission is found, as 
well as the gain in noise stability obtained with this transformation. Both 
the mean square and the information criteria of optimality are used. 


It is shown in the theory of potential noise stability [1] that the receiver which identifies its received 
signals with their most probable parameter values is ideal for weak fluctuating noise. The noise stability of 
such a receiver was investigated in detail by V. A. Kotel'nikov for the case of equiprobable transmission of all 
parameter values in some range. In practice, the parameter probability distribution law is frequently unknown. 
In this case, the receiver is generally chosen on the basis that all parameter values have equal probability of 
being transmitted. It is advantageous to consider the noise stability of an ideal receiver, designed for a known 
first probability distribution law of the parameters, this law differing from a uniform distribution, and to compare 
its noise stability with the noise stability of a receiver constructed for the case of an unknown first distribution 
law. This allows us to judge what a knowledge of the distribution law can give us in terms of increasing the 
noise stability of reception. 


Noise Stability of an Ideal Kotel'nikov Receiver for a Nonuniform Parameter 
Probability Distribution 








Let the parameter A, varying with a certain probability distribution P(A ) within the limits of +1, be trans- 
mitted by means of a signal A(A, t) along a channel with weak fluctuating noise of intensity 0, in cps~*/*, By 
using the methods of potential noise stability theory, we find the probability that, if the signal x(t) is received, the 
parameter A lies within the limits A" < A < A‘ + dA. This probability equals 





exp {- Ls [a (t)—A(d’, t)]?+ InP on} dr 
PIV <i <V4+a= 7 = P,(d’)dr. 


\exp -- 4 [z(t) — A (A, t)]? + In P ov dd 
-1 








Here, T is the interval during which the parameter value is taken to be unchanged and 





T 
(z(t) — AQ, OP — ¥\ [x(t) —A(QA, sede. 
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The function P,@') characterizes the probability of various errors in the reproduction of the parameter 
when x(t) is received. It can be written in the form 





Ps(’) = keoxp {— 3 [2 Q— A, OF + In Par}. 


Here, ky is a coefficient which depends on x but not on A‘ or t. For a given x(t), the most probable value 
of the parameter A,,, must minimize the following expression 





{[2@—4Q, OF—inP@)} = min. 1) 


xn 


With weak noise the errors on the receiving end are small and it is therefore possible to consider only the 
probabilities of reproducing parameter values close to the most probable value. This allows us to discard the 
higher terms in the expansion of function A(A, t) in the neighborhood of A,,,, i.e., to consider that 


A(h, t) = A(han, t) + Ar Qrens #) (OK — Ran): (2) 


Here, 


; 0A (A... t)- 


xn 


By using (2), as well as the expression 


( P.Q’)d’ =A, (3) 


—oc 


we can define k, and, consequently, the function P(A"). 


The magnitude of the mean square error equals 


B= | den) POR)’. (4) 


It should be noticed that, when (3) and (4) are integrated between infinite limits, an error is admitted but, 
for weak noise, P,(A") falls so rapidly with increasing (A‘-A,,,) that the integration between infinite limits does 
not entail any essential error. Solving the problem posed in its general form gives rise to mathematical dif- 
ficulties. We consider the particular case wherein the parameter is distributed normally, symmetrically about 


zero: 
»” 


2 
PQ’) =ge "° A>a>—t. (5) 





The probability that, if x is received, parameter A" was transmitted will equal 





J , ’ Ken ’ x" 
Ps (’) rs ks — [- a [Ax (hens t) (A or han)? Poe » (A al hen) —_ =} . (6) 
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From condition (1) we easily obtain 


27 Ren 
at ltt) — Ahem Ar Qan, ) = 8 
f 








By substituting (6), in (3) and determining k,’, we find the error probability function: 


P, (X’) = eN V2 C—O Ayn) PN | 
4 
where 


2 
xn 


= ne 
2,2 ” 


an 


r= 
p= “ou A’ (hen; t) + 
We find, from formula (4), the magnitude of the mean square error: 


a? “ei Us 
coy eee Or 3 
0 4B 





Here, 5, is the mean square error given by the ideal receiver when the signal statistics are not taken into 
account; 


ee ici cee : 
8T A? (An, ¢) ” 





When the statistics are taken into account, the error is decreased by a factor of By: 


B *o 1+4 *o 
cats te + re (8) 


It follows from (8) that an appreciable gain in noise stability is obtained only for 4) < 459, i.e., for very 
small Ag», since 5, is small in the case of weak noise. From the practical point of view, the case Ay < 45, is 
rarely realized, Ordinarily, Ay >> 455, and then dy * 5p. 


Thus, in the majority of practical cases, the ideal Kotel'nikov receiver provides no increase in noise 
stability when the parameter's first distribution law is taken into account. An essential disadvantage of a receiver 
which takes the first distribution law into account is that its operating mode depends on the ratio of the spectral 
density of the noise to the parameter's dispersion. For some definite ratio, o7/)2 , this receiver will provide 
minimum error, but if this ratio is changed, the receiver will not provide this minimum. In particular, if the 
receiver is designed for the reception of signals with a definite noise level, it will then give an error when the 
transmission is noiseless, this error being of the same order of magnitude as the gain which the receiver provides 
when the statistics are taken into account as compared to a receiver which does not take the statistics into account. 
We shall illustrate this with one example. 


Let the signal be transmitted by means of amplitude modulated dc, i.e., A(A, t) = UA, and let the para- 


meter obey a normal law, i.e., 
a2 


2n2 


PQ\)=qe °. 
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The ideal receiver reproduces the most probable parameter value which satisfies condition (1). By computing 
the defivative of P(A) with respect to A for A = Axn, We obtain the equation defining A xn? 








2TU » 
= ~ on [x (t) — Um dan] — — ==: (). 
x 
0 
From whence, 
2 | x (t) 
hen 4 o } pos 
M2 2TU2, 
By taking into account that 
o2 
«tates adi 9 
sTu3, ss ” 


we obtain the following expression for the znost probable value of the parameter A for the noise level defined by 
equation (9): 


ee EEN 


If a signal without noise acts on this receiver, i.e., x(t) = U,,Atrans, it then turns out that the received 
value, Aree, will differ from the transmitted value, A+,ans: 


For Aq >> 459, the error of reproduction will equal the ratio (46 o/ ro}. i.e., the gain given by the receiver 
when the statistics are taken into account. 


The analysis given shows that, in the case when the dispersion of the transmitted random parameter is many 
times larger than the admissible magnitude of the mean square error, the use of the knowledge of the first prob- 
ability distribution law in the ideal Kotel'nikov receiver provides a negligible gain in noise stability and gives 
rise to an additional error when the noise level changes. 


Nonlinear Parameter Transformation 








We now consider the influence of the method of modulation on noise stability from the point of view of 
using signal statistics. When a parameter A is transmitted through a channel with fluctuating noise to an ideal 
receiver, the mean square error is defined, as is well-known [1], by expression (7). 


The mean square error, averaged over-all values of parameter will be 





jo \ Pi) dh. (10) 


1 [~ “ t) y dt 


1 





ee | 
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It follows from (10) that the mean square error depends on both the function P(A) and the function A(A, t). 
Obviously, for each given function P(A) one can find the optimal function A(A, t) which would minimize the 
mean square error. In the majority of modulation systems, the function A(A, t) depends on d in such wise that 


([24a-or dt = const. 


0 


In this case, the error is completely determined by the function A(A, t), and does not depend on P(A). If, 
however, there exists a dependency 





SRT a7 0 


0 


then, according to (10), the resulting error will depend on the parameter's probability distribution [the function 
P(A]. There now arises the problem of finding the optimal function A(A, t) for a given P(A). 


@A (i, t) 


2 
aL ] dt =f(h) can be obtained either by the use of a special “nonlinear” 


T 
The dependence \ [ 
0 


modulation or by a nonlinear parameter transformation prior to modulation. 
We shall consider the second method in more detail. 


In the transmitting device, before it is applied to the modulator, let the parameter A be transformed to 
some auxiliary parameter y = y(A). For transmitting the parameter y along the channel, such a form of modulation 
is used for which @A(y, t)/8y does not depend on y, and the error in transmitting parameter y equals 5. 


Then, by taking into account that 9A(y, t)/ 8A =[8ACy, t)/ dy] [8y/ 4A], we find that the mean square 
error, averaged over all A, equals 








1 1 
7 , o? [ro ar =8 (rae 
a (| 24 ay — “1 
dy 
0 
where 
, _ Oy 
tio 7 





32 4 
eas ree te Sa 
To P() (11) 
7y adh 
\w 


The influence of the nonlinear transformation can also be explained without the ideal receiver. Figure 1 
shows the nonlinear function y = y(A) and the probability density P(A). Let a quite small error Ay, constant over 
the entire range of values of parameter y, arise when the parameter y is transmitted. 
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ation 


(11) 


1 
over 





For small A, the errors Ay and Ad are related by the relationship 
lA): P(A) Ay = Ady’. 


The mean square error of reproduction of parameter A equals 
‘y 


eS fee Gain? eek . 


y? 


‘ 
' 
1 
i 


a | 


Since Ay = 5, and does not depend on A, we then oltain the 
same result as in (11) for the ratio of the errors. 


= 


Obtaining the greatest gain in noise stability by taking statistics 
é into account thus leads to the search for the optimal function y = f(A) 
+1 a which would minimize the integral 








rf-----------> 


S 
o's oo 
= 


1 
P(r 
Fig. i ;e™ \ aa dh. 


1 
The problem thus posed is equivalent to the variational problem of finding the minimum of the integral 


Xe 


I (y) =\F (h, y, y’) dr. 


x 


Here, F is a continuous function, differentiable with respect to y' for all values of y’ and differentiable with 
respect to the arguments A and y in some region ofthe OAy plane, and y = y(A) is a function which is continuously 


differentiable on the interval A; < A < Ag. In our case, the function F(A, y,y") = P(A)/ y*(A) satisfies all the 
conditions enumerated above. 


The solution of such a variational problem leads to the differential equation 
—-_ Fy =0, (12) 


where 


In the case considered, Fy = 0, Fy" = —2P(A) Cy"), and differential equation (12) has the form: 


oF OD (y’)? — 3P (d) (y')ty” =0. 





By assuming that y’ © 0, which always holds in our practical problem, and by integrating the separated 
variables, we get 


PQ) d+ ¢. (13) 


I 

Oo 
° 

5 
“_o-~ 
> 
at 
<= 

I 

is) 
) 

cs» 

i 
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The coefficients cy and c, are found from the conditions that y = 41 for A = 41: 
2 


amet ae ; , Cc; => — 1. (14) 
\V Pia) da 
—1 


Co = 


By substituting (13) and (14) in (11), we find the expression for the gain in noise stability: 


oe 2 


1 


[ p*/s (n) a}" 3 


ine" 





It should be mentioned that an analogous solution is obtained when one seeks the optimal nonlinear quan- 
tizing of the signal [2]. 


We now consider two particular cases: a truncated normal probability distribution and a step-function para- 
meter distribution. 


The Truncated Normal Probability Distribution 
Let P(x) satisfy equation (5), where q is defined by the condition 





Then, 


x» 


3 
y=cVge “°. 


The derivative of the optimal nonlinear function has the form of a normal law. The nonlinear transforma- 
tion function sought will be the integral of the probability 


The gain in noise stability equals 


Bu, — 0-475 @'/ (4 | Ao) 
N ho @M*(4/2,)’ 





where 





@ (z) = 7 Je 2 du. 


Figure 2 gives the gain as a function of 1/Ay. For 


<1 By=0,35 = 


1214 





t) 


ra- 





It is clear from the figure that an essential gain in 
noise stability (greater than a factor of 2) is only obtained 
for 1/A5 > 6, i.e., in the case when the parameter dispersion 


3 Z is many times smaller than the entire range to be trans- 
V4 mitted. 











ae ~~ The Step-Function Probability Distribution 
aed 


¥ The step~function probability distribution (Fig. 3) 
may be written in the form 

















Ro for AIS 


BERS 





P(Q)=) » 
¥ for 1>|A|> 
7" 2 2 8 See Ue lhe 


(15) 











f 
@ 
f 
a 





























S 





If the parameters a and g are given, then the coef- 
Fig. 2. ficient R can be determined from the condition 




















P1A,)-y(A) ; 
P Pia) from whence 
Ss > ap an aa 
. ee FS: 
R= sayER7): 
1 
; In accordance with (13) and (14), 
aVi 
‘ y, => for |A| <1, 
VB +a-1 rs 
Rk - , S 1>|raj>+ 
y , haath | Yu = for a 
0 Va a VB+a—1 


Pig. 3. The optimal nonlinear transformation for a step- 


function probability distribution is thus a piecewise linear 
transformation with breaks at the points A = +1/a, where 
the slope of the line y = f(A) in the central portion is * Vg times larger than the slopes of the extreme portions. 


The gain resulting from such a transformation equals 


( — i) 
By =a ps ~- 
(yv5+a—1) 


For a >> 1 and g >>a’, B,, = a, i.e., the gain equals the ratio of the length of the entire range of the para- 
meter to be transmitted to the length of the range of the probable parameter values. 





(16) 


It is easily shown that the ratio of the mean square error when transmitting parameter values in the range 
|A| > 1/c to the error in the range | A] < 1/a equals 


2 3 
fa oo" VE. (17) 
Ay Yi 
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Example. Let a = 5 and g = 1000. Then, the ratio of the errors equals Ay/A, = § 1000 = 10 and the gain 
in noise stability equals B,, = 3.5. 


Thus, with a step-function probability distribution, the gain in noise stability obtained is the larger, the 
larger is the ratio of the probability densities in both ranges of the parameter and the narrower is the range of 
probable parameter values. It follows from (16) that the gain cannot be greater than the ratio of the length of 
the entire parameter range tothe length of the range of its most probable values. It is clear from (17) that the 
gain is achieved by decreasing the error in the range of probable parameter values and increasing the error when 
low-probability parameter values are transmitted. 


The case considered above of a step-function probability distribution can be of practical value in systems 
of telecommunications of parameters if there is a certain range of parameter values which corresponds to normal 
operation of the object and there is a range of abnormal operation where the accuracy requirements in this range 
are quite low. 


A Nonlinear Transformation Optimal from the Point of View of Obtaining the Great- 
est Quantity of Information 








We previously considered the question of a nonlinear transformation optimal from the point of view of 
minimizing the mean square error. We now consider this same problem from the information-theoretical point 
of view. 


In a channel with noise, let there be transmitted some auxiliary parameter y(A), and let A* be obtained by 
the inverse nonlinear transformation of the received signal y,...(), distorted by the noise. For greater clarity of 
the results, we will assume that the parameter A varies within the limits [- a/2 to +a/2]. The quantity of infor- 
mation in A* about A equals 


J=H(’)—H, (0). 


Hence, 
a/2 
H(’) = — \ P(X’) log P(X’) dh’ 
—a/2 
and 
a a/2 
Hy, (’) =— i P (A) dah \ Py (X’) log P, (d’) dy’. (18) 

—a/2 —a/2 


For weak noise, we can assume that the error in the reproduced parameter is distributed normally. Due to 
the nonlinear transformation, the dispersion of this error is not constant over the range of values of A: 


de 
y’ (A)° 





8 Q) = 


With this, 


1 _(@-.p 


P) (’) —Vmne @. (19) 


By substituting (19) in (18) and by taking into account that the limits of the second integral in (18) can, without 
significant error, be replaced by infinity, we obtain 


7 
H, (X’) =logV2ne— \ PQ) logy’) dr. 


—a/2 
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18) 


to 


(19) 


hout 





The maximum amount of information will correspond to the minimum value of the quantity H)(A"). By 
solving the problem posed using the variational method described above, we find that the minimum occurs for 


y’ (4) = aP (a). (20) 
The maximum amount of information is 
a/2 a/2 
Jmax = — \ P (d’) log P(X’) dX’ + P (d) log P (4) dk + log V 2ne t : 
—a/2 —a/2 


For weak fluctuating noise, ignoring marginal values, we can assume that P(A) © P(A"), Then, Imex = 
= log V 2we(a/5,). If the nonlinear transformation is not present, i.e., y"(A) = 1, then the quantity of infor- 
mation equals 


J=— ( P (d) log P (4) dk — log V 2ne 8,. 


—a/2 


The optimal nonlinear transformation thus allows one to transmit a quantity of information which is greater 
by the quantity 


AJ = loga + ( P (4) log P (A) dh. 
—a/2 


It follows from (20) that the transformation which is optimal from the point of view of quantity of infor- 
mation is that which provides a uniform probability distribution of the auxiliary parameter transmitted along the 
channel with noise. This transformation is different than in the case when the criterion was the mean square error, 
where it is nonlinear to a greater degree for the same law, P(A). This corresponds physically to the fact that, in 
the case of the transformation which is optimal from the information criterion, the rarer values of the transmitted 
parameter are transmitted with a larger error than in the case when the criterion was the mean square error. 


Indeed, with the information criterion, the ratio of the errors for different 4 is expressed by the formula 


(s4] _ P (As) 
A (Xs) Jing P(A)’ 


and for the mean square criterion 





[4031 VPU) 
nh 


The possible increase in the quantity of information for the truncated normal distribution ['cf, (5)] equals 








b 
aa ite: ey ee aoa 
ASen = oe ou) + Tiwi phe” o6 Ve. 
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The graph of this function is shown in Fig. 2. 
For 1/A, > 3, 


>3 AJin = 10g0.49 (;.)- 


=|- 


In the case of the step-function distribution of probabilities (cf, (15)], the quantity of information can be 
increased via a nonlinear transformation by the amount 


—1 
Al, = — log? +P + span loss. 





For g>>a and a >> L, we will have AJ,, = log a. 


It is obvious that, for the nonlinear transformation which is optimal from the information criterion point 
of view, the mean square error will be greater than for the transformation which is optimal from the point of 
view of minimum mean square error. As computations show, this increase can be very significant. For example, 
in the case of a step-function distribution with g = 1000 and a= 10, the quantity of information can be increased 
by log,10. With this, the mean square error is increased by a factor of 100. This is due to the fact that the low- 
probability parameter values will be transmitted, with such a transformation, with an error 1000 times greater 
than for the most probable value. In practice, the choice of one criterion or another depends on the require- 
ments imposed on the system for transmitting data. In certain cases both criteria considered here might be 
invalid, since the transmission of low-probability values with low reliability might be inadmissible in some 
operating conditions. 


SUMMARY 


1. If the relative dispersion of the parameter is much larger than the required accuracy, the use of the 
known first distribution law in constructing an ideal Kotel'nikov receiver provides little gain in noise stability. 
An essential disadvantage of the ideal receiver which takes the first probability law into account is the dependence 
of its mode of operation on the ratio of the noise spectral density to the parameter dispersion. 


2. For any parameter probability distribution there may be found an optimal nonlinear transformation 
which minimizes the mean square error when the transmission is along a channel with weak fluctuating noise. 
In the cases where there is essential nonuniformity, the nonlinear transformations permit significant decreases 
in the mean square error. 


3. To transmit the maximum quantity of information along a channel with weak fluctuating noise, the 
. initial parameter probability distribution must be transformed to a uniform one. This means that the less 
frequent values will be transmitted with larger errors than the frequent ones. With this, the mean square error 
increases, 
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ANALYSIS OF ROYER MULTIVIBRATOR OPERATION IN 
INDUSTRIAL TELEMETRY DEVICES 


A. M. Pshenichnikov 


(Moscow) 


The conditions are investigated under which multivibrator oscil- 
lations are generated, the frequency of these oscillations being deter- 
mined by the magnetic polarity reversal time of the core. The functional 
dependencies of the circuit's output frequency and voltage and its power 
requirements on the input voltage, the circuit parameters and the core 
material and transistor characteristics are determined. 


A multivibrator circuit proposed by Royer [1] whose frequency is determined by the core's magnetic 
polarity reversal time is shown in Fig. 1. This circuit has been widely used in telemetering and analog computer 


devices as well as in power supplies for transforming dc to ac. Papers have recently appeared treating the use of 
Royer multivibrators in power supplies [2, 3]. 


In the present paper we investigate this circuit from the point 
of view of its use in telemetry devices. We consider the conditions 
under which oscillations are generated and we determine the circuit's 
output frequency and voltage and its power requirements as functions 
of the circuit parameters, the transistor characteristics and the core 
material characteristics. 





For our analysis we approximate the hysteresis loop by straight- 
line seqments (Fig. 2). The magnetic permeability equals jp» on 
segments 1-2 and 3-4 and, on segments 2-3 and 4-1, is assumed 
initially to be » = oo. In the sequel we shall determine the effect of 
a finite value of 1 on multivibrator operation. 











When a dc voltage is applied to the circuit's input, the collector 
current of one of the transistors (PT -1, say), due to the nonidentity of 
the transistor characteristics, will be greater than that in the other 
(PT-2). Thanks to this, the voltage induced in winding W, of transformer 
? T will be so directed that minus is applied to the base of PT-1 and plus 

to the base of PT-2. This leads to a further increase in the collector 
4 3 current of PT-1 and a decrease in collector current in PT-2. As a result, 
PT-1 conducts completely and PT~-2 is cut off. With this, the trans- 
former core initially changes its magnetic polarity along arm a-2, and 
0 4 then along arm 2-3 (Fig. 2). The emf induced in winding W, falls 
sharply when point 3 on the magnetization curve is reached, The cur- 
/ 2 rent in transistor PT-1 then decreases and, in winding W,, an emf will 
~ Fin |? be induced to the opposite sign. As the result, PT -1 is cut off and 
PT-2 conducts. Reversal of core magnetic polarity occurs along arms 
Fig. 2. The core's magneti- 3-4 and 4-1. Thereafter the circuit operates analogously. The multi- 
zation curve. vibrator frequency is determined by the time taken to reverse magnetic 
polarity along arms 2-3 and 4-1, since p >> pio. 


Fig. 1. A Royer multivibrator circuit. 
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I. Conditions for Oscillation Generation 





Until oscillations are generated, the core's magnetic state is defined by point a on the hysteresis loop (Cf., 
Fig. 2). In order to simplify the analysis of the conditions for circuit excitation, we shall assume that, with fluctua- 
ting oscillations, current i, (Cf., Fig. 1) increases while current i,’ remains equal to zero. Thus, the problem reduces 
to the determination of the conditions for self-excitation of an amplifier with positive feedback, the equivalent 
circuit of which is shown in Fig. 3. In this circuit, re, rp and r, denote, respectively, the equivalent impedances 
of the emitter junction, the base and the collector and r,,, denotes the equivalent junction impedance in the circuit 
with grounded base. The quantity a *1,,,/r, [4]. The equations for the circuit of Fig. 3 can be written as follows: 


O = in (re + ro + Ra) + tere —Ws 10°, (1) 


O = in (re — rm) + bc (Pe + re — Pm + Ri) + Wa Se 10°, (2) 


where i}, and i, are, respectively, the base and collector current increments. 


To determine the law of variation of the increment of current i,, we determined the quantity d@/dt. When 
the circuit is switched in, the flux changes from point a along arm a~2 of the magnetization curve (Fig. 2) and 
its magnitude can be expressed in terms of the current increment in the following way: 


D = — Dp, + AHS (Wk — W,) ir, (3) 





where 2 is the mean length of the core's magnetic lines and k = i/i;, is the gain. By differentiating equation (3) 
we find that 


dm 0,4 di 
See (Wak —W,) (4) 





After transformation, we solve the system of differential equations consisting of (1), (2) and (4), obtaining 


I (ry + ret Re) (re +o — hm + 21) — Pe (rg — Pm) 


7 eS 5 fel 5 
iy = tbo oxp 0,4 mpoS (Wek —W )(Wi(r, + 7. — hm + Ai) + Wer, t, () 





where i,, is the initial deviation of the base current. The expression obtained determined the law of variation of 
small deviations of current from its value in the equilibrium state, and defines the system's stability, the system 
being unstable when the exponent is positive. Since it is always the case that r, > I, > Ig, then for 


o (6) 
k> W, 
the size of the deviation grows without bound with increasing t and the system is unstable. For k < W,/W, the 


size of the deviation tends to be zero with increasing t, and the equilibrium state is asymptotically stable [5]. 


The condition for system instability is found after the value of k obtained from equations (1) and (2), is 
substituted in expression (6) 


WwW; 
_ W, "ot tet Rs) +r, — ar, > 1. (1) 





r +, + Ri +r, (1—a)]} 


It is clear from this inequality that to facilitate self-excitation, it is necessary to choose triodes with 
maximal a and to decrease impedances R, and Rg. To choose an admissible ratio, W,/W,, we solve the above 
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inequality for W2/W;. Ordinarily, W;/W, is not chosen smaller than 0.5. Since Te << 1,(1~a), then 


Ww 
resp [Rit re (1 —a)}. 


After simplification of expression (7) and some transformation, we get 





(Way ar. W, Ry, +r,.(i—a) 


quel) cme cree waar 8 
\W, rot+r.t+ AW, rot+ret+ Rs <0. (8) 


By solving the quadratic inequality thus obtained, we find the permissible values of W/W, 





ar, — V (ar? — 4(r, +7. + Re) [Aitr.(i—a)) WwW; c 
2(r, +7, + Re) WwW, 
ar. + V (ar, —4 (rp, +7. + Re) (Ait+r, (i—a)] (9) 
< 2(7, +7. + Re) 2 











For a triode with a = 0.9, Rc = 3° 10° ohm, R, = 200 ohm, r, = 10 ohm for Ry = 500 ohm and R, = 5000 ohm, we 
have 0.03 <W, W, < 53, For a triode with a = 0.98, r,, = 10° ohm, mr, = 400 ohm, r, = 10 ohm for Ry = 500 ohm 
and Re = 5000 ohm we have 0 < W,/W, < 172.5. 


Thus, as the transistor characteristics are improved, the range of admissible values of W,/W, increases. 
However, there is no practical necessity to broaden the range of the ratio obtained from the first transistor. 


Expression (7) is a necessary, but not sufficient, condition for the generation of oscillations by a Royer multi- 
vibrator, This is due to the fact that, although a small deviation of base current (5) will be increased if inequality 
(9) holds, it is not known whether or not, for given circuit parameters and source voltage E, the collector current 
reaches a value equal to the core’s coercive force, i.e., whether or not equation (4) will be valid on the entire 
segment a-2 of the magnetization curve (Fig. 2). For equation (4) to be satisfied on segment a~2 it is necessary 
that the following inequality hold: 


0,4n (i, W, —i, W,) 


(10) 


ws 





/\ 


After expressing i, and i}, in formula (10) in terms of the circuit parameters (Cf., Figs. 1 and 2), and after neglecting 
second-order terms, we obtain the second necessary condition for multivibrator operation 


E H Jar,r, 
> 04x W,[mr.a — m* (AR; + 7, — 7a) — ry — Te fa] , 





(11) 


where m = W,/Wg. 


Thus, the magnitude of the minimum voltage necessary for circuit operation is proportional to the coercive 
force of the core material, to the mean length of the core's lines (mean diameter) and to the emitter junction 
impedance, In modern transistors the latter quantity varies comparatively insignificantly (from 20 to 40 ohms). 
The number of turns of magnetizing winding Wg is given by the condition that the required frequencies be obtained 
and from the power requirements of the circuit (Cf., below). 


It is clear from an analysis of formula (11) that the quantity E has a minimum for 


ar, 
—— . 12 
2(fe—%e% + Ai) ( ) 





My 
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For the same coefficient m (for m = m,) the size of E decreases 
the increasing a, but this relation holds, basically, for values of 
m Close to my. Variations in the base and collector impedances 
has little effect on the size of E. 








Fig. 3. Equivalent circuit for one of the Il. Designing the Multivibrator Circuit 
amplifiers with positive feedback of the 
Royer multivibrator. 





The problem of designing a complete multivibrator circuit 
(Fig. 4) is the determination of the dependence of the output 
frequency and voltage and of the multivibrator power requirements 
on the input voltage, on the circuit parameters and on the characteristics of the core material. 


For simplifying the computations we shall assume that the multivibrator's output frequency is determined only by 
the cores magnetic polarity reversal, and does not depend on the transient response connected with the switching of 
the transistors, The latter assumption is valid since, for the frequencies considered here, the processes related to 
the switching of the transistors occur much more rapidly than the core’s magnetic polarity reversal. In the multi- 
vibrator circuit of Fig. 4, in contradistinction to the basic circuit (Fig. 1), there have been added regulating re- 
sistors Rs and Rg, made by tapping off a portion of the magnetizing winding with number of turns equal to Wy’, for 
the purpose of regulating the output frequency, and by the inclusion of an additional supply voltage Ey. With our 
assumptions taken into account, we consider the circuit of Fig. 5, equivalent to half the multivibrator. In this cir- 
cuit, the transistor is replaced by key K. The relationship between the circuit parameters and the time of flux 
growth in the core can be obtained from the following system of equations, set up for the circuit of Fig. 5 in 
accordance with Kirchoff's laws: 


E ss iR, +- ink, + W, ~ 108, 


(13) 

i,Ry + W, 10 = i,R,, (14) 

et, +h, (15) 

W, 10 + Ri =0, (16) 
E, = i,R, —W, 2 10", (17) 

H,+H,+H,+H,+H, =H.. (18) 


Here, Hf is the field strength induced by the eddy currents, H, is the field strength induced by current i 
flowing in winding W,; H, is the field strength induced by current i flowing in winding W,; H,' is the field strength 
induced by current i, flowing in winding W,’; Hs is the field strength induced by current i; flowing in winding Ws. 
Equation (18) holds as a consequence of the equality of the total field strength induced by currents i, fy, ig, i, 
and the eddy currents during the entire transient response and, on the other hand, the coercive force (with reversal 
magnetic polarity along portion 2-3 of the magnetization curve of Fig. 2). 








By solving the system of equations (13-18) we obtain, after some transformation, the differential equation 








0,4 -10-8 (Wo (1 + 2) (1+ @ + + 2cn + cn® + acn*) Ws. We 
I | cn R, [1 + a)? + a] >? + Rr d® + (19) 
_ {0,4nW (i + a) (1+ n+ an) 0,42W,F 
+ 0,8-10°% Byte dz = {OE ee a dt, 





where R = Rg + Rg, @ = Rg/Ry, n = We/Wy’, c = R/R, 7 is the specific conductivity of the core material, B,,, is the 
saturation inductance, the quantities z and dz are shown on Fig. 6 [6]. As the flux changes from — @,, + H, p95 to 
m + Hg poS, where S is the core’s cross-sectional area, the time changes from 0 to T/2 and z changes from 0 to 

d/2, where d is the thickness of the Permalloy strip. 
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Fig. 5. Functional schematic for investi- 
gating multivibrator core magnetization 
polarity reversal. 











To determine the multivibrator's output frequency, 
Fig. 4. Multivibrator with frequency tuning we integrate expression (19) 


by changing of the ratio Rg/Rg. 
0.4nW,(1+ a)(1+n+ an) . 0,4nW,2, 











po 28 niki {i +e+ecqy §—~—IRk, — Me 
=7= 5 - : 20 
1 1,6n-10-8@,, Sethe et mae ts. HS Wal 4m 10-°B vd? (20) 
l en*R, [(1 + a)* + ac} Re Ry ™ 


Formula (20) shows that, for a multivibrator constructed according to the circuit of Fig. 4 with the core 
material's coercive force and the eddy currents taken into account, and with the magnetic permeability equal 
to jig ON portions 1-2 and 3-4 of the magnetization 
curve and equal to » = o@ on portions 2-3 and 3-4, the 
pulse frequency is proportional to the input voltage and 








ae inversely proportional to the saturation inductance and (to 
a first approximation) to the core's cross-sectional area 
( &, = BS). 
1 
az 1. Effect of the Eddy Currents on the 
Se t Ry Frequency 





The effect of the Foucault currents on the frequency 





























a can be defined as the deviation of the output frequency f;, 
me rr determined without the eddy currents being taken into 
—d 7 ke account, from the frequency f determined with the latter 
L— taken into account 
Lo @ 4 5% — 21 100%. 


Fig. 6. For taking into account the 
effect of the eddy currents during 
magnetic polarity reversal, 


By substituting, instead of f, and f, the corresponding 
expressions obtained from formula (20), we obtain, after 
some transformation, 





{w2 (1 + a)(4 + a+ ¢ + 2cn + cn* + acn®) w? wil 





\ on, (+ a) + ac] +R, + Fy, 












It is clear from expression (21) that the deviation of the 





&% frequency from the true value when one ignores the eddy 
currents is proportional to the conductivity of the core 

_ ‘ material, the square of the strip's thickness and the mean 

4h length of the core's lines, and is inversely proportional to 
the core's cross-sectional area. This deviation does not 

Jb depend on the circuit's input voltage. 

J? | For cores constructed of alloy 65NP, with a strip 

aah thickness of 0.05 millimeters, the deviation is about 0.5%, 

74} 2. Effect of the Load Impedance 











To determine the effect on multivibrator operation 
of the load impedance, we estimate the deviation of the 
device's output frequency for a given load as compared to 
the frequency, f,, when the device is unloaded 


51% con pon 100%. 
1 














0s tk 6 wh f, 1. .kilohms If, instead of f; and f,, we substitute the corresponding 
expressions from formula (20), and after some 
Fig. 7. Deviation of the frequency for a given transformation, where 
reduced load impedance on the multivibrator's : 
output from the unloaded frequency. 1) a = 0, Rite = Ri(F) ; 


m=1;2)a=0,m=5; 3)a=o,m=1;4a= 
=o, m= 5, n= 5, c = 2, R = 1000 ohm, R = 500 








ohm. 
we obtain ye 100% ae 
: bs ° fata ett eee te m?*) - 
L re en*R, [(1 + a)? + acj + Rs} 


Thus, as the load impedance R, varies, the deviation of the multivibrator's output frequency from the un- 
loaded frequency varies along a hyperbola and depends neither on the core's geometry nor on the circuit's input 
voltage. Figure 7 shows 5;% as a function of Ry ye. 


It is clear from the curves of Fig. 7 that even in the worst case (a = a, m = 1) it is possible, for Ry. = 
= 15,000 ohm, to assume, with sufficient design accuracy, that the multivibrator operates on an infinite load (free 
running). With an output amplifier the multivibrator's load impedance can easily be brought to the value cited. 


3. Effect of the Slope of the Hysteresis Loop's Lateral Arms on the Dependence on 
Voltage of the Multivibrator's Output Frequency 


Formula (20) was derived from the condition that the lateral arms of the hysteresis loop were vertical. This 
is defined by equation (18). If the slope of the lateral arms of the hysteresis loop (Cf., Fig. 8) is taken into account 
then, in the system of equations consisting of (13-18), equation (18) must be replaced by the following relationship: 








® + Hy ps 


H, + Hy + H,+H2 + H,=—~;—. (23) 


If we neglect, for simplicity of the computations, the action of the eddy currents and the effect of the load 
resistance, which is possible if the conditions enumerated above hold, then we can solve the system of equations 
(13-17) and (23) to obtain 
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0.4%. 10-8 Piatt Re ast acn*) Wid 1 © (24) 
l en* Rf, ((1 + a)? + ac} +7} or +35 + 
__ 0,4nW, (1 + a) (1 + 2 + an) oe 0,40WiE, H 
nlR, {(1 + a)? + ac} r lRy - 


By solving equation (24) we find that 








©. —pSH, 1 aie 
© (t) = MpS —| MpS + Om + ty" Je MuQ . (25) 
where 
RE 0,4nW. (i + a) (1+ n+ an) 0,4nW,E, 
M= nlR, [((1 + a)? + ac} sa — H,, 
Sto. 04 US Ot 493 et ae w? 
veo om Ry (+ ap + ad] Ry) 


The time necessary to reverse the core’s magnetic polarity from point 2 to point 3 on the magnetization curve 
(cf, Fig. 8) is found from the formula 





MpsS —® i 

7 5 —On— type 

EPpy eRe ©, — v5 * sg 
MpS + ®,, + po se 


If we take into account that 1 >> 19, M > H, and B,, = 4,,,/S, we can write the expression for the output 
frequency in the form 


i 








j= i aa TG (27) 
2080 [I0(1+ aie)—I(1— are) 
Since By/My < 1 we obtain, after some transformations, 
i 
= 28 
(~ leat) (Fay es 
1+ 7 \ie +5 (mm FoF et Me 
where f,,, is the multivibrator frequency for y = a, and a ( =). is the general term of the 


series, 


The deviation of the multivibrator's output frequency when the magnetic permeability is » as compared 
to the frequency when yp = oo is defined by the expression 





bo = 100% -—— 


4 
: 100%. (29) 
oo 4 + 











1225 




















Fig. 8. Magnetization curve. 
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Fig. 10. Dependence on the applied voltage Fig. 9. Dependence on the magnitude of the 
of the error in determining 5, from the sim- input voltage of the deviation of the multi- 
plified formula, for H, = 0.2 oersted and vibrator's output frequency when the magnetic 
p = 50,000, permeability is » and the coercive force is H, 

from the frequency for an ideal rectangular 
Figure 9 shows the dependence of the deviation, hysteresis loop. 


53%, on the applied voltage for various p and He, 

computed from formula (29) for a multivibrator containing a core whose mean line length is J = 11.3 centimeter, 
B,,, = 10,000 gauss, for Wy = 500, Wg = 900, R = 1000 ohm, Ry, = o, Ry = 10,000 ohm, c = 2, a = 0,n = 4.5 andE = 0.4 
volt. With adequate design accuracy, and taking into account that B,, M,,< 1, we can determine the deviation 
from the simplified formula 


, B?, 
54% = — soap 100%- (30) 


Figure 10 shows the dependence of the error in determining 5,% from formula (30) instead of formula (29) 
for p= 2 


t, — bs 
—. 100% (31) 


for H, = 0.2 oersted and p = 50,000. 
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As is obvious from the graphs, the maximum value of n does not exceed 6.5%. Thus, the size of 5,, with 
sufficient accuracy, can be considered proportional to the square of the residual inductance and inversely pro- 
portional to the square of the magnetic permeability and the square of the applied voltage. 


In calculating the error in a telemetry system, it is desirable to know the error reduced to a nominal scale. 
This error can be calculated from the formula 


E 
bore %= By % (32) 


where Ey is the nominal input voltage to the circuit. 


It is clear from the analysis provided that for large residual inductances and for p < 10° there arises an essen- 
tial nonlinearity in the relationship between the output frequency and the input voltage. When alloy 65NP is used, 
the error from this nonlinearity does not exceed 0.1%, while for alloy 50NP it is less than 1%, This figures are 
valid only for well-annealed cores. 


4. Influence of the Impedance of the Magnetizing Circuit and of the Coefficient on 





a Circuit Operation 











For stable circuit operation, it is important that the operation be independent of the variations in the magneti- 
zing circuit impedance R,, a contribution to whose magnitude is also made by the emitter-collector impedance of 
the conducting triode. As shown by analysis of circuit operation, the dependence of output frequency on the 
impedance R, will be at its greatest when a = co. However, even in this case, as shown by an analysis of formula 
(20), the output frequency will not depend on R, if the following three inequalities hold: 


0,4 0,4nmE, A, 

IR, Es ih, + Ww,’ (33) 
A i m? 
R, > wR R, > R° 


It follows from an analysis of the inequalities in (33) that, for independence of circuit operation from 
variations in impedance Ry, one should decrease Ry as much as possible, increase impedance R, and choose a 
core material with the minimum coercive force. The coefficient m = W,/W, should be chosen with simultaneous 
consideration being given to the conditions of decreasing the minimum voltage necessary for generation of oscil- 
lations and the independence of circuit operation from variations in impedance Ry. These requirements contradict 
one another, and so a compromise solution must be reached. 


Figure 11 gives the experimentally obtained de- 
f, cycles a-« pendence of the frequency on the coefficient c for a 
184 f multivibrator with a = 0 and with a = oo (curves | and II). 
The same figure gives the dependencies of the fre- 
16+ 1 quency on the coefficient c as determined from formula 


(20), namely, curves III and IV. 
/4} 
a0 The discrepancies between the computed and the 


tab experimental data do not exceed 2 to 3%. 


It is clear from the graphs of Fig. 11 that a varia- 
tion of the coefficient a from 0 to oo, i.e., a translation 
of the potentiometer slider from one extreme position to 
the other, can change the multivibrator frequency by 
6 ‘ : . , +8 to 10%, This is particularly important for mass 

0 ! 2 J 4 de production, when the characteristics of all the transmit- 
ting telemetry devices must be identical and not depend 
on the variations in core material characteristics. 


0 


8 








Fig. 11. Experimental and computed curves of 
the multivibrator frequency as a function of the 
coefficient c for a = 0 and a = o. 
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5. Multivibrator Drain on the Voltage Source to be Measured 








The average current drawn by the circuit is 


bt 


TR 
ly =F \ idt, (34) 
) 
where i is the value of current defined by equations (13-18) 


En(i+a)—Wa(i+an +an)F 10-8 
i= ‘ (35) 


cmtbet Fes) 





By substituting (35) in (34) we obtain, after integration, 


E (1 +a) _ 4W,,f (1 +n -+ an) 10-* 
Ri(1+e+755) mR; (1 +0+75-) 








Tay (36) 


In formula (36) the frequency value substituted for f should be computed from formula (20) with account 
taken of the load impedance and the Foucault currents. Figure 12 gives the values of average current drawn by the 
multivibrator as computed from formula (36) (curve 1). The discrepancies of the computed data from the experi- 
mental (curve 2) are about 5 to 7%. 


It is clear from an analysis of formulae (36) and (20) that the 
current drawn by the circuit is proportional to the input voltage and 








. milliamperes the coercive force of the core material. The basic method of 
reducing the drainof the multivibrator on the transducer amounts to 

3+ decreasing the voltage E necessary for obtaining a given nominal 
frequency. However, the circuit must, with this, operate stably with 

2h a minimal frequency. Obviously, the minimum voltage at the cir- 

a cuit's input must satisfy condition (11). Starting from this condition, 

ls one must decrease the core's diameter and increase the number of 
turns, We. 

iy ote oe a ee me 


6. The Multivibrator's Output Voltage as a Function 
Fig. 12. Current drawn by the multi- of the Input Signal 
vibrator as a function of applied volt- 
age. 1) is the computed function and 
2) is the experimentally determined 
function. 








With an open output circuit, the voltage at the multivibrator's 
output is determined from the expression 


d® 1 
Eu, = Wa * 10% (37) 
The magnitude of d@/dt can be found from equation (19). Thus, the multivibrator's output voltage, in a first 
approximation, depends linearly on its input voltage. The discrepancies between the computed and the empirical 
characteristics do not exceed 15%. 


7. Speed of Response 


When a multivibrator is used in telemetering or computing devices, an important quantity is the time taken 
to establish its output frequencies after a change in the circuit's input voltage. Since the output frequency is deter- 
mined by the time taken to reverse the magnetic polarity of the multivibrator's transformer core, the transient 
response duration does not exceed half a period of the newly established frequency corresponding to the new voltage 
value, 
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SUMMARY 


1. The multivibrator's pulse frequency is strictly proportional to the circuit's input voltage only in the 
case when the magnetic permeability on the hysteresis loop's lateral arms (segments 2-3 and 4-1 on Fig. 2) equals 
oo. Otherwise, there is a nonlinearity error. This latter error is about 0.1% when well-annealed cores of alloy 
65NP are used, and about 1% when cores of alloy 50NP are used. 


2. Eddy currents decrease the multivibrator output frequency but do not distort the linear relationship be- 
tween frequency and input voltage. 


3. The variation in multivibrator frequency when the load resistance is lowered, as compared with the free- 
running frequency, is inversely proportional to the reduced load resistance and does not depend on the core geometry 
or on the circuit's input voltage. 


4. With the proper choice of parameters, the operation of the circuit depends very little on variations in 
the impedance of the magnetizing circuit, a contribution to which is made by the emitter-collector impedance 
of the conducting transistor. 


5. For tuning purposes, the multivibrator's output fréquency can be varied by 8 to 10% by moving a poten- 
tiometer slider (by changing the coefficient a from 0 to a). 


6. The current drawn by the multivibrator from the transducer is proportional to the input voltage and to the 
coercive force of the core material. 


7. The time to establish the output frequency, when the signal at the multivibrator's input changes, does not 
exceed half a period of the frequency corresponding to the new voltage value. 


8. The core material must possess minimal coercive force, a high value of maximum magnetic permeabi- 
lity and of saturation inductance and an adequate specific resistivity. Alloy 65NP fully satisfies these require- 
ments. 


Received February 25, 1959 
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FREQUENCY ELEMENTS OF REMOTE CONTROL 


I. N. Lisitskaya, V. N. Mikhailovskii, K. D. Nadtochii 
and V. N. Okhotskaya 


(L*vov) 


Certain types of frequency relays are considered. Data 
are provided on relays and frequency generators developed at 
the Institute for Machinery (mashinovedeniya) and Remote 
Control of the AN USSR. 


The majority of dispersed objects subject to remote control (irrigation systems, gas and oil fields, gas pipelines 
etc.) are generally interconnected by two-conductor telephone lines. When such objects are subjected to remote con- 
trol, it is reasonable to make use of frequencies above and below the audio range, since this makes it possible to have 
parallel operation of the telephonic communications link and the remote control apparatus. It is more efficient (eco- 
nomically) to operate in the sub-audio frequency range (f = 300 cycles), where the signal damping in the conductors 
is relatively small, which permits the use of apparatus over tens and even hundreds of kilometers while still guarantee- 
ing noise stability. One of the principal problems which arises with remote control (dispatcher control) of dispersed 
objects is the choice of method and means of selecting an object. 


By using the very promising frequency-combination method of selection with time relaying [1] for a small 
number of working frequencies (n = 10), well spaced in the chosen range (f = 300 cycles), and by using r= 3 non- 
repeating frequencies for each point, one can provide reliable selection of the necessary number of objects (N = 720) 
joined by parallel lines in which there is noise from such sources as magneto ringing of the telephone apparatus and 
lightening discharge. To implement this, it is necessary to have a highly stable characteristic matching of the trans- 
mitting generator and the frequency filter (relays). This matching must amount to this, that the possible bias, 6, in 
the transmitting generators frequency, fk: must be less than half the transmittance band of the filter, and the pass 
bands of filters with neighboring frequencies do not overlap. 


Naturally, one tries to arrange matters so that the operating range of signal voltages and the external temper- 
ature are as large as possible, the filters’ input impedance is as high as possible Rj, = 1 kilohm), the dimensions 
and the power dissipation are as low as possible, and the output power is sufficient for reliable operation of the 
simplest and cheapest small-dimension telephonic relay, for example, type RKM (P & 0.1 watt). 


The conditions just posed may be satisfied quite adequately by relays and master generators with electro- 
mechanical filtering elements, and even by crystal-triode relays with RC-filters. The most promising of the 
electromechanical filters are the reed filters. Slide-wire filters are difficult to build and have very low effi- 
ciencies, and tuning~fork filters are complicated to prepare and their quality is unjustifiably high for the con- 
ditions considered, which require only a strict stabilization of the master generators’ frequencies. 


Today, even the best of the well-known electromechanical frequency relays with reed filters possess certain 
design complexities, relatively long operation and release times, and sizeable dimensions. In this paper we shall 
attempt to investigate the possibilities of creating frequency elements with improved parameters. 


1. Contactless Frequency Elements Based on Reed Filters 








An electromechanical filter, as is well-known, works on the principle of electromechanical resonance of 
a steel lamina oscillating in the field of a constant magnet under the action of the variable magnetic flux of an 
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excitation winding. The resonant frequency of the electromechanical filter, not being the frequency of the 
lamina's free oscillations, depends not only on the elasticity Py and mass m of the lamina, but also on the magni- 
tude of the constant magnets’ inductance in the air gap, By, and in the gaps 4, and 4, (in the lower and upper 
portions of the magnetic system) and on the ampere-turns, Iw, of the excitation winding. 
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The expression for determining the resonant frequency of an electromechanical vibrator can be obtained 
from the equation [2.3] 





d?« ky dx B; Bolw | _ Bolws .. 
oe Es [a KK eM oe CTE (2) 


where x is the deviation of the lamina from its mean position, kg is thedamping coefficient, K and M are coef- 
ficients which depend on the parameters of the electromechanical vibrator , s is the cross-sectional area of the 
air gap and I is the current in the excitation winding. 


It is easy to find the expression for the electromechanical vibrator's resonant frequency from equation (1): 





4 B Bolw 


h = in “8~cTE eee 





M. (2) 


Expression (2) shows that increasing the induction B, and the current in the excitation winding, I, leads to a 
decrease of the resonant frequency , When B, and 1 are decreased, the resonant frequency of the electro- 
mechanical vibrator increases, approximating *o the value of the lamina‘s natural frequency of oscillation (fr -+ 
+ fy for By + 0 and 1 + 0). Change of the air gap also leads to a variation in f,. When the gap is increased the 
frequency f; increases, and tends to fy, but the dependence of the frequency on the magnitudes of By, I and 6 is 
decreased. Figure 1 gives the analytically obtained functions f;/fy = ¢@2/5n) for Uj, = const, and f,/fyy = 
= ¢_(Uj,/Un) for 5, = const (where 5 n is the size of the magnetic gap for which the efficiency n is maximal, 
and U;,,/Uy is the ratio of the input signal's voltage to the voltage taken as nominal), which coincide well with 
the experimentally determined characteristics, shown by the dashed lines. 


For carrying out experiments we prepared specimens of the electromagnetic filters for frequencies of 60- 
300 cycles, the laminas having a thickness of 0.35 mm, a width of 88 mm and lengths of 30 to 40 mm, With 
this, the size of the gap was 5:n © 0.5 mm, The filter dimensions were 25 X 35 X 60 mm. 


The characteristics given in Fig. 1 show the decrease in the dependency of the electromechanical vibrator's 
frequency on the instability of the gaps and the input signal levels as the gaps 6, and 5, are increased. Curve 1 
was obtained for 5; = 0 and curve 2 for 5; = 5yN. To increase the stability of the electromechanical filter's 
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resonant frequency (which is particularly essential for a 
frequency generator) it is necessary to guarantee that it 






































operates on the section of the curve for which 5, / dyy > 1. 
w, > a As a result of the investigation of several generator 
Deahe circuits, we chose the circuit shown in Fig. 2, consisting 
of a transistor amplifier with a feedback path through a 
+ uy selection element — an electromechanical vibrator. The 
+ out feedback mode is established by the choice of Rg, Rg and C. 
-< The proper choice of resistors Ry and Ry allows high temper- 
Wy > eid ature stability to be achieved, whereby a variation of the 
LL | ambient temperature from +20 to +60°C causes a variation 
a sic of voltage U,,,, which does not exceed 45%, and the 











variation of the generated frequency is less than 0.01% per 
Fig. 2. degree C. If necessary, the temperature stability of the 
frequency can be increased by replacing the ordinary spring 
steel,of which the vibrating laminas were made by steel with a low temperature coefficient. 


The influence on the zenerator's output voltage of the method of mounting the filter on the chassis was 
avoided by using a rubber shock absorber with a corresponding distribution of the filter's mass. 
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Fig. 3. 


The best of the contactless electromechanical frequency relay circuits investigated turned out to be the 
ones shown in Fig. 3,a and b*. In the circuit of Fig. 3,b a cathode follower was introduced in order to match the 
impedances of the electromechanical filter's winding and the output amplifier. 


The input power of the frequency relay necessary for reliable operation of relay RKM, connected at the out- 
put, (Rwi = 700 ohm, P,,,, = 0.15 to 0.2 watt) was 5 milliwatts. The relay's input impedance was Ry, % 3000 ohms, 
The nominal input voltage was Uys = 5 to 7 volts, the relative working frequency band, 2 AL (100%)/ Ff, ay’ Was 
8 to 12% when U, = 1.3 Uy , and 3 to 5% when Us = 0.7 Uys. 


When the ambient temperature varied from 0 to 40°C there was an insignificant (< 0.5%) narrowing of the 
working frequency band. There was no relay mistuning with variations in temperature. 


The working range of the frequency relay's input signal voltage and, consequently, the reliability of the 
selection system operation under variation in communication channe] parameters can be broadened by a factor 
of from 5 to 10 by connecting a signal amplitude limiter (for example, one composed of four crystal diodes) at 
the relay's input. In the frequency range AF = 45 to 315 cycles, we succeeded relatively easily in positioning 
nine or ten frequency channels. The physical construction of the filters was essentially simplified by the replace- 
ment of the ordinarily used magnetic circuit apparatus by two plates of ordinary soft steel of 1.5 mm thickness. 
This permitted a sufficiently broad pass band and a small damping time to be obtained. 


* This circuit was investigated by Yu. Yu. Sikachevskii. 





2, Electromagnetic Contact Frequency Relay 





The functional schematic of an electro-mechanical contact 
frequency relay is shown in Fig. 4. The contacts of the electro- 
mechanical filter relay are connected to a single-stage semicon- 
ductor amplifier, whose load is the winding of a small-size tele- 
phone relay (RKM). Capacitor Cy is used for smoothing the pulsa- 
ting current in the base circuit, and capacitor C, is used for smooth- 
ing in the collector circuit and for increasing the relay's operate 
time (to increase the relay's stability to noise pulses). 
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In normal operation, the voltage on the frequency relay's 
contacts is about Uy, = 0.8 volt and the current flowing through 
the relay contacts is I, = 0,5 milliampere. With such voltages 
and currents, contact erosion due to electrical phenomena does 
not occur [4, 5]. The working frequency bands cited earlier for the contactless relay are obtainéd with input 
signal levels of Us =(1 +0.5)Uyy ,; the nominal value of signal voltage can be taken as equal to Uy * 5 volts. 


Silver was used as the relay contact materials. The relay's operate band varied relatively little (< 0.5%) 
in the temperature range from 0 to 40°C. 


With the nominal signal voltage Uys, the operate time of the frequency filter together with the RKM relay 
connected at the output lay within the limits of 0.1 to 0.3 second. Thanks to the retardation of operation of the exe- 
cutive relay, the frequency relay had low sensitivity to noise pulses and shocks. When dc voltage pulses, U, = 10Uy . 
were applied to the frequency filter's winding, the frequency filter did not operate. The power for reliable relay 
operation was approximately 5 milliwatts, the input impedance was R;, = 3 kilohms, the dimensions were 30 x 30 x 
x 60 mm. The bracing of the relay was standard, analogous to the bracing of an RKM telephone relay. 


A description of an analogous relay with somewhat different construction is provided in [6]. 


3. Crystal Triode Frequency Relays 











The circuit for a frequency relay realized with plane crystal triodes (Fig. 5) consists of a resonant RC-ampli- 
fier with a twin-T bridge in the negative feedback path. As is well known, the operation of a twin-T bridge in 
amplifier circuits tends to provide a no-load mode and a given voltage since, in this case, the maximum Q factor 
of the circuit is reached. In order to make the operation of the twin-T bridge approach the no-load mode, the 
bridge output in the circuit shown was loaded by a triode connected in the circuit with a grounded collector. The 
output triode, KT4, to whose collector is connected the winding of an RKM-type relay, is a power amplifier. Its 
input impedance is low and therefore it is connected to the output of amplifying triode KT, via cathode follower 
KT; whose input impedance is slightly shunted by the output of triode KT. If a sensitive output relay is employed 
(polarized, magnetic, etc.), the necessity of triodes KT; and KT, diminishes, The winding of such an output relay 
can be connected in the collector circuit of KT, via a dividing capacitor and a succeeding rectifying device. 






































Fig. 5. 


Analytical investigations of the circuit showed that its gain K and its Q-factor Q at the quasi-resonant 
frequency depend in a complicated way on the circuit parameters and the twin-T bridge elements: 
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Te, Th» Tj, and tp, are parameters of the crystal triodes, Ww») is the quasi-resonant frequency and gp is the input con- 
ductance of triode KTs. 


For the chosen frequency range (60 to 300 cycles) and m = n, the optimal value of the Q factor is obtained 
from m = 0.5 to 0.7, y = 0.7 to 1, gs = (1 to 2) X 10 mho,. The remaining parameters of the twin-T bridge are 
chosen in accordance with the condition that it be tuned to the given quasi-resonant frequency [7]. Experimental 
data on the relay prepared at the IMA showed the practical independence of the width of the operate band from the 
temperature in the range 0 to +35°C, adequate stability of the resonant frequency (drift of less than 0.5%) with 
variations of 250% of the supply voltage and variations of the magnitude of the input signal, U, $(1 +0.5)Uy ¢. 
The width of the operate band increased with increased supply voltage and magnitude of input signal. 


The relay's input impedance was Rj, = 40 to 50 kilohms, the operate power was Pop © 2 microwatts. For 
the nominal value of the input signal voltage, the width of the operate band was from 3 to 5% of the resonant 
frequency. The relay dimensions were 25 x 35 x 70 mm, 
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USING TRANSISTORS TO’INCREASE THE EFFICIENCY OF 
REVERSIBLE DC MAGNETIC AMPLIFIERS 


M. A. Rozenblat and G. V. Subbotina 


(Moscow) 


A new circuit for a push-pull magnetic amplifier with dc load is 
described, the high efficiency of this circuit being obtained by the use of 
transistors as switches which permits the harmful effects of the halves 
of the circuit on each other to be practically eliminated. 


The paper provides results of the investigation of a magneto-semi- 
conductor amplifier built in accordance with the high-efficiency cir- 
cuit presented here, and the data are also given for a high-stability 
measuring amplifier constructed on the same basis. 


The well-known push-pull magnetic amplifiers with de loads are distinguished by their common disadvantage — 
low efficiencies, not exceeding 50% of that in the most effective circuits, even when there is no loss in the windings, 
rectifiers or cores. The reason for this phenomenon comes down to the shunting interaction of the two halves of the 
push-pull circuits which leads to an undesirable cross distribution of currents. To overcome this limitation one con- 
nects a ballast resistor at the amplifier output, and 50 to 83% of the output power is lost in this resistor [1]. 


New possibilities for improving the characteristics of magnetic amplifiers, push-pull amplifiers in particular, 
are revealed by the joint usage of them and transistors [2-5]. ? 


In the present paper we give the results of the development and investigation of a push-pull magnetic 
amplifier based on a new circuit.® In this circuit, the transistors act as switches which, in dependence on the 
sign of the signal, permit the load to be switched to the output of that one of the two single-ended amplifiers 
whose output current is increasing. With this we succeeded in almost completely using the power of the single- 
ended amplifiers when they operate in the push-pull circuit with a single dc load. 


Increasing the efficiency of a push-pull amplifier is accompanied both by an increase in the power yield 
and by an increase in the power yield and by an increase in the extent of the linear portion of its characteristic 
and in its power gain. With this, the amplifiers inertia is not changed, which gives a significant increase in its 
dynamic energy factor. 


Below, we provide the theory of operation of, and the results of experimental investigation of, a magneto- 
semiconductor amplifier, as well as giving the parameters of a high-stability measuring amplifier based on the 
circuit considered here. 


Circuit and Theory of Operation of the Amplifier 





The circuit of the magneto-semiconductor amplifier (Fig. 1) is an ordinary differential push-pull magnetic 
amplifier with internal feedback in which the rectified current from the single-ended arms is applied to the load 
via emitter-collector junction transistors T1 and T2. 


*M. A. Rozenblat. Magneto-semiconductor amplifier [in Russian]. Patent application (Avtorskaya zayavka) 
623120/26, 
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The special feature of the circuit is that there is applied to the triodes, which are connected by a common 
emitter, a controlling voltage which is proportional to the output voltage of the corresponding single-ended 
magnetic amplifier. The magnitude of the controlling current is established by a divider made up of resistors Ry 
and Rp. 


With the application of a control current of a definite sign to the magnetic amplifier's input, a large cur- 
rent flows from one arm to the load, generating a significant voltage between the emitter and base of the cor- 
responding triode and opening it, while the base current of the second triode is sharply dropped and the degree of 
opening of this triode is decreased. 
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Fig. 2. Graphs for explaining the theory of operation of the push-pull 
magnetic amplifier with triodes. 


Thus, the second triode will present a large impedance to the first arm's current which both limits the 
harmful current shunting and increases the effectiveness of the entire circuit's operation. Then the signal current 
changes, the roles of the triodes are changed and a current flows in the load in the opposite direction with the 
same result as far as circuit effectiveness is concerned. 
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If, with minimum current in an arm, the corresponding triode is completely cut off, then the maximum cur- 
rent in the load of the push-pull amplifier will equal the maximum current of its single-ended arm. 


A graphic presentation of the operation of the amplifier can be obtained if we consider the graphs of the 
output voltage of a magnetic amplifier whose cores have rectangular hysteresis loops. We assume, for simplicity, 
that in the absence of a signal the saturation angles of both single-ended amplifiers isO , = 6, = 1/2, thereby 
obtaining the curves of the output voltages Uy and U, which are given in Fig. 2,a. In this case, the magnetic 
amplifier'’s output voltage U; equals zero, independently of the presence of the triodes. 


When a signal is applied to the amplifier's input, the saturation angle 6, of the first single-ended amplifier 
decreases while that of the second, 6, increases by the same amount (Fig. 2,b). In this case, during the interval 
0= wt = 4; there are, at the outputs of both amplifiers, small, practicaliy identical, free-running voltages which, 
with the proper choice of triode mode of operation, are insufficient to open the triodes. Initial triode collector 
current I, flows through the emitter-collector junctions, and the resulting load voltage (current) equals zero. In 
the interval 9; = wt = 6, a voltage is applied to the input of the first triode (in the emitter-base circuit) which 
suffices to open the triode, while the second triode remains closed as before. Therefore, current 1, at the first 
amplifier's output is, for all practical purposes, not shunted through the second amplifier's rectifier, but falls 
completely on the load. 


In the third interval, 6, = wt= @, after the core of the second amplifier has been saturated, both triodes 
conduct. Since the output voltages (currents) have identical instantaneous values, the load voltage in this interval 
will be zero. Therefore, the output voltage U; (or current) of the amplifier has the form shown in the lower graphs 
of Fig. 2,b. The forms of the curves for voltages Uy, U; and U; with maximum signal at the amplifier's input are 
shown in Fig. 2,c. In this case, the second triode is cut off during all the periods and the load current is virtually 
equal to the current at the output of one single-ended amplifier. 


In an actual amplifier the divider parameters Ry and Rj, must be so 
chosen than the triodes do not limit the magnitude of the current at the out- 





















































1, milliampere ‘ puts of the single-ended amplifiers in the interval of saturated cores, 
108 6=< wt = . This condition will be met if, for @ = wt, the collector 
90 ] current is 
i | ic = Ten + q—z ie >i, 
10 
60 where i is the instantaneous value of the current at the output of a single- 
» ended magnetic amplifier, I,, is the initial collector current, ip is the base 
| current corresponding to amplifier output current i and established by means 
40 of resistors Ry and R},, « is the current gain of the triodes with a common 
base. It should however, be borne in mind that, in a single-ended circuit, 
a PIER Y the load current is i, = i~i,. Therefore the choice of the triodes’ mode of 
20 operation, as determined by R, and Rp, must also guarantee, for maximum 
| efficiency, that the value of base current ip is much less than the current i. 
Oe ] In practice; it suffices if i,, < 0.1i. Implementing these two conditions is 
ae: easily done, since the current gain of triodes with common emitters, [a/(1 —a)], 
. ad ye @ 4S, volt is ordinarily not less than 15 to 20 and, in the best types, is considerably 
greater. 


Fig. 3. The function 1, = f(Up) 
for transistor 313B. The presence of a magnetizing current (voltage) at the output of the 
second single-ended amplifier in the interval 6; swt = 9, can somewhat 
open the corresponding triode, which leads to a certain lowering of the load 
current of the push-pull amplifier. Ordinarily this effect is negligibly small due to the nonlinear character of 
the function I, = f(Up) in the region of small values of Uj, (Fig. 3). 


It is clear from the theory of operation we have been considering that the choice of parameters and design 
of the push-pull magnetic amplifier reduces, essentially, to the choice of parameters and design of the correspond- 
ing single-ended magnetic amplifier. If the triode parameters are properly chosen they have practically no effect 











on the operation of the single-ended amplifier. It is not difficult to choose the necessary type of triode once the 
amplifier’s output parameters are given, The triode’s maximum current equals the maximum current of the 
single-ended amplifier. The maximum voltage on the triode occurs when the output current of the correspond - 
ing single-ended magnetic amplifier is a minimum, and equals its maximum output voltage. Therefore, the 
triodes are chosen in accordance with the conditions that Ic agm = Itmax 294 Ucadin = U1 max: 


We note that the suggested method of connecting triodes is completely applicable to other push-pull 
amplifier circuits with-de outputs, including circuits with external feedback paths, used for increasing zero stabi- 
lity [1], for center point rectification circuits, etc. 


EXPERIMENTAL DATA 


Experimental investigation of a number of amplifiers based on the circuit of Fig. 1 showed that the push-pull 
amplifier’s maximum load current is less than the current at the output of the corresponding single-ended amplifier 
by not more than 5-to 10%, which was due to the partial shunting of current i, through the triode's emitter-base 
circuit and to the effect, cited above, of the free-running current of the second single-ended amplifier (Fig. 4, b 
and c). We note that a further approximation of the push-pull amplifier's load current to the load current of the 
single-ended amplifier can be attained by using a nonlinear impedance for R}, (for example, a thyrite type). In 
this case, a decrease of voltage at the single-ended amplifier's output increases the magnitude of impedance R;, 
which leads to an additional limitation of the current through the given triode. The use of nonlinear impedances 
and of other well-known means for decreasing magnetic amplifiers’ free-running current are advantageous in cases 
when the single-ended amplifiers have a small range of load current variation. 


Figure 4 shows the characteristics of one and the same push- 
pull amplifier: a is for an ordinary circuit with ballast resistors 
l,, ma and b is for the suggested circuit, with the ballast resistors re- 
placed by triodes. In the second case, with identical load impe- 
dance, the maximum load current is increased by a factor of 
approximately 2.4 and the power is increased almost sixfold. 


The amplifier has the following parameters: cores of 
alloy 65NP with thickness of 0.15 mm, lateral cross-sectional 
area of 1 cm’ and mean diameter of 4 cm, a supply voltage 
of 9 volts at 50 cycles, w., = 300 turns, w,, = 3500 turns, Rip = 
= 600 ohms, R; = 20 ohms, triode type 313B. The power gain of 
the circuit of Fig. 1 was 100,000. 


It is clear, from the amplifier's theory of operation, that 
the triodes operate in a noncritical mode and do not, therefore, 
have much effect on the amplifier characteristics or the 
stability of the null. This was verified experimentally. Replace- 
ment of the triodes without any special selection having been 
made had virtually no effect on the amplifier characteristics. 
Ha The choice of resistors Ry and R), is also not critical. 


The increased dynamic energy factor of the amplifier due 
to the decreased power loss at its output, and the virtual absence 
of triode effect on the amplifier's null stability allowed it to be 
used for implementing a dc measuring amplifier. For this pur- 
pose, a negative electric feedback path was connected from the 
amplifiers output to its input. With this, the positive feedback 
path was increased by several turns for the purpose of increasing the open-loop gain. By doing this we were able, 
with a single-ended amplifier whose parameters were given above, to obtain, with no special selection of circuit 
elements and with a gain of 2000, a constancy of gain to within 0.2% when the supply voltage varied by 10 to 15% 
and the load impedance varied by a factor of three. With this, the amplifier's time constant equalled three periods 
of supply voltage. The null drift of the amplifier over an 8-hour period of continuous operation corresponded to a 
signal power of 1.5 x 10°" watts. With a stabilized voltage supply, the accuracy of the amplifier's operation 
increases. 


Fig. 4. Magnetic amplifier characteristics: 
a is for a push-pull circuit with ballast resis- 
tors, b is for a push-pull circuit with triodes 
and c is for a single-ended amplifier. 
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SUMMARY 


The method suggested for increasing the efficiency and the energy factor of push-pull magnetic amplifiers 
with dc outputs is based on the use of transistors operating as controlling switches. With this, the triode (transis- 
tors) parameters are not critical and, if properly chosen, have little effect on the single-ended amplifier character- 
istics, 


The suggested method of constructing push-pull dc magnetic amplifiers can also be used both for low-power 
high-sensitivity amplifiers and for power output stages. The maximum power of such amplifiers is determined by 
the parameters of the triodes. 
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INVESTIGATION OF THE OPERATION OF A MAGNETIC AMPLIFIER 
WITH SELF-SATURATION ON A THREE-PHASE LOAD 


A. L. Pisarev 


(Moscow) 


The steady-state operation is considered for a magnetic amplifier with 
self-saturation operating on an active three-phase load for amplifiers with free 
and with suppressed ac components in the control windings. 


Methods are given for constructing the control characteristics, and the 
question of the voltage on the rectifiers is considered for both cases. 


The use of magnetic amplifiers with self-saturation for controlling three-phase loads has a number of 
advantages as compared with the use of ordinary saturated coils. 


The basic advantage is the significant increase in gain which, in a number of cases, allows one to do without 
intermediate amplifiers, and to reduce significantly the size and weight of the entiredevice. The presence of 
rectifiers in the load circuit complicates the device but, as will be shown below, the quantity of these can be 
small, 


























Fig. 1. 


The circuit of a magnetic amplifier with self-saturation and a three-phase load is shown in Fig. 1. 


For the analysis of the circuit we make the following assumptions. 








ut 


. The circuit is acted upon by a symmetric three-phase system of sinusoidal voltages. 





2. The amplifier has an active load. 







3. The rectifiers have a negligibly small forward impedance and an infinite back impedance. 


4, The dynamic hysteresis loop of the core material is approximated by the parallelogram with the para- 
meters shown in Fig. 2. 































5. The active impedances of the amplifier's load windings 





Vie are assumed to be negligibly small. 
oT a E 6. The supply voltage (phase) is chosen in accordance with 
A 8 the relationship 






ae "Fi Wi, w-1078, (1) 











The basis of expression (1) are given in [1]. 











7. The magnetizing current is negligibly small in comparison 
pH with the current in the saturated portion of the period. 





It is known, from amplistat theory [2], that their properties 
depend essentially on whether or not ac components can be found in 
the control winding circuit. In three-phase circuits, just as in single- 
phase ones, two extreme modes of amplifier operation can be distin- 
guished: operation with free ac components in the contro! winding 
and operation with suppressed ac components. A magnetic amplifier 
(Fig. 1) contains six reactors, two per phase. Each reactor is denoted by a double subscript. For example, the 
subscript aa denotes reactor a phase a, etc. The reactor control windings may be either individual for each 
reactor or common to each pair of reactors (as shown on Fig. 1). 








Fig. 2. 


If the control windings of all phases are connected in parallel to a source of control voltage with low internal 
impedance, then the mode of amplifier operation with free ac components in the control windings will occur. 


If the control windings are connected in series to a control voltage source with high internal impedance 
then the mode of amplifier operation with suppressed ac components in the control windings will occur. It should 
be mentioned that, when the phase control windings are connected in series, even with a very low impedance in 
the control circuit, the magnetic amplifier's mode of operation will approximate to the mode with suppressed ac 
components in the control windings, since the even current harmonics, not multiples of three (second, fourth, 
etc.), cannot flow in the control circuit(their sum is zero). 


We shall consider further these two modes of magnetic amplifier operation. 


1. Amplifier with Free ac Components in the Control Windings 








Similarly to single-phase systems, the following conditions will be valid for each pair of reactors. 


1. Due to the presence of a very low-impedance contro] winding-loop shunting the pair of reactors, the 
variations in the reactors’ magnetic flux will occur simultaneously and with identical velocities. For the reactors 
of phase a, for example, the following relationship is valid: 
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2. When one of the phase reactors is saturated, the magnetic flux of the other reactor maintains an invari- 
able value 





Paa = V,. (3) 





hap = const for 











In [1] there was considered the operation on a three-phase load of a magnetic amplifier with series-con- 
nected load windings. It was shown that, for an amplifier with free ac components in the control winding, three 
modes of operation of the load circuit could occur, depending on the size of the saturation angle ®,, these modes 
being differentiated by the form of the load current curve and by the characteristic equation of the circuit's 
electrical equilibrium. The boundaries of these modes are defined by the following values of the saturation angle: 


for the first mode 





TSS ei (4) 
for the second mode 

7 << FS; (5) 
for the third mode 

0<%&< >. (6) 


We now show that when there are free ac components in the control windings, the operation of the load cir- 
cuit of the magnetic amplifier with self-saturation which we have been considering is analogous to that of the 
amplifier without self-saturation which was considered in [1]. 


In fact, the two conditions cited above on the magnetic flux variations hold in both cases. The closed loop 
of the control windings, in both cases, guarantees that equality (2) holds. The difference amounts to this, that in 
the interval of magnetic flux variation the network voltage is divided equally between the two reactors of a given 
phase in the magnetic amplifier with series-connected load windings while, in the analogous interval for the given 
case, all the net voltage is applied to each of the reactors. 


Moreover, when one of the phase reactors is saturated in the magnetic amplifier with series-connected load 
windings the load current flows in the load windings of both reactors, while in the analogous interval of the given 
case, the load current flows only in the winding of the saturated reactor. With this, the load circuit of the second 
reactor is cut off by a rectifier. In the given case, the reactors operate alternately. 


These special features however, do not disturb the analogy of the reactors’ operation as electrical circuit 
elements, since the same three modes of load circuit operation occur in the case under consideration as in the 
case of a circuit with series-connected load windings. The form of the curve for current in the load is identical 
for both circuits. Figure 3,a shows the curve of the load current for phase a for the case corresponding to the 
mode described by (5). This same curve can be seen on the oscillogram of Fig. 4. 


Based on the analogies of the physical processes in the reactors, the equations for the electrical equilibrium 
of the given circuit, set up in segments defined by the magnetic states of the reactors, are identical with the 
equations which define the operation of a magnetic amplifier with serles-connected load windings in the analogous 
modes [1]. The only difference is that, in the equations defining the increments of the reactors’ magnetic fluxes, 
the denominators will contain W, instead of 2W,. The validity of this substitution is obvious and requires no 
explanation. 


Figure 3,b shows the curves of the magnetic fluxes of reactors aa and ag. The expressions for these mag- 
netic fluxes are defined as follows. In the interval Osa0 ~ <6 < Q. ox (Fig. 3) reactors bg and ca are saturated 


and the magnetic flux of reactors ba and cg are invariant. The magnetic flux of reactor aa increases and that of 
reactor ag decreases. 


The following equations are valid: 


tse — (ic — in) R = 0, (1) 
Uca — Caa +4- i-dt = 0, 


le —— lp. 
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Fig. 3. Relationships for the magnetic amplifier with free ac 
components in the control windings. 


The last two equations are correct if assumption 7 is 
taken into account. 


'L We find from equations (7), (8) and (9) that 





u 
€Caa = Uca + + . (10) 


For 6 = O00 — 5 Paa = D,, Pap = D,, and 


Pe Me / \ / therefore the expressions for the instantaneous value of 
magnetic flux will be 


Fig. 4. Load and control current oscillograms 
for a magnetic amplifier with free ac compo- 
nents in the control windings. 
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Figure 3,c represents the magnetizing component of the load current of phase a (in a scale different from 
that of the current in the saturated portion of the period). 


In the case considered here, the processes which occur in each reactor pair in the interval of magnetic flux 
variation are qualitatively similar to the processes in a single-phase amplistat [2]. If the slope of the lateral 
sides of the material's magnetization characteristic is not taken into account, the magnetizing component of the 
load current will be represented by a rectangle with ordinate 2Hci o/W, and an abscissa determined by the duration 
of the interval of magnetic flux variation. With this, the control current is defined by the expression 
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However, tests showed that not taking the slope of the magnetization characteristic's lateral sides into 
account leads to a significant discrepancy between calculated and experimentally determined amplifier control 
characteristics. 


We consider the half-period of supply voltage in which the magnetic flux Poa increases while magnetic 
flux Qap decreases. The magnetizing current flowing through the load winding of reactor aa (Fig. 3,c) includes 
the current for magnetizing reactor ag (its ordinate ishatchedin the diagram) and the current for magnetizing reac- 
tor ac itself. With this, the segments 1], J,’, ly and 1j on the current curve correspond (in a changed scale) to the 
segments J; , Ig, lg and lg on the magnetization characteristic. The instantaneous value of the control current in 
this case, as for a single-phase amplifier [3], obeys the equality 


: i. 
lea = — luap Ww ° (13) 


In the following half-period the reactors change places. The control current, icg, of phase a is shown in 
Fig. 3,d. 


ae For illustration, Fig. 4 gives the oscillograms of 


the phase voltage Up, the load phase current i, and the 
phase control current ic. Since, in the three modes 
defined by (4), (5) and (6), in the interval of magnetic 
flux variation in the reactors, the physical conditions 
are different, then the expressions for the control current 
in these modes will also be different. 


The average value of the control current will be 
defined by the following expressions (their derivation is 
given in Appendix I): 


for the first mode (4) 


¥s = Sr {(%— *)_ 





(14) 
— v(— 1.05 + 1.38 sin 0,— 0.4 cos 6,) | ; 
for the second mode (5) 
Ie, 
ma 
oe l 
Fig. 5. Computed and experimentally Ip = lF ~ (5 + 0.6 cos 9,— 0.7 sin ®,) | . (159 
y - 


determined control characteristics of a 
magnetic amplifier with free ac compo- 
nents in the control windings. 
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for the third mode (6) 


A ae m 
le = yt 1% — ¥ (1.73 sin 9, — 0.730,)}. (16) 





For the enumerated modes, the mean value of the load current is expressed as follows (the derivation of the 
formulae is given in Appendix II); 


for the mode of (4) 








t= ae [t toon (0. + &) an 

for the mode of (5) 
h = Voom sin (9, + 3): (18) 

for the mode of (6) 
I = m4 }- cos 95). (19) 


Thus, in expressions (14) and (17), (15) and (18), and (16) and (19), Ic and I, are related to each other via 
the parameter 9,. 


Computation of the magnetic amplifier's control characteristics is carried out by determining the values of 
Ic and I; for a number of given magnitudes of 8. 


When the magnetizing component of the load current must be taken into account, this can be done approx- 
imately by assuming that 1, =1,. For 6, = 59/6, it is easily seen than Iy = 21.W./W, (since with this A@ = 
2, and, instead of partial hysteresis cycles, there occurs a symmetric cycle and, consequently, the equation {,.= 


= Zines is valid). 


By drawing the quantity I,, (segment AB on Fig. 50) and then joining point B with point O, we obtain the line 
BO whose ordinates are approximately defined by the magnetizing component of the load current I. 


Figure 5 gives the working portions of the computed control characteristic, 1, for a magnetic amplifier oper- 
ating on a three-phase active load with phase control windings in parallel. Below is given the data for the ampli- 
fier for which the characteristic of Fig. 5 resulted. 


Core data: the cores are toroidal; Zl, = 11.55 cm; q, = 0.685 cm’; core material was 65NP; dynamic coer- 
cive force was H, = 0.26 ampere- turns per centimeter; v = 0.27. Circuit data: R = 285 ohms; Up = 29 volts; 
W,, = 1500; W,, = 200. The rectifiers were type DG-Ts24. 


Figure 5 also gives the experimentally determined control characteristics (curve 2). 


2. Amplifier with Suppressed ac Components in the Control Windings 








As was already stated, this mode of magnetic amplifier operation can occur when the control windings of 
the individual phases are connected in parallel and supplied from a source with a significant inte.nal impedance. 


We consider the processes occurring in the reactors of one of the phases (for example, phase a). With the 
assumptions that were made regarding the smallness of the load windings’ active impedance and the rectifiers’ 
back impedance, we can easily show that, in this case, conditions (2) and (3) also remain valid. 
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A pair of reactors of one of the phases can be likened to a single-phase magnetic amplifier with self- 
saturation and suppressed ac components in the control winding, the properties of which are very well known [2, 
4, 5]. If, at some moment of time, reactor ag ceases to be saturated then, starting with this moment, the reactor 
begins to act like an inductance. The current in its winding retains its previous direction. For this interval, the 
following equation is valid: 





d 
Uag — WLw ne —i,R =0. (20) 
Here, Ua is the voltage between points a and 0 (Fig. 1). 
d 
Part of voltage ug), equal to Uaa = Wyo et , is also applied to reactor aa, forcing magnetic flux 


aa ‘o increase and rectifier aa to begin switching. For the load windings loop, if the active impedance of 
this loop is neglected, the following expression will be valid 





dag qIP iq 
where 

CL — dP 

a: 


Equality (2) holds during the entire interval of reactor flux variation and, therefore, the intervals of magnetic 
flux variation of both reactors coincide in time. 


When reactor aq is saturated, the load winding of reactor a g is short-circuited and, consequently, @ ap = 
= const. It thus turns out that, with the assumptions made, conditions (2) and (3) are met in magnetic amplifiers 
with either free or suppressed ac components in the control windings. The difference amounts to this, that in the 
first case the meeting of the condition on the smallness of the load windings’ active impedance and the rectifiers’ 
back impedance is optional, since the closed loop of control windings of each pair of reactors necessarily guarantees 
that conditions (2) and (3) will hold. In the second case, it is mandatory that this condition be met. 


In practice, this condition is met in relatively powerful magnetic amplifiers, where R >> r and r is the 
active impedance of the load winding. The case of finite magnitudes of r and r;, where r, is the rectifier's for- 
ward impedance, will be considered qualitatively below. It is obvious that when conditions (2) and (3) are met, 
the processes occurring in the magnetic amplifier's load circuit in the given case will be completely identical 
with those in the amplifier with free ac components in the control windings. 


The relationships between the average value of the load current and the saturation angle, (17), (18) and 
(19), remain valid. 


We turn now to the determination of the relationship between the control field strength Hy and the magni- 

tude of the reactor's magnetic flux. It is known from the theory of single-core magnetic amplifiers [5] that if a 
fall in the reactor's magnetic flux occurs under the action of an external voltage (with switching rectifiers) then 
the relationship between the instantaneous values of the resulting magnetizing force and the reactor’s magnetic 
flux is given by the descending arm of the dynamic hysteresis loop of the reactor (EAD on Fig. 2). In this case, 
the difference between the three-phase and the single-phase amplifiers reduces to the difference in the laws of 
variation of the external voltage applied to the reactors. This circumstance does not affect the character of the 
relationship between the resulting field strength and the magnetic flux. 


It was shown above that, with the assumptions made, during the entire interval of change of the reactors’ 
magnetic flux, these changes occur under the action of the external voltage ugs [Cf., equation (20)] and therefore 
there can be no additional fall of the magnetic flux after the cessation of current flow in the load circuit. At 
the moment when the current flow in the load circuit ceases, the resulting reactor field strength becomes equal 
to Hcy and the reactor's magnetic flux becomes equal to $;. In the case considered, therefore, the relationship 
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between Hcp and 4, is defined by the descending arm of the reactor's dynamic hysteresis loop. 


The relationship between the magnetic flux increment, 46 = $, —@,, and the saturation angle can be found 
for the three modes of amplifier operation defined by (4), (5) and (6) by integrating the expressions for the voltages 
applied to the reactor (Cf., Appendix Il). These relationships have the following forms: 


for the mode defined by (4) 


A® = |! — cos (6, + =|: (22) 


for the mode defined by (5) 


1.517, +108 n\ 
AD = — "Was a (0,42); (23) 
for the mode defined by (6) 
1,50, - 108 
A®D = —_ (1 — cos 9,). (24) 
1° 


To construct the amplifier's contro] characteristics, it is necessary to give a number of values of ®,, to 
determine for them the corresponding values of load current by use of formulae (17), (18) and (19), and then to 
etic determine the magnitude of A® by formulae (22), (23) and (24). Then, knowing the configuration of the limit- 
ing dynamic hysteresis loop for the core material (or an approximation to it), one determines the corresponding 
values of H., and I, from the computed values of A@. 
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" Fig. 6. Characteristics of a magnetic amplifier Fig. 7. Relationships for a magnetic ampli- 

, with suppressed ac components in the control fier with suppressed ac components in the 

windings. control windings. 

As an example of this construction, Fig. 6 gives the control characteristics, curve 1, of the magnetic 

*: amplifier considered in the previous example. In this example, the load impedance is 400 ohms. The experi- 
mentally determined characteristic, given by curve 2, is shown for purposes of comparison. Figure 6,b gives the 

? limiting dynamic hysteresis loop for material 65NP which was used in the construction of the control character- 

fore istics. With lower values of source voltage the magnetic amplifier's control characteristic has an instability in 
the region of small load currents, the nature of this instability being the same as in the case of a single-phase 

1 amplifier [4]. To illustrate this phenomenon, curve 3 of Fig. 6,a gives the control characteristic of the same 

) amplifier with Up = 25 volts. 
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It is clear from Fig. 6,a that the experimentally determined contro] characteristics are somewhat lower than 
the computed ones, particularly in the upper portion of the operating segment. This is explained by the fact that, 
in the construction of the calculated characteristic, an idealization was assumed, amounting to the assumption 
that r +1, = 0. As was shown above, this assumption excludes the possibility of a fall in reactor magnetic flux 
after cessation of rectifier switching. In actual amplifiers such a fall always occurs. This may be shown by the 
same methods as are used for single-phase amplifiers [4]. As an example, we consider the second mode of oper- 
ation, defined by (5). For any magnetic amplifier in steady-state operation, the area bounded by the curve Wywdgf 

/d@ in the interval of increasing magnetic flux equals the area bounded by this curve in the interval of falling 
magnetic flux. It was shown above [Formula (11)] that in the interval of reactor magnetic flux variation in the 
mode defined by (5), the voltage 

ta 


Uqgg= 1.5ug>W, — + 


On Fig. 7, hatched area ABCD corresponds to the growth of magnetic flux #,.,. If we setr+ 1, = 0, the 
fall of magnetic flux oa) ac Will correspond to area EGHK. Starting at the moment 9s, and until the end of the 
interval of reactor ag conductance, u,,:= 0. Magnetic flux Paw cannot vary. If, however, we assume that r + 
r, = 0, then magnetic flux pf ac an fall in this interval, since ugg* = ug, (r + rr)/(R+1r+t;). This voltage 
must balance the emf of reactor aa, e,,, = —W, wd Pace AO. The additional fall of magnetic flux ag is 
denoted by the hatchedarea in the interval KL. 


Thus, if it is assumed that in the entire interval from point E to 9,,, the fall of magnetic flux "ae occurs 
under the action of an external voltage (as was done for the idealized case) then, since 9s,, = 9,,,,+ 1, it is 
found that magnetic flux aq falls by a larger amount than it increases (due to the hatched area in the interval 
KL), which is impossible. Consequently, one can only assume that magnetic flux f ao, falls in the interval E 0,4, 
in accordance with a different law than it increases in the interval AOsaq. Indeed, at the point 0,, where the 
relationship Wywd Pau /4® = ugp becomes valid (current ia passes through zero), rectifier ac stops switching and 
magnetic flux p ea begins to fall under the action of the magnetizing force (mf) of the control winding along 
line MN. With this, obviously, area MHN equals the hashed area in the interval KL. 


The lowering of the reactor control characteristic obtained experimentally, as compared with the computed 
characteristic, is also explained by the presence of a relaxation drop in the reactor's magnetic flux. It is clear 
that with an increase in the ratio (r+ r,)/(R + r+ t,) the discrepancy between the experimental and the computed 
characteristics will increase. A quantitative accounting of the effect of the finite impedance r + r, on the course 
of the control characteristic is hampered by the fact that there has not yet been presented an analytic expression 
for the fall of magnetic flux during the time of action of an mf. 


3. Voltage on the Rectifiers 








For a correct choice of rectifiers, it is necessary to know the magnitude of the back voltage on them in the 
two magnetic amplifier circuits considered above. 


We now analyze the operation of the rectifiers in a magnetic amplifier with free ac components in the 
control windings. Let us consider the operation of the rectifiers of phase a. 


In the interval of reactor magnetic flux variation, 9,,.,—"/3<® < 9a, (Fig. 3), the following equation 
holds 


1Pag 
dj ~° 





Ugo — ial? = — Ww (25) 


Simultaneously, equation (2) is valid. 


Rectifier aa switches and rectifier ag does not switch, since the change in magnetic flux Gp is connected 
with the reactors’ contro] circuit. There is a negligible back voltage on rectifier ag due to the actual finite 
impedance of the control winding loop. 


Pag 


WLo— <Ww n° 


0 





dG..e 
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In the saturation interval, On acr* 0< 9 sack 2n/3, rectifier aca continues to be switched but, on rectifier 
ag , there appears the back voltage 


Uap = la(r +r), (26) 


whose amplitude is easily computed if the magnitude and form of current i, are known. If there is residual 
inductance I, in the reactors, the magnitude of ya, takes the form: 





U rap = la V (r +r, + (oL,)?. (2°) 


The voltage of rectifier ag is shown on Fig. 3,e. On Fig. 8,a.there are given the oscillograms of phase cur- 
rent and voltage on one of the rectifiers, which verifies what has been said. 


Mey A 


me int 














Fig. 8. 


In the other modes of magnetic amplifier operation [the modes defined by (4) and (6)], everything so far 
stated remains in force. Thus, the maximum back voltage on the rectifiers in the case considered occurs with 
the maximum current in the load circuit. 


As an example, we give a comparison of the computed and the measured voltages on the rectifier in the 
circuit whose parameters were given in Section 1. 


The measured circuit impedances were r = 22 ohms, r, = 2.5 ohms (for a current of 180 milliamperes); 


the load current for I, = 0 was IL max = 180 milliamperes. According to the computations, Ur max = 4.4 volts 
while, according to the measurements, Ur max = 4.5 volts. 


Thus, in the given case, the maximum amplitude of the back voltage on the rectifiers was about 0.106 
U 
p max: 


Rectifier operation is different in a magnetic amplifier with suppressed ac components in the control 
windings. In the interval of increasing magnetic flux aq (Fig. 7), rectifier aa is switched. It remains switched 


in the interval of reactor aa conductance, and also in the interval E 6, where a fall of magnetic flux v1) aq Occurs 
and the following relationship is valid 


Poa 
Wy, 1) 7 > Uaa’- (28) 


Switching of rectifier aa ceases for 9 = 6). In the interval 6, <9 < 9,,, the following back voltage is 
applied to it 


d d 
Uae = Yaa — Ww Pet = tay — Wi, oe (29) 
In the interval of reactor ag conductance, rectifier aa again begins to be switched, since voltage u,,+ for 
® = Osa, sharply falls from wag to U go(r + £,)/(R + 1 + t,) (Fig. 7) and relationship (28) again becomes valid. The 
oscillograms of the phase current and voltage on the rectifier are given in Fig. 8,b, With small values of control 
current, the voltage on the rectifier can have a very different character. 


In this case, the fall of magnetic flux aq: under the action of the control winding's mf, either is concluded 
in the interval 0, < @ <9,,, and then, at the end of this interval, Tl /d® = 0, or proceeds so slowly that the 
following inequality is valid at the end of this interval 
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- + Tr IPoa 
Van Rr +r > Wie | (30) 









Then, in the interval of reactor ag conductance (0 sag on Fig. 7), a back voltage equal to 

































r+r df 
Uraa = Uaog rt —_ Wi, @® a (31) 
’ r 





remains applied to rectifier aa (Cf., the oscillogram of Fig. 8,c). 


As the ratio (r+ r,)/(R + r+ r,) increases, such a character of voltage on the rectifiers remains for larger 
magnitudes of control current. It should be mentioned that, in an amplifier with free ac components in the control 
windings, the magnitude of the maximum amplitude of back voltage on the rectifiers is approximately propor- 
tional to the ratio (r+ r,)/(R + r + t,) while, in the given case, its magnitude is virtually independent of this 
ratio (when the condition R >> r + 1, holds). This becomes understandable if one takes into account that here the 
quantity U; max is determined in accordance with (29), where none of the quantities in the right member depends on 
the ratio (r + r,)AR + r+ rp). 


The absence of any analytic expression for the fall of magnetic flux comes under the action of a constant 
mf does not permit the computation of the magnitude of the back voltage on the rectifiers for various modes of 
amplifier operation. Experiments have shown that, in the mode of operation defined by (6), the back voltage on the 
rectifiers in amplifiers with suppressed ac components in the control windings attains a maximum for saturation 
angles of the order of 1/6. The amplitude of the back voltage on the rectifiers can reach a value of about 
0.7U5 max- Thus, one can conclude that the back voltage on the rectifiers in amplifiers with suppressed ac 
components in the control windings significantly exceeds that in amplifiers with free ac components in the control 
windings. 


By comparing the properties, resulting from our static mode analysis, of the two circuits considered of 
magnetic amplifiers with self-saturation working on three-phase loads, we can conclude that the circuit with free 
ac components in the control windings possesses a larger maximum gain and requires a smaller number of rectify- 
ing elements, due to the lower back voltage on the rectifiers. The circuit with suppressed ac components in the 
control windings possesses a better linearity of control characteristic. With the proper choice of the magnitude of 
the supply voltage, this characteristic can be considered to be virtually linear. 


Appendix I 





For deriving the functional relationship between control current and saturation angle, we use an idea sug- 
gested in [3] for a single-phase amplifier. 


We consider a pair of reactors for any of the phases (for example, phase a). If, in some half-period of 
supply voltage, the magnetic flux of reactor aa increases while the magnetic flux of reactor ag decreases, the 
instantaneous value of the magnetizing component of reactor ag 's load current is determined from the expression 











H p 
ioe = ee ( 1—+'a"}. (32) 
By taking (13) into account, we get 
, H, |, Pap 


The mean value of control current is expressed by the formula 
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& 
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= = | ied =— aH [o, — 4) — Ss | Fant |. (34) 
Oy "1 


Here, 9; is the angle at which the interval of reactor magnetic flux variation begins. 


We consider initially the mode of reactor operation defined by (5). For this mode, the magnetic flux Pag 
is governed by expression (12). By substituting (12) in (34), and taking into account that u, = U,;,sin 0 and that, 
in the given case, 0, = 6, — 1/3, we obtain 








0 
H,! 108 
ee kn :. (« —he { 50 tin 0) a0 | (35) 


After integration, and replacement of U,;p by a function of & in accordance with (1), expression (35) takes 
the form: 





Hi, l. bid & 4 ‘ 
We [5 -»(F +08 cos 6, — 0.7 sin 6, )]. (36) 


For the mode of operation defined by (6), magnetic flux @,, has an expression analogous to (12) with the 
sole difference that the lower limit of integration will be zero (0; = 0). Then, by substituting the expression for 
Pag in (34), we obtain 


6 t) 

Aol oy v : 10° | ; 

= aw, | %-o, \ (%.—We| 1.5 Uy, sin 0) a0]. ° (37) 
0 0 


By proceeding with (37) analogously to what we did in the previous case, we obtain, for the mode of 
reactor operation defined by (6), 


H, |, 
I, We [9, —v (1.73 sin 6, — 0.73 8,)}. (38) 

For the mode of operation defined by (4), the reactors’ magnetic flux varies, not in accordance with one 
law, as in the previous cases, but in accordance with three different laws, as a function of the value of the para~ 
meter (cf, Equations (20), (25) and (29) in[1]). Therefore, in the right member of (34), instead of 


9 
\ Pap d0 we can write the sum of the increments of magnetic flux 
"1 
Tt 
% a 1/2 6, 
SA, = \ P_p49 + \ Pap49 + \ Pap. 
"A o, — - 4 n/2 


The limits of integration correspond to the limits in which equations (20), (25) and (29) of [1] are valid. By 
substituting in (34) the total increment of magnetic flux, TAG, s , which is computed in an analogous fashion to 


the previous cases by using relationships (2), (5), (7), (20), (25) and (29) from [1] (the explicit expression for LA@ag 


is not given here for lack of space), we obtain the following expression for the control current 











































H, l, % . 
l= — We [(o,- “-\ —v(— 1,05 + 1,38 sin 0, — 0.4 cos 8,) | (39) 





Appendix II 





As an example, we give the derivation of the expression for the load current in the mode of reactor oper- 
ation defined by (5) (Fig. 3,a). In the interval 0, < 9 < 94. + /3, the current of phase a is defined by the 
formula 


ee (40) 


and, in the interval 0,,,,+ 1/3 <9 < 9,,,,+ 29/3, by the formula 


ta = a (41) 
The mean value of load current is found from the expression 
, % +5 +> 
I == izR — \ u,, a9 a \ u,, 40 => 
6, ® +7 
Tr a 
%+5 %®+5 (42) 


=zAl - ( V3U,,8in (0 —*) d0 + \ V5u,,sin(@—7) a0 |. 


j 


By carrying out the indicated computations, we get 


g — oem sin(,+ 5). 





The values of load current for the modes defined by (4) and (6) are computed analogously. 


Appendix III 





We have previously shown the identity of the processes in the load circuits of amplifiers with suppressed ac 
components in the control windings and the processes in amplifiers with free ac components in the control wind- 
ings. Therefore, the equations obtained in [1] can also be used in the given case with the absence, referred to 
above, of the factor 2 in the denominator of the expression for A@ We now determine the increment of magnetic 
flux of reactor aa. For the mode of operation defined by (4), this increment is determined as the sum of three 
increments: 


A®,, = AD, + AD, + AD... (43) 


since the growth of magnetic flux in the interval 1/6 <9 <9, occurs in accordance with three different laws. 
For #/6 < 8 < 0,—#/3, expression (20) is valid, for 9,-w/3 <® <#/2, expression (25) is valid and, for 
a/2 <@ < 6,, expression (29) is valid, all three of these expressions being from work [1]. 


By substituting the values of the corresponding voltages in these expressions, we obtain 








> 1.5U ,,-108 1.5 U,,-108 < 
AD, ho sin @ d6= OZ ERiew (0,— =). (45) 
aut 

6, 

vee V3" m’ 10° | A 5x ¥3U,,-10° 1 T (46) 

A®,,.= “> sin (0 — =) d§ -- — ee (5 +cos (0,+%)|. 

" 
2 


If we substitute (44), (45) and (46) in (43), we obtain, after simplifying the expression 


AD,4 = saa [1 — C06 (0, + =) |: (22") 


For the mode defined by (5), equation (11) is valid, if one takes into account that the upper limit will be 9,: 


6, 
15U,-10* | 1.577, .408 
a eS i on ae Lad 3° 
A®,, Wee \ sin 6 d6 Wie cos (0, +3 ). (23") 
6, -> 
For the mode defined by (6), equation (41) of [1] is valid: 
15u,-10° "* 1.50, -108 
AD. = West \ sin 6 d0 = — (1 — cos 6,). (24°) 
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THE USE OF STATISTICAL METHODS FOR DETERMINING THE 
CHARACTERISTICS OF OBJECTS 


SURVEY 


Yu. P. Leonov and L. N. Lipatov 
(Moscow) 


The various statistical methods for determining the character- 
istics of linear objects are considered. 


Recently, in connection with the necessity of increasing control accuracy, the problem of the accurate 
determination of objects’ differential equations has arisen. In particular, several authors believe [1] that lack of 
ability to determine the characteristics of objects with sufficient exactitude is one of the basic causes limiting 
the accuracy of control. The best-known methods of determining objects’ characteristics, involving the appli- 
cation to their inputs of signals of given forms (sinusoids, step functions, etc.) have essential limitations. On the 
one hand, they entail the necessity of varying the system's mode of operation. On the other hand, there is ordi- 
narily “noise” in objects (particularly in industrial objects) which either essentially limits the accuracy of measure~- 
ment or makes it completely impossible. Most recently there have been developed statistical methods for deter- 
mining the characteristics of objects, these methods being, to a significant degree, free of the disadvantages 
enumerated above. In this paper we consider various well-known methods of determining object characteristics 
by statistical means and evaluate them on a comparison basis. In section I we analyze the application of statis- 
tical methods for determining frequency characteristics. Although the use of statistical methods ordinarily 
requires a large volume of computation, they allow the frequency characteristics of systems to be computed in 
cases where the ordinary methods turn out to be completely inapplicable. In section II we consider methods of 
determining the equations of objects without disturbing their operating conditions, and we do this for systems with 
one input, with many inputs, with feedback and for systems which are nonlinear. 


I. The Use of Statistical Methods for Obtaining Objects’ Frequency Character- 








istics [3] 


1, Preliminary Remarks 








The ordinary method for determining amplitude-phase characteristics consists of the following. To the 
object whose stability is to be investigated a harmonic stimulus (D sin wt or D cos wt) is applied; then, if the 
system is strictly linear, after the natural oscillations of the system have been damped down, there remains an 
output signal of the form B sin (wt —-9). The input and output signals are written on one and the same diagram 
(Fig. 1). One then computes the relative amplitude, A = B/D, and the phase shift, @, i.e., the modulus and 
phase, respectively, of the investigated object's frequency characteristic for the given frequency w. 


This widely-known method has one essential disadvantage. If the object to be investigated contains an 
internal noise source, or if additional noise is applied to the object's input, measurement of the output signal 
becomes either completely impossible or can be carried out only for low frequencies. In this case, the output 
signal has the form shown in Fig. 2. The noise is particularly powerful in reducing measurement accuracy at 
high frequencies, when the amplitude of the output signal becomes comparable with the noise amplitude. The 
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frequency range for which satisfactory results are obtained can be widened by increasing the amplitude of the in- 
put signal. 


However, all actual objects possess “saturation” and the last measure can, therefore, only partially amelio- 
rate the situation. It has been established in practice that the highest frequency for which the amplitude of the 


output signal of a “noisy” object can be measured is obtained when the amplitude of the output signal equals the 
mean square value of the noise. 


To overcome the difficulties stated above, we consider statis- 
tical methods of determining frequency characteristics. One of 
Output them (the amplitude method) makes it possible to determine, when 
a le am there is noise in the object, the imaginary and real parts of the 
frequency characteristic. The second (the method of null phase) 
is based on phase measurements. These methods require the 
mechanization of the computational work. However, they are 


nput applicable in cases where the ordinary methods are completely 
signal inapplicable. 








6! Time 
—_ _— 













Fig. 1. 2. The Amplitude Method 





This method is based on the following relationships. Let there 
be two signals, f(t) = Dsin wt and f,(t)=B sin (wt+ &)+ n(t), 
where n(t) is a stationary random function (the object's “noise”) 
and M {n(t)} = 0, M being the symbol for mathematical expectation. 


It is then easily verified that 


Object's output signal Time ‘ 
R, (0) = lim Ah fs(t) foit)dt=+ 4. BDeos®, (1) 
; T+o T ; 2 
Sinusoidal input signal Correspondingly, 
Fig. 2. T 
R, ,(0)= lim a f(t) fe(t)dt = + BD sin 6, (2) 


Taco 7 
0 


where f,(t) = Dcos wt. 


On the basis of (1) and (2) we have 





Z Jf ” >» Ri» (0) 
A= pr V Ri.(0)+ #3.(0), 0 = are R50)" on 


where A = B/D is the value of the amplitude characteristic for frequency w and @ is the value of the phase cha- 
racteristic at the same frequency. 


We note that R, 2(0) and Rs 9(0) are the values of the cross-correlation functions Ry 2( T) and Rs 9(T) for T = 0. 
It is easily seen, moreover, that Ry (0) and Rs 90) are proportional, respectively, to the real and imaginary parts 
of the system's frequency characteristic vector at the frequency w. 


It is completely obvious that computing Ry (0) and Rs (0) exactly is impossible because of the impossibility 
of implementing the passage to the limit as T-> o. In actuality, computation is not done by formulae (1) and (2) 
but by the formulae 


T 


Ry. (0) = -\ fall) fa(e) at, (1a) 
0 
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T 
Ry g(0) =| fat) fale ae, (1b) 
0 


The quantities Ry (0) and Rs 2(0) are estimates of Ry 2(0) and Rs 2(0). These estimates are the more accurate, 
the longer the time of integration T. In particular, it is advantageous to take T = k 2n/u, where k is a sufficiently 
large positive integer. We note that R, , and Rs, depend on T. However, if the interval of integration is chosen 
_ sufficiently large, R, z and Rs do not, for all practical purposes, depend on T. This serves as a criterion for the 
correctness of the choice of T. 


It is essential to remark that the computations of integrals (1a) and (1b) are necessarily carried out for the 
function f,(t) for t > t, where ty is a time after which the natural oscillations of the system have virtually ceased. 
Since these computations can be carried out on an analog computer, it is convenient to compute the estimates 
of Ry 2(0) and Rs 2(0) from the formula 


t —(t—x) 
Ri. (0) = a i \ *<.. fx (2) fa (x) dx, (4) 


To 3 


where T, is the computer's (model's) time constant. If T, is significantly greater than the period, 2n/w of the 
input signal, then R, 9(0) is virtually independent of t fort > 5Tp. 

The basic scheme necessary for obtaining the frequency characteristic by the amplitude method is shown in 
Fig. 3. 


3. The Method of Null Phase 





This method is based on the following relationships. If there are two signals, ft) = Dycos(wt + a) and 
f(t) = Bsin(wt +6) + n(t), where n(t) is a stationary random function and M{n(t)}} = 0, it is then easily found that 


T 
i i 
Rg o(0) = lim J fat) fa(t) dt = = DsB sin (0 — a). (5) 
0 
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Fig. 3. 1) Sine and cosine generator, P f, 
2) object and 3) computer. 











Fig. 4. 1) Sine and cosine generator, 
We now assume, as we earlier did, that f,(t) is the 2) object and 3) computer. 
output signal of an object at whose input the signal f,(t) = 
= Dsin wt has been applied and at whose output there occurs the noise n(t). Then, if the phase, a, of signal ft) 
is made equal to 9 then, in accordance with (5), we obtain 


Ry. (0) = 0 (6) 


for 
a = 60, (7) 
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Thus, by knowing the established phase a for which Rg 2(0) = 0, we obtain the phase of the frequency 
characteristic for frequency w. In order to obtain the modulus of the frequency characteristic, we initially estab- 
lish a value of a for which Rg (0) = 0, and then we so vary it that the algebraic sum of the phases becomes equal 
to 90°, 


Then on the basis of (5), we have Rg2(0) = D,B/2, from whence 


Ry 2 (0) 


DD, (8) 





B 
Au =2 


The scheme for obtaining the frequency characteristic by the null phase method is shown in Fig. 4. A smooth 
variation of the phase a in the limits (0, ¢#/2) is attained by means of potentiometer P and switch S. The value 
of a as a function of the ratio of the potentiometer arms, Pp, can be obtained from the diagram of Fig. 5. The 
formulae for a have the forms: 


for S in position S; 





a= —arctan( >?) : (9) 


for S in position S, 





¢=— ++ are tan 7): 


It is easily seen that, as P}, varies, there is a corresponding change in the amplitude of the signal fgt) = 
= Dy coswt + a), where D; = mp,D. The correction factor mp for taking account of the variations in amplitude D,; 
has the form: 


m, =V 1 —2P,+2P> (10) 


The curve for m), is shown in Fig. 6. 
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as a / 
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Fig. 5. Fig. 6. 


Because of the variations in amplitude Dj, it is convenient, in computing the amplitude characteristic, to 
use the formula 


A=2 R,, (0) 
m, . 
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4. _A Practical Circuit for Obtaining Frequency Characteristics by Statistical 
Methods 





The quantities Ry (0), Rs2(0) and Ry2(0) can be easily computed by means of an electrical model (analog) 
for solving differential equations. The functional scheinatic of the analog for determining the magnitude of 
these quantities is shown in Fig. 7. The feedback paths around the amplifiers allow the time constants to be 
varied within wide limits. Ordinarily, one chooses T,; = Tz. 
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Fig. 7. Mis a multiplier, A is an operational amplifier, and 
r is the potentiometer ratio. 
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Fig. 10. Fig. 11. 


Figure 8 shows the typical traces obtained in the determination of Ry,9(0) and Rs 2(0) by the amplitude 
method. If the amplifiers in the circuit of Fig. 7 are used as integrators then the values of Ry 2(0) and Rg 9(0) equal 
the mean slopes of the line obtained in this case. Characteristic traces are shown in Fig. 9. For very low fre- 
quencies of the sinusoidal input signal, a trace such as shown in Fig. 10 is obtained. With this, it is necessary to 
produce the average lines as is shown on the diagram. A typical trace of the determination of the magnitude of 
Rg.9(0) by the null phase method is shown on Fig. 11. 


For the methods considered, it is necessary that the sinusoidal and cosinusoidal signals have one and the 
same frequency and amplitude. These signals can be obtained either by means of a generator for sin wt or by 
means of a circuit set up on the model of Fig. 12. It is certainly convenient to obtain these functions by means 
of a model (analog). This can be done by solving the equation 


oT +e nl. | (11) 
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on the model. In practice, in order to obtain generation of a model, it is necessary to set up an equation with a 
negative damping coefficient, of the form 


d*z dz 
ae — 2ko a +o%'z =0, k <1. (12) 


In order to compare the statistical methods with the ordinary ones, a second-order equation was set up on 
the model (Fig. 13). A noise generator at the object supplied random noise. At the object's output there was 
connected an amplifier with variable gain, so that it would be possible to carry out measurements for small values 
of amplitude at the output. The measurements were made by three methods: a) by the ordinary methods, b) by the 
amplitude method (statistical), c) by the null phase method (statistical). 


The amplitude and phase characteristics for these three methods are shown, respectively, in Figs. 14-16. 


When the frequency characteristic is obtained by the ordinary method, the conclusion can be drawn, based 
on Fig. 14, that the accuracy in determining the characteristic falls sharply when the output amplitude is made 
less than, or equal to, the mean square value of the noise. Measurement accuracy is good only up to a frequency 
of 4 oscillations a minute (0.07 cycles), for which the amplitude of the output signal equals the mean square value 
of the noise. The exact amplitude and phase characteristics were determined theoretically on the basis of a well- 
known equation, which was set up on the model. 


Figure 15 gives the exact curve and the experimental points of the characteristic as obtained by the ampli- 
tude method. Accurate measurement was possible up to a frequency of 95 oscillations per minute. Thus, the 
upper limit on the frequency for which measurement could easily be implemented was 20 times higher than this 
limit with the ordinary method. 


The most accurate was the null phase method, however, it requires somewhat more time for the measurements 
to be made. Figure 16 gives the theoretical curve and the experimental points for this case. Accurate measurement 
here was also carried out up to a frequency of 95 oscillations per minute. However, the accuracy of the null phase 
method is, in practice, higher than that of the amplitude method, since the null phase method is a compensating 
one. The superiority of the statistical method is particularly clearly seen on Fig. 17, where the input and output 
signals are given with the presence of noise, and without it. It is clear from Fig. 17 that the sinusoidal signal at 
the output is completely hidden in the noise. It is impossible, not only to determine its amplitude and phase, but 
even to establish whether or not there is a sinusoidal signal. On the right is shown the output function which was 
separated from the noise by statistical processing. 


The examples considered of determining an object's frequency characteristic in the presence of noise proves 
the effectiveness of employing statistical methods. The statistical methods possess the following advantages. 


1. It is possible to obtain the frequency characteristic of an object with an internal noise source and with a 
significant noise level in the output signal. 


2. The characteristic of a “noisy” object can be obtained with a sinusoidal signal of low amplitude at the 
input. This, in particular, is important for objects with large gains where, even for small input signals, saturation 
manifests itself. This is also important for the experimental investigation of linearized nonlinear systems where, 
in a small neighborhood of some value of the input signal, the system can be considered as linear. 


The disadvantages of these methods include the following: 
1) the large amount of time necessary for taking measurements; 
2) the use of an electronic model with two multipliers (these latter must have high accuracy); 


3) to obtain the characteristic by the statistical methods, it is necessary either to have the model directly 
at the object to be investigated or to have a trace of its input and output which can be used the necessary number 
of times (for example, on magnetic tape). 
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Il, Determination of Objects’ Characteristics While They are in Operation [4-9] 
1, Preliminary Remarks 











The statistical methods for determining the frequency characteristics of "noisy" objects are very convenient 
when it is possible to apply sinusoidal signals to the input of the object under investigation. However, it is fre- 
quently required to determine a system's equation in its normal operating conditions without disturbing any con- 
nections or applying signals of special forms. With this, there is noise in the system which does not permit the 
registering, at the system's output, of the signal which is the reaction to the input signal. Under these conditions, 
the methods given in the previous section can not be used. 


We will consider at this point methods which are useful for this case. As was stated in the introduction, we 


shall present formally only the ideas of these methods without giving any rigorous mathematical foundation for 
them. 


2. Determining the Weight Function of a Single Input System [4] 





We consider a linear object (Fig. 18) at whose input there is a signal y(t) and at whose output there is a 
signal x(t). In the object itself there is a noise source which engenders an equivalent noise signal at the output. 
It is assumed that the statistical characteristics of the noise are unknown. In particular, it is impossible to obtain 
a trace (registration) of a realization of the object's noise n(t). Moreover, it is not known what effect the object's 
noise mt) has on the signal x(t). However, it is known beforehand that the object is linear and, consequently, is 
completely characterized by its weight function k(t).* Under these conditions, it is necessary to determine the 
object's weight function k(t) on the basis of traces of the realization of random 
| nt) functions at several points of the system. To solve this problem, we can proceed as 
follows. We imagine a signal at the object's output in the form 





Wt) ty 











t t 
z(t) = | k(x) (¢— a) de + (ky (2) mt — 2) de, (13) 
Fig. 18. 5 } 


where k(t) is the system's weight function between input and output and k,(t) is the system's weight function 
between the point where the noise arises and the output. Formula (13) is valid for linear objects. If we take 
into account that the weight function of a stable object has the property that | k(t)| > 0 for t-+o, we can always 
find a sufficiently large value of T such that | k(t)| * 0 fort =T. By using this remark, we can write formula 
(13) in the form 


Tj 


T 
z(t) =| k(x) y(t — =) de + (ky (2) m(¢—a) de, (¢ >7). (14) 
0 0 


In the sequel we shall assume that signals y(t) and n(t) are stationary random functions, Then, if we con- 
sider the value of the function x(t) after the natural motions of the system have been virtually damped down, we 
can consider x(t) as also being a stationary random function. Formula (14) is obviously valid for every realization 
of the random functions considered. It follows from this that, when there is noise in the object, it is impossible 
to determine the function sought, k(t), by registering the signals at the object's input and output, since the function 
k,(t) and the realization of the random function n(t) are unknown. We now assume that we know the random 
function z(t) which is correlated with function y(t) [and, consequently, is also correlated with function x(t)] but 
is not correlated with the object's noise n(t). Then, by using (14), we can write the following equation for the 
correlation functions: 


? 
R,, (t} = \ k(x) Ry, (t— t) ae, (15) 
0 


* We assume here that the object is a system with constant parameters. However, the method is easily applied to 
linear systems with variable parameters. 
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Ti 
1 
where R,, (t) = lim FT. \ w(u)z(u—t)du is the cross-correlation function of x(t) and z(t) and R(t) is the 
T,-»©o ! Ps yz 


cross-correlation function of y(t) and z(t). With the conditions assumed, cross-correlation function Ryg(t) = 0 and 
the second integral in equation (14) reduces to zero. In particular, in many problems it can be assumed that n(t) 
is not correlated with the input function y(t). Then, we can set z = y and rewrite (15) in the form 


T 


R,,, (t) = \ k(t) Ry, (t — 2) dx, (15a) 
0 


T; 
| 
where Ay, (t) = lim = \ y(u)y(u—t)du is the auto-correlation function of process y(t). However, it is 
yy T1000 T; ; 


necessary to bear in mind that y(t) is not the only random function which permits relationship (15) to be obtained. 


Equation (15) gives the solution of the problem posed. By solving it, we can determine k(t), the function 


sought. 


Thus, for the solution in the general case it is necessary to have the traces of the realizations of three functions, 
z(t), y(t) and x(t) but, if it is assumed that y(t) and n(t) are uncorrelated, then only the traces of two functions are 
necessary since one can set y(t) = z(t). We note that the function x(t) cannot, obviously, play the part of z(t). 


In conclusion we can make several remarks as to the computation of weight function k(t). First of all, it 
is impossible to compute Ryy(t) and Ryylt) exactly in substituting them in equation (15a) because of the finite 
duration of the traces. Only estimates, Ryy(t) and Ryylt) of these functions can be obtained. The question of the 
accuracy of determining k(t) because of errors in determining the functions R,,(t) and Ryylt) is beyond the limits 
of this work. It is quite obvious, however, that the more accurately Ryy(t) an Rxy(t) are computed, the closer k(t) 
will be to the solution k(t). The errors in Ryy(t) and R,,(t) due to the finite duration of the traces of the realiza- 
tions of x(t) and y(t) can be easily determined by very well-known methods [2]. 


3. Determination of the Weight Function of a Multi~Input Linear System [4] 








The case frequently occurs when it is necessary to determine the weight function of a system with several 
inputs and several outputs. The method considered for determination of the weight function is easily extended 
to this case, 


A system with several outputs does not give rise to any additional difficulties 
since each output can be considered independently of the others. Consequently, a 
We system with n outputs can be considered,as n independent systems have the same 
ee inputs. Therefore, we can limit our consideration to systems with several inputs 
Wf and one output. Initially we consider a system with two inputs (Fig. 19). Such a 
Fig. 19. system can be characterized by two weight functions, k(t) and k(t). We will assume, 


as before, that k(t) = ky(t) = 0 for t = T, where T is the system's memory.* For a 
two-input system, we can obviously write the following expression for the output signal: 





ky (t) ys (t — 1) dt. (16) 


at | 


T 
z(t) = | k(t) y(t — 2) de + 
0 


From this we obtain the system of two integral equations: 
T T 
Ry (t) = \ k(t) Ryy (t — 2)de +f be (4) Ryy, (t — 2) de, (17a) 


0 9 


*The “memory” is the interval OT in which the weight function k(t) is not identically zero. 
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T T 
R(t) = \ k(x) Ry (t— sat +4 ky (x) Ry, (t— 2) dz. (17b) 
0 


0 


If Ryy y(t) = © then, to determine k(t) and k,(t) we obtain two separate equations analogous to equation (15). 
If, however, y and yy are correlated, we can then separate from signal y,(t) a component z(t) which is not cor- 
related with y(t). Indeed, we have the equation 


-+co 
rn (t) = 2(t)-+ | ka (t) y(t — 2) de. (18) 


co 


On the basis of (18) we obtain, for the correlation functions, 


4-00 
Ryy, (t) = Rye (t) 4 \ kat) Ryy (t — 2) at. (19) 


—oo 


It is clear from this last expression that, if there exists a solution kg*(t) of the integral equation 


+00 
Ry, ) = { ka (t) Ry, (t — +) dz, (20) 


—oo 


then, by setting kg(t) = kg*(t), we obtain, from (19), Ry,(t) = 0. The presentation of signal y,(t) in the form of 
(18), where z(t) is not correlated with y(t) can be used to reduce the solution of the system of integral equations 
(17a) and (17b) to the solution of three separated equations. In fact, on the basis of (18) we have 


+co -+-09o 
Ry, (t) = yy (t) + 5 J ky (71) ky (42) Ry, (t — t+ 4) dtd te = we 
oe 
=R,,(t)+ \ kel) Ry, (t+ 24, 
+00 
22 
Rais (t) = R,, (t)+ \ keg (*) Rox (¢ — t)dt. (22) 
If we substitute Ry, yot) from (20) and Ry, (4) from (21) in (17b), we get 
T +0o 
Ryx(t) =| | (1) a (4a) Ryy (8 + te — a) dtd t+ 
0 —co 
T T +0o (23) 
+ \ ky (™1) RL, (t — 1) dt + \ ky (%3) ks (Ta) RR, (8 4 te — Ta) Ot Te. 
0 0 —co 
By substituting (17a) in (22), we obtain 
+00 T 
Ryx (t) = Ry, (t)+ \* (1) ana [he (ta) Ryy (t+ 1 — t) + ‘ees 


+- ky (t2) Rus, (t + — T:)] dts. 
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By comparing (23) and (24) we obtain 


. 
R,, (t)= \ ky (1) R,, (t — t1) dy (25) 
0 


If we interchange the roles of y and y,;, we can obtain an analogous equation of k(t). For using (25), 
functions R,,(t) and Rz,({t) are obtained from (21) and (22). We note that equations (17a) and (17b) as well as 
the method for determining the weight functions k(t) and k,(t) remain in force if, in the system (Fig. 19), there 
is a source of noise which is not correlated with y(t) and y;(t). 


4. Determining the Characteristic of a System with Feedback [4, 5] 








Let there now be a complex multi-loop system (Fig. 20) consisting of many elements, each of which can 
contain a source of noise. As before, we assume that it is impossible to obtain directly the exhaustive character- 
istics of the noise arising in any element of the system. Assuming the elements to be linear, we must determine 
the system's weight function on the basis of the traces of realizations obtained at definite points of the system. 
As for an object which is not shunted by a feedback path, the question arises as to which points of the system must 
be tapped for traces of realizations and how many of such realizations are necessary to compute the weight func- 
tion. To solve this problem we use the same method as in the case considered in section 2, namely, we register 
the realizations of three signals: y(t) at the element's input, x(t) and some auxiliary signal z(t) at the element's 
output. With this, as before, signal z(t) must be correlated with y(t) 
but not with the noise n(t) arising within the element. If it is then 
necessary, for example, to determine the weight function of element 
3 (Fig. 20), integral equation (15) must be solved for the latter. How- 
ever, the essential difference between objects shunted by feedback 
paths and open-loop systems consists in the choice of the function 
2t). In particular, for an element shunted by feedback paths it is 
no longer possible to choose the input function y(t) as the function 
z(t). Due to the feedback, y(t) will necessarily be correlated with 
the noise arising in the element investigated and, consequently, it 
cannot be eliminated by Equation (14). For the scheme of Fig. 20 
it is possible, in many cases, to choose function 1 (t) or Z,({t) as z(t). 
In particular, the method to be presented is applicable to the deter- 
mination of the weight function of an object in a single-loop system 
(Fig. 21). As signal z(t) one can frequently use the signal at the 
wt) ait input of the system (Fig. 21), but not at the input of the link being 

investigated. 














We now consider, as applied to the single-loop system (Fig. 21), 
the case when it is impossible to find a function 2(t) which satisfies 
Fig. 21. the necessary requirements. It is then possible to determine the 
object's weight function by carrying out an experiment twice for dif- 
ferent controller parameters. 


By denoting the functions of the first experiment by the superscript 1 and those of the second by the super- 
script 2, we car easily obtain the following formula [6]: 


2*Gt (jo) — 165, (io) 


(26) 
D,, (j0)'GX, — Veg (iw)P’EL, (jo) 


® (jo) = 





where (jw) is the object's frequency characteristic and po and $c,({jw) are the controller's transfer functions 


for the first and second tunings respectively. The function xx Jw) is defined by the formula 


1G... (jo) = GF, (/o)'Gz, (io), 
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where 1G jw) is the spectral density function of the output signal for the first tuning of the controller. By using 
this method we can also easily obtain the spectral density function of the object's noise. 


Thus, by knowing the statistical characteristics of the signals x(t) for two trials and by knowing cq and @cq, 
we can determine the object's transfer function and the spectral density of the noise Gpy for the case of the system 
of Fig. 21. However, the method considered possesses two disadvantages. First, intervention in system operation is 
required, since it is necessary to vary the controller's parameters. Second, the duration of the experiment is signi- 
ficantly increased and, during this time, the object's characteristics may change. 


5. Determining the Best Linear Approximations of Nonlinear Systems [4, 5, 9] 





Up untill now we have only considered linear systems, which are completely characterized by their weight 
function k(t). We note that the methods presented in section II can easily be extended to the case of nonstationary 


input signals. In this case we obtain more complicated integral equations for determining the weight function k(t), 
which will be a function of the variables t and T: 


“ 
R,, (t, u) = \ k(t, t) Ry, (u, t—*) dr, 
0 


where 
Ry (ty, ts) = M {Y (ty) Y (ta)] 


and 
R yy (tas ta) = M [X (t)) ¥ (ta). 


The method of determining the weight function of linear systems can be applied for a constructed linear 
system, the output signal of which is close (in the sense of a minimum second error moment) to the signal at the 
output of the given nonlinear system. Indeed, let the output signal of the nonlinear system be x(t) and the input 


signal y(t). We can then determine the equivalent linear system for the condition that the second error moment, 


o*, be a minimum, namely, 


T 
t a 
ochesh eget - ) dx} 


(27) 
By applying the well-known methods of variational calculus to determine the h(t) which minimizes the 
functional o? after computation of the mathematical expectation, we obtain the integral equation 
T 
Ryx (t) = \ h(t) Ryy (t—~ 1) de, (28) 
0 


It is easily seen that the function h(t) which satisfied (28) minimizes the function o*. In certain problems, 
the approximation just described of nonlinear systems by means of linear systems close to the nonlinear (in the 
sense cited) is very convenient. In particular, many nonlinear systems admit linearization. This means that, for 
small disturbances in the neighborhood of some value of the signal, they behave as linear systems. Thus, the 
quality and stability of such systems can be determined on the basis of a linear model with the condition that 
the disturbances be small. For this, the random disturbances must be small. 


If the random disturbances in the system's operating processes are large, the equivalent linear system, ob- 
tained by processing such traces, is a statistically linearized system. This means that all the nonlinearities in the 
nonlinear system are replaced by equivalent linear elements whose coefficients depend on the mean square value 
of the noise. With this, for each value of noise dispersion in the input signal, the coefficients are chosen from 
the condition that the signal at the system's output best approximate, onthe average, the output signal of 
the linear system. 
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With a sinusoidal input signal, 2 statistically linearized system coincides with a linear system computed 
on the basis of a harmonic output signal. In this case, statistical linearization is equivalent to ignoring all the 
higher harmonics. Despite the convenience of representing a nonlinear system by an equivalent linear system, the 
latter does not completely reproduce all the properties of the nonlinear system and can therefore not replace it jn 
all cases. 


6. Several Remarks on the. Solution of the Basic Integral Equation [4, 5, 7, 8] 





The solution of the basis integral equation poses no difficulties. One method of solving it consists of 
transforming equation (15a) to a system of linear algebraic equations by replacing the integral by a sum. We 
then obtain 


> a,,k, = by (v=1,...,), (29) 
ust 
where k,, == k (t,,), b, = Ray (tye Gy = Ryy (ty — 1,,)- The system of linear algebraic equations can be solved 


either on a general-purpose computer or on a special-purpose machine. The method of transformin; an integral 
equation to a system of linear algebraic equations was used in work [4]. A second method of solving equation (15a) 
amounts to the finding of a continuous solution by the method of successive approximations. This method was 
investigated as applied to optimal systems in[7, 8]. This method is applicable, without change, to the finding of 
an approximate solution of equation (15a). To determine the successive approximations, as was shown in [7, 8], 
one must use the formula 


= 
Kents (t) =k, (t)—a@ (| Roy (t — t) k, (t) dt — Row it) | (n = 0,4, ...). (30) 
0 


The initial (null) approximation, k)(t), in formula (30) is chosen arbitrarily. The number a must be chosen 

from the condition that « <= 1/Amax where Anax is the largest eigenvalue of the kernel, the correlation function 

t): In practice, one can take a = 1/2 in the majority of problems. We note that the algorithm of (30) is a 
variant of the method of steepest descent. However, the computations connected with (30) are much simpler 
than those used in the algorithm for the method of steepest descent which has been applied until recently. This 
simplification results from the fact that, in (30), a is chosen as a = const whereas, in the usual formula for the 
method of steepest descent, it is necessary to compute a painfully at each step. When formula (30) is used, con- 
vergence is somewhat slower than with the formula for steepest descent, but this decrease in speed of convergence 
is very significant and is completely compensated for by the choice of a constant a. 


Finally, for solving equation (15a) formally, the method of imaging (transforms) may be used. 
By using the two-sided Laplace transform, we obtain, on the basis of (15a), 


Rey (P) == K (p) Ryy (— p*), (31) 


where RyyfP): K(p) and Ry?) are the Laplace transforms of, respectively, the functions Rxy(t), k(t) and Ryyft). 
From (31) we have 


Ry (P) 
ss — . 32 
K (p) Ryy (— p*) ( ) 
It is essential to note that, in (31), K(p) will be the transfer function of a physically realizable system. 
However, in the author's opinion, formula (32) is impossible to use for direct calculations. Indeed, in practical 
cases, there are no functions Rxy(p) and Ryy(—p*), but estimates® of these functions, Ry yp) and Ryy(-p*). More- 
over, the processes for which the corresponding estimates are computed are only asymptotically stationary. All 


this leads to the circumstance that, instead of formula (32), one uses the formula 


*The estimates Rx yp ) and R paN 4 the functions Rxy(P) and Ry") are called random functions close, in some 
sense or other, to are and 
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; Rey (p) 
K (p) = = (33) 
R yy (— P*) 


where the superscript bar denotes an estimate of the corresponding function. It is obvious that, from (33), one can 
obtain a K(p) which corresponds, in particular, to a physically unrealizable unstable system. In this inheres the 
difficulty in a direct usage of formula (33). This difficulty can be avoided if, for solving the problem, we use, 
not formula (33), but a formula from the theory of optimal systems. It is well-known [2] that in optimal systems 
theory one solves an equation of the same form as in (15). However, the essential difference between the solution 
of the problem for optimal] systems and the solution of equation (15) consists in the fact that, in the first case, the 
functions Ryyl t) and R yt) are given, and it is required to determine the weight function which lies in the class 
of physically poboonce | e systems, while in the second case it is known beforehand that the system is physically 
realizable and the function X(t) is obtained as the result of passing the random function Y(t) through the system. 
It is obvious that the method of solving the equation for the optimal system is also useful for finding the weight 
function which satisfied equation (15a), With this, the convenience of using the results which relate to optimal 
systems consists in this, that with realizations Rxy(t) and Ryyt) of the correlation functions R,,{(t) and Ryyft), A 
solution is obtained which corresponds to a stable and physically realizable system, The realizations R, (t) and 
Ryyft) are computed by very well-known methods on the basis of realizations of the processes X(t) and Y(t). 


SUMMARY 


The statistical methods cousidered in section II have great potentialities as compared with the ordinary 
methods for finding frequency characteristics. In particular, the most important distinguishing feature of these 
methods is the possibility they offer of determining the characteristic of a system without disturbing the normal 
conditions of its operation. In connection with this, such methods are also applicable for the investigation of 
living organisms. 


On the basis of our analysis of the statistical methods considered, we can formulate the following basic 
problems related to the estimation of system operators under normal operating conditions: 


1) estimating the operators of linear systems with distributed parameters; 


2) estimating the operators of linear systems shunted by feedback paths without disturbing normal system 
operating conditions; 


3) estimating the operators of nonlinear systems, in particular, determining the nonlinearities and defining 
their characteristics; 


4) defining the accuracy of the estimates of system operators obtained. 
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